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ABSTRACT 


Calibration  devices  suitable  for  providing  known  amounts  of  nitric  oxide 
(NO)  at  temperatures  ranging  from  300  K  to  2000  K  and  pressures  of  0.5  atm 
(50.7kPa)  to  2.0  atm  (203kPa)  are  described  with  their  design  considerations. 
Methods  for  confirming  nitric  oxide  concentrations  are  given.  The  spectro¬ 
scopic  theory  for  the  absorption  of  ultraviolet  radiation  in  the  y(0,0)  band 
of  nitric  oxide  is  reviewed.  Experimental  values  for  oscillator  strengths 
and  broadening  parameters  for  NO  with  various  collision  partners  are  provided. 
Experimental  results  confirming  the  adequacy  of  a  computer  spectral  model  and, 
hence,  the  calibration  are  presented  along  with  the  details  of  the  model. 
Finally,  the  results  of  an  empirical  calibration  of  an  infrared  gas  correlation 
spectrometer  are  given. 
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I.  INTRODUCTION 


Since  Johnston  (1971)  and  Crutzen  (1970,  1972)  Independently  suggested  that  the 
injection  of  nitric  oxide  (NO)  into  the  upper  atmosphere  could  significantly  diminish 
the  ozone  (O3)  concentration,  an  accurate  knowledge  of  the  amount  of  NO  emitted  by 
jet  aircraft  has  been  a  serious  concern  to  those  Involved  in  environmental  studies. 

This  concern  intensified  when  McGregor,  Seiber,  and  Few  (1972)  reported  that  NO  con¬ 
centration  measured  by  ultraviolet  resonant  spectroscopy  were  factors  of  1.5  to  5.0 
larger  than  those  measured  by  extractive  probe  sampling  with  subsequent  chemilumine¬ 
scent  analysis.  These  initial  measurements  were  made  on  a  YJ93-GE-3  engine  as  part 
of  the  Climatic  Impact  Assessment  Program  (CIAP)  which  was  one  of  four  studies  (CIAP, 
NAS,  COMASA,  COVOS  (see  References))  commissioned  to  determine  the  possible  environ¬ 
mental  consequences  of  high  altitude  aircraft  operation,  especially  supersonic  air¬ 
craft.  After  those  studies  were  initiated,  economic  factors  strongly  favored  the 
production  and  operation  of  subsonic  aircraft.  Nevertheless,  since  the  subsonic  air¬ 
craft  fleet  is  large  and  does  operate  as  high  as  the  lower  stratosphere,  interest  in 
the  causes  of  the  discrepancies  between  the  two  NO  measurement  methods  continued. 

Few,  Bryson,  McGregor,  and  Davis  (1975,  1976,  1977)  reported  a  second  set  of  measure¬ 
ments  on  an  experimental  jet  combustor  (AVCO-Lycoming)  where  the  spectroscopically 
determined  NO  concentrations  were  factors  of  3.5  to  6.0  higher  than  those  determined 
by  the  probe  method.  In  this  set  of  measurements,  optical  data  were  obtained  not  only 
across  the  exhaust  plume  but  also  in  the  sample  line  connecting  the  probe  with  the 
chemiluminescent  analyzer.  The  sample  line  optical  data  seemed  to  agree  with  the  chemi¬ 
luminescent  analyzer  data;  hence,  it  was  suggested  that  the  discrepancies  were  due  to 
phenomena  occurring  in  the  probe.  These  results  stimulated  a  third  set  of  measurements 
involving  ultraviolet  spectroscopy  (Few  et  al,  1976a,  1976b),  infrared  gas  correlation 
spectroscopy  (D.  Gryvnak,  1976a,  1976b)  and  probe  sampling  on  an  Allison  T-56  combustor. 
The  measured  ratios  of  the  ultraviolet  to  the  probe  values  typically  ranged  between 
1.5  and  1.9  depending  on  the  data  reduction  procedure.  The  ratios  of  the  infrared  to 
the  probe  values  varied  between  1.1  to  1.5  also  depending  on  the  method  of  data  reduc¬ 
tion.  In  addition  to  these  engine  and  combustor  data,  evidence  supportive  of  the  ac¬ 
curacy  of  the  ultraviolet  spectroscopic  method,  i.e.,  calibration  data  and  model  pre¬ 
dications,  was  presented  by  McGregor,  Few,  and  Litton  (1973);  Davis,  Few,  McGregor 
and  Glassman  (1976);  and  Davis,  McGregor,  and  Few  (1976).  Nevertheless,  it  was  still 
not  possible  to  make  a  judgment  on  the  relative  accuracy  of  the  spectroscopic  and 
probe  methods.  The  most  significant  reasons  for  this  were:  the  complexity  of  the 
spectroscopic  theory  and  computer  model  required  to  infer  concentration  from  optical 
transmission;  the  inadequate  treatment  of  probe  use;  and  the  incomplete  exhaust  tem¬ 
perature  and  pressure  data  which  are  necessary  for  a  valid  comparison  of  the  methods. 
Recently,  Oliver  et  al  (1977,  1978)  as  part  of  the  High  Altitude  Pollution  Program 
has  ranked  these  discrepancies  as  a  major  and  a  continuing  source  of  uncertainty  in 
atmospheric  model  predictions. 


the  systematic  errors  associated  with  both  the  optical  and  probe  sampling  techniques 
for  measuring  NO.  To  accomplish  this,  the  study  was  divided  into  three  parts.  The 
first  was  devoted  to  calibrating  the  ultraviolet  and  infrared  spectroscopic  methods. 
This  entailed  developing  procedures  which  could  provide  known  concentrations  of  NO 
over  a  wide  range  of  temperatures  and  pressures,  and  also  reviewing  and  correcting 
the  ultraviolet  spectroscopic  theory  used  in  the  engine  and  combustor  measurements 
cited  above.  The  second  part  of  this  study  was  focused  on  sample  extraction,  transfer, 
and  analysis  by  chemiluminescent  instrumentation.  The  sampling  methods  were  used  on 
three  successively  more  complicated  combustion  systems  starting  with  a  flat  flame 
burner  and  culminating  with  a  jet  combustor.  The  results  are  presented  in  TASK  II 
Report:  Probe  Methods  .  In  the  third  part  of  this  study,  optical  measurements  were 
made  on  the  same  three  combustion  systems  operated  at  the  same  conditions  used  for 
the  probe  measurements.  The  results  of  the  optical  and  probe  measurements  were  com¬ 
pared  and  are  given  in  TASK  III  Report:  Comparison  of  Optical  and  Probe  Methods. 

This  report,  i.e,  TASK  I,  covers  the  results  of  the  first  part  of  this  study. 

As  indicated  above,  the  optical  calibration  entailed  several  steps  which  will  be 
briefly  described  below. 

The  first  was  the  selection  and  development  of  devices  which  would  provide  known 
concentrations  of  NO  from  room  temperature  to  approximately  2000  K.  Known  concentra¬ 
tions  at  pressures  from  subatmosphere  to  two  atmospheres  were  also  required  but  not 
over  the  above  temperature  range.  The  selection  of  the  devices  was  governed  by  the  fol 
lowing  criteria.  NO  decomposition  through  either  hetero-  or  homogeneous  processess  was 
to  be  minimal  ( £  10%)  relative  to  the  discrepancies  cited  above.  Temperature  and  con¬ 
centration  distributions  were  to  be  known  along  the  optical  path  and  stable  during 
' '  measurement.  The  optical  path  and  the  concentrations  were  to  be  such  that  unam¬ 
biguous  optical  data  would  be  obtained. 

The  second  part  of  the  calibration  process  started  with  the  translation  of  a 
spectroscopic  computer  model  developed  and  used  at  Arnold  Research  Organization 
(ARO)  by  the  authors  cited  above  (McGregor,  Davis,  Few.  Glassman) .  This  model  was 
generously  provided  by  J.  D.  Few  and  H.  N.  Glassman.  After  encountering  serious 
deficiencies  in  that  model  (see  Appendix  B) ,  a  new  model  was  developed  which  was  used 
to  process  both  resonant  line  and  continuum  source,  ultravilet,  spectroscopic  data.  In 
addition,  broadening  parameters  and  oscillator  strengths  were  experimentally  determined 

The  third  part  of  the  calibration  process  was  the  actual  verification  of  the  new 
model  over  the  temperature  range  given  above  using  a  hollow  cathode  resonant  lamp. 
Measurements  were  also  performed  by  ARO  personnel  with  a  capillary  discharge  lamp. 

These  results  are  described  in  an  internal  report  by  Few,  Lowry,  McGregor  and  Keefer 
(1979).  A  review  of  that  report  is  included  in  Section  IV. 

Finally,  an  empirical  calibration  of  an  infrared  gas  correlation  spectrometer 
was  performed.  The  results  of  that  calibration  are  presented  by  D.  Gryvnak  of  Ford 
Aerospace  in  Appendix  A. 
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II.  CALIBRATION  DEVICES 


A.  Design  Considerations 

II. A. 1  NO  Decomposition,  Equilibrium  and  Kinetic  Considerations 

The  complexity  of  the  problem  of  providing  known  concentrations  of  nitric 
oxide  at  elevated  temperatures  and  given  pressures  may  be  demonstrated  by  the 
paradoxical  predictions  of  equilibrium  theory  for  the  reaction 

2 NO  £  N2  +  02  (!) 

Assuming  infinite  time  is  available  and  an  initial  mixture  of  NO  in  nitrogen, 
equilibrium  values  of  NO  may  be  calculated  using  mass  balance  and  the  equilib¬ 
rium  constant  Kpet*  =  p^  x  pQ  /p^fjo  =  5.06  x  10-^  exp (+21 , 700/T)  (determined 
from  JANAF  tables).  In^Figurl  1,  where  ratios  of  equilibrium  to  initial  values 
are  plotted,  it  is  shown  that  as  temperature  decreases  from  2000K,  the  amount 
of  decomposition  becomes  greater.  Equilibrium  data,  in  fact,  predict  that  NO 
in  N2  mixtures  successfully  used  in  calibrating  analytical  instruments  should 
not  be  useful.  Clearly  equilibrium  information,  by  itself,  is  inadequate  since 
kinetics  play  the  dominant  role.  Indeed,  the  problem  of  decomposition  is  most 
severe  at  elevated  temperatures. 

Nevertheless,  Figure  1  suggests  that  calibration  procedures  and  apparatus 
must  be  selected  carefully  if  problems  caused  by  equilibrium  conditions  are  to  be 
avoided.  Fortunately,  the  rate  of  destruction  of  NO  is  limited  by  chemical 
kinetics,  which  in  turn  is  dependent  on  local  temperature,  pressure,  and  the 
presence  of  other  reactive  species  and/or  surfaces.  In  fact,  at  room  tempera¬ 
ture  the  homogeneous  reaction  rate  is  so  slow  that  pure  NO  can  be  stored  under 
pressure  for  indefinite  periods  with  negligible  conversion  to  N20,  N02  and  02 . 
Alternatively,  much  NO  can  be  lost  in  short  periods  of  time  in  the  presence 
of  a  catalyst  or  a  reactive  surface  such  as  copper  or  even  other  gaseous 
species  such  as  oxygen  where  NO  may  be  lost  via  the  reaction. 

NO  +  NO  +  02  -*■  2N02  (2) 

Furthermore,  it  is  clear  that  only  through  careful  consideration  of  the  rate 
of  chemical  kinetics  (both  homo-  and  heterogeneous)  can  wide  variations  of 
temperature  and  concentration  conditions  be  provided  for  NO  calibration.  A 
literature  review  of  both  the  homogeneous  and  heterogeneous  reactions  (e.g., 

Wise  and  Freeh  (1952),  Yuan  et  al. (1959),  McCullough  (1975X  Kaufman  and  Kelso 
(1955),  Winter  (1971),  Shelef  and  Rummer  (1969),  and  Fraser  and  Daniels  (1958)), 
suggests  three  temperature  regimes  for  the  decomposition  of  NO  at  atmospheric 
pressure:  (1)  below  1000K,  heterogeneous  kinetics  dominate,  (2)  above  1400K 

homogeneous  reactions  dominate,  and  (3)  between  1000  and  1400K  a  transition 
region  exists.  The  size,  i.e.,  temperature  window  of  the  transition  region 
for  a  given  experiment  depends,  of  course,  on  the  reactivity  and  surface  area 
of  exposed  walls,  the  concentration  of  NO,  and  the  identity  of  any  third  body 
(diluent)  if  present. 
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Of  the  several  reports  on  the  homogeneous  (thermal)  decomposition  of  NO, 
perhaps  the  most  comprehensive  is  the  recent  work  by  McCullough  (McCullough 
(1975),  McCullough  et  al.  (1977)).  In  this  work,  mixtures  of  NO  in  Ar  (100 
ppm  -  5%)  were  passed  through  a  packed  ceramic  flow  reactor  at  temperatures 
between  1500  and  2100K  at  atmospheric  pressure.  In  addition  to  presenting 
a  detailed  chain  mechanism  for  the  decomposition  of  NO,  this  work  provides 
experimental  data  for  NO  thermal  decomposition  at  the  concentration  levels  of 
interest  in  the  present  study,  i.e.,  100-5000  ppm.  Assuming  that  the  effective 
dead  volume  in  McCullough's  reactor  is  2.75  cm^  (estimated  from  this  published 
data  and  using  McCullough's  decomposition  data),  the  time  for  5%  loss  of  the 
initial  NO  can  be  computed.  The  results  are  shown  in  Fig.  2  for  several 
temperatures  and  concentration  levels.  A  reduction  of  _<  10%  during  test  time 
is  considered  to  be  the  objective  of  the  calibration  system. 


II. A. 2.  High  Temperature  Calibration 


These  data  suggest  that  an  apparatus  capable  of  providing  a  calibrated 
source  of  nitric  oxide  up  to  2000K  must  be  able  to  heat  the  NO  within  several 
milliseconds  and  continually  replace  the  decomposing  NO  with  freshly  heated 
NO.  Otherwise,  the  stationary  nitric  oxide  would  undergo  thermal  decomposition 
during  the  period  of  calibration,  which  for  the  infrared  (IR)  system  nay 
be  5  to  10  minutes.  Even  for  the  ultraviolet  (UV)  measurements,  a  minute  or 
so  is  required.  Since  a  flow  system  appears  necessary  in  order  to  continually 
replace  the  heated  NO,  a  static  heated  cell  would  clearly  be  impractical. 
Several  other  calibration  devices  can  immediately  be  eliminated  for  high 
temperature  use.  For  example,  NO  can  be  premixed  with  an  inert  carrier  and 
then  heated  by  flowing  it  through  a  hot  furnace  or  heat  exchanger.  This 
technique,  unfortunately,  is  not  feasible  since  the  time  required  for  heating 
of  the  gas  will  be  on  the  order  of  seconds.  Several  methods  of  preheating  the 
carrier  gas  and  then  injecting  cold  NO  prior  to  the  optical  measurement  were 
also  considered,  e.g.,  a  plasma  torch.  However,  in  these  cases  characteristic 
times  of  diffusion,  i.e. ,  mixing,  were  found  to  be  limiting.  A  shock  tube  has 
the  advantage  of  elevating  the  gas  temperature  nearly  instantaneously  (micro¬ 
seconds)  but  insufficient  time  is  available  (only  several  hundred  milliseconds) 
for  calibration  before  a  reflected  or  rarefraction  wave  changes  the  gas 
temperature  and  pressure. 


The  most  viable  alternative  which  would  both  heat  the  NO  within  several  milli¬ 
seconds  and  continually  supply  freshly  heated  NO  is  a  lean  ^/Oj/Ar  flame  seeded 
with  NO.  Argon  rather  than  nitrogen  is  selected  as  the  diluent  since  it  would 
prevent  additional  formation  of  NO  via  the  Feldovich  (1947)  reactions. 


This  technique  has,  in  fact,  been  used  successfully  at  UTRC  (Seery,  et  al 
(1979))  in  l/10th  atmosphere  flames  using  molecular  beam  sampling  and  by  Kaskan 
and  Hughes  (1973).  This  previous  work  indicated  that  NO  would  not  be  lost  in 
the  reaction  zone  of  a  lean  H2/O2  flame.  To  investigate  the  conservation  of  NO 
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In  a  post  flame  region  where  optical  measurements  would  be  made,  decomposition 
reactions  and  their  rates  were  examined.  These  calculations  predict  that  NO 
decomposition  should  be  small  immediately  above  a  lean  H2/O2/ Ar  flame  at 
atmospheric  pressure.  Due  to  the  very  wide  flammability  limits  of  a  hydrogen/ 
oxygen  flame,  it  is  expected  that  a  flat  flame  burner  could  be  used  as  a 
calibration  device  over  the  temperature  range  1000-2000K.  In  addition, this 
device  is  particularly  suitable  since  in  Tasks  II  and  III  of  this  program,  a 
flat  flame  burner  using  methane  for  fuel  is  needed. 

II. A. 3.  Low  Temperature  Calibration 

Given  that  a  ^at  flame  can  be  employed  for  NO  calibration  from 

1000-2000K,  the  lower  temperature  range,  300-1000K  with  possibly  some  tempera¬ 
ture  overlap  remains  to  be  considered.  A  literature  review  indicates  that 
heterogeneous  reactions  in  the  absence  of  high  concentrations  of  reactive 
radicals  dominate  below  HOOK.  Unfortunately  agreement  on  overall  reaction 
rates  and  heterogeneous  mechanisms  is  generally  so  poor  that  Shelef  and  Rummer 
(1969)  remarked  in  their  summary  of  catalytic  reduction  of  NO, 

"It  is  truly  so  that  catalysis  is  an  art  rather  than  a 
science  and  although  those  working  in  catalysis  are 
impelled  to  offer  plausible  mechanisms  for  the  catalytic 
process  based  on  kinetic  and  other  data,  these  attempts 
often  fall  very  short  of  exact  mechanistic  descriptions 
of  surface  processes  involved.  We  can  indicate  broadly 
those  solid  surfaces  that  will  be  of  use  in  the  catalytic 
reactions  of  NO  but  it  will  require  many  more  experiments 
before  we  can  hope  to  formulate  the  theory  that  will 
enable  us  to  predict  in  advance  the  catalytic  behavior 
of  a  given  surface  for  these  important  reactions." 

Conclusions  generally  differ,  for  example,  with  regard  to  the  experimental 
reaction  order  for  which  values  of  zero  (Yuan,  et  al  (1959),  and  Fraser  and 
Daniels  (1958)),  one  (Winter  (1971)),  and  two  (Wise  and  Freeh,  (1952))  have 
been  obtained.  Here,  the  reaction  order,  n,  is  defined  by  the  ovetall  rate  law. 

dtN0-l  — ' k[NO]n  (3) 

dt 


where  k  is  the  effective  rate  constant  which  in  general  depends  on  temperature, 
reaction  order,  surface  material,  and  active  or  available  surface  sites.  The 
reaction  order,  in  turn,  may  depend  on  NO  concentration,  temperature,  surface 
material,  and  available  surface  sites.  For  example,  Wise  and  Freeh  (1952) 
show  that,  by  doubling  the  surface  to  volume  ratio,  the  overall  decomposition 
rate  doubles  at  900K.  [NO],  in  Eq.  3,  represents  the  concentration  of  nitric  oxide. 


To  reduce  reaction  rates,  it  is  clear  that  total  surface  area  should  be 
minimized  and  that  a  relatively  inactive  surface  should  be  used.  Porous 
surfaces  should  be  avoided  due  to  relatively  high  total  surface  areas.  One 
measurement  (Fraser  and  Daniels,  (1958)  using  nitrogen  absorption  and  tte 
Brunauer,  Emmett  anJ  Teller  (BET)  theory  (Brunauer,  et  al  (1938)),  indicates 
that  the  total  surface  area  of  typical  metal  oxides  are  approximately  thirty 
times  their  macroscopic  surface  areas.  Polished  quartz  or  glazed  ceramic, 
of  course,  would  provide  ideal  surfaces  due  to  their  smooth  surfaces  and 
well-known  low  catalytic  activity.  Considerable  experimental  data  are 
available  on  ceramic  materials  such  as  alumina  (AI2O3)  or  zirconia  (Zr02)I 
however,  the  surfaces  are  not  glazed  and  NO  pressures  are  typically  much 
higher  (approximately  10-500  torr)  than  the  pressures  of  interest  in  this  work 
(1000  ppm  of  NO  in  one  atmosphere  of  N2  is  equivalent  to  .76  torr).  Extra¬ 
polation  of  this  decomposition  data  to  the  desired  conditions  of  the  present 
work,  i.e.,  low  initial  NO  concentrations  and  small  surface  areas  must  be 
questionable,  especially  in  light  of  the  comments  by  Shelef  and  Kummer.  In 
fact,  direct  application  of  their  equations  (Fraser  and  Daniels  (1958),  Mori, 
at  al.  (1973),  and  Yuan,  et  al.  (1959)) to  the  proposed  conditions  would  predict 
loss  of  essentially  all  NO  within  a  matter  of  seconds. 

Wise  and  Freeh  (1951)  measured  the  decomposition  of  NO  in  a  quartz  tube 
at  high  NO  concentrations.  They  assumed  two  competing  mechanisms,  homogeneous 
and  heterogeneous,  for  which  they  determined  a  reaction  order  of  two  for  both 
mechanisms.  Their  overall  equation  may  be  written 

iifU-  --(k!  ♦  k2)  1«0J2  <4' 

where  k-^  *  1.4  x  10^  exp  (-10770/ T)  cm^-mole-l-sec-^ 

k2  ■  3.1  x  10^-5  exp  (-41270/T)  cm^-mole”^-sec  ^ 

and  where  [NO]  represents  the  NO  concentration  in  moles/cc.  kjand  k2  are  the 
hetero-  and  homogeneous  rate  constants,  respectively,  for  a  quartz  vessel  having 
a  surface  to  volume  ratio  of  0.49  cm”  .  These  equations  suggest  that  it  takes 
approximately  one  day  at  1200  K  and  54  years  at  800K  to  decompose  5%  of  an 
initial  NO  concentration  of  1000  ppm.  However,  one  must  be  very  hesitant  in 
accepting  these  conclusions  at  face  value.  For  example.  Wise  and  Freeh  per¬ 
formed  their  experiments  at  the  equivalent  of  250,000  to  750,000  ppm  NO  in  one 
atm  (101  kPa  of  N^.two  and  one-half  orders  of  magnitude  higher  than  the 
concentrations  of  interest  in  this  work.  In  addition,  these  authors  do  not 
report  sufficient  information  on  their  experimental  technique  or  sufficient 
experimental  data  to  have  confidence  in  their  results  (see  Appendix  C) . 

Finally,  no  other  experimental  results  on  NO  decomposition  over  quartz  were 
found  in  the  literature  to  confirm  or  contradict  these  results. 
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If  Wise  and  Freeh's  conclusions  are  accepted,  then  a  static  quartz  cell 
could  be  used  up  to  1100  K  only,  of  course,  if  extreme  care  is  taken  so  as  to 
avoid  exposure  to  any  other  materials  within  the  heated  cell  and  no  gaseous 
impurities  such  as  oxygen  are  present.  Since  the  previously  published  report 
is  in  question,  a  simple  alternative  might  be  to  construct  a  heated  quartz  cell 
and  flow  mixtures  of  NO  in  a  carrier  through  the  cell.  Although  this  approach 
was  considered,  it  was  decided  that  uncertainties  in  temperature  profiles 
(radial  and  axial)  and  difficulties  in  probing  together  make  this  approach 
undesirable. 

Several  other  systems  were  considered  where  the  carrier  gas  is  preheated 
and  then  the  NO  injected  in  a  mixing  zone  just  upstream  of  the  optical  measure¬ 
ment.  All  of  these  designs,  however,  were  subject  to  limitations  by  characteris¬ 
tic  diffusional  rates  and,  therefore,  could  not  be  recommended.  An  alternative 
system  was  conceived  in  which  the  nitric  oxide  and  carrier  gas  is  premixed  and 
then  heated  prior  to  the  optical  measurement  up  to  the  temperatures  of  interest. 
The  gas  is  heated  using  a  quartz  bed  heat  exchanger.  Quartz  is  used  to  mini¬ 
mize  catalytic  activity.  This  concept  appeared  to  have  several  advantages. 

First,  uniform  mixing  of  NO  with  the  carrier  would  not  be  a  problem.  Second, 
a  design  could  be  selected  to  be  compatible  with  the  flat  flame  burner  and, 
consequently,  much  of  their  support  facilities  would  be  identical  and  could 
be  used  for  either  configuration.  Furthermore,  the  flowing  gas  would  be  easily 
accessible  to  probe  measurements  of  concentration,  and  temperature.  This  cali¬ 
bration  device,  entitled  the  flowing  gas  heater  (FGH) ,  seemed  to  be  most  practical 
for  the  lower  temperature  range  (300-1000  K) .  Details  of  design  and  operating 
characteristics  are  described  in  a  following  section. 

At  room  temperature,  a  static  cell  should  also  be  used.  With  this  device, 
concentration  and  temperature  profiles  are  well-defined  and  a  variety  of  experi¬ 
ments  can  be  performed  including  broadening  measurements  using  several  different 
gases,  and  comparison  of  results  from  certified  gas  mixtures  and  from  mixtures 
made  on  site.  This  device  is  described  in  Section  III.A.l. 

In  summary,  the  approach  for  calibration  of  the  optical  measurement  for 
nitric  oxide  includes  the  use  of  three  devices:  1)  a  static  cell  at  room  tem¬ 
perature,  2)  a  flowing  gas  heater  for  the  temperature  range  300-1000  K,  and  3) 
a  flat  flame  burner  with  a  flame  useful  in  the  temperature  range 

1000-2000  K.  Items  two  and  three  were  constructed  with  similar  geometry  so 
that  either  could  be  installed  in  a  central  facility  which  allows  for  optical, 
thermocouple,  and  probe  access.  This  facility  is  provided  with  a  gas  handling 
system. 
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B.  Description  of  Devices  and  Performance 


II.B.l.  Central  Gas  Facility 

The  central  facility  includes  a  gas-handling  system  and  a  water-cooled  shroud 
to  contain  the  gas  flow.  The  former  consists  of  necessary  plumbing  that  blends 
various  mixtures  of  gases  for  feeding  both  the  flat  flames  and  the  flowing  gas 
heater.  Separate  blending  facilities  were  constructed  for  the  main  and  buffer 
gas  flows.  The  stainless  steel  shroud  has  mounts  for  thermocouple  and  gas 
sampling  probes  and  provides  for  optical  access. 

A  diagram  of  the  gas  handling  system  is  depicted  in  Fig.  3.  The  gases  used 
in  these  experiments  were  argon  (Ar) ,  nitrogen  (N2) ,  oxygen  (O2) ,  methane  (CH4), 
hydrogen  (H2) ,  carbon  dioxide  (CO2) ,  nitric  oxide  (NO),  and  two  mixtures  of 
nitric  oxide  and  argon.  The  mixtures  were  approximately  10  percent  NO,  90  per¬ 
cent  Ar,  and  25  percent  NO,  75  percent  Ar.  For  each  mix  ratio,  several  bottles 
were  purchased  and,  for  each  bottle,  the  vendor  (Scientific  Gas  Co.)  specified 
the  ratio  and  certified  the  percentages  to  within  +  2  percent  of  their  stated 
values.  The  concentrations  were  verified  to  within  approximately  3  percent,  using 
both  mass  spectrometr ic  (MS)  and  chemiluminescent  analysis  (CA) .  In  the  latter 
case,  since  NO  concentrations  of  10  percent  and  25  percent  cannot  be  easily 
measured  using  a  standard  CA,  the  mix  ratios  were  verified  by  diluting  the  gas 
using  the  gas-handling  system  and  comparing  results  using  the  same  analysis 
techniques  with  gravimetrically-prepared  mixtures  of  NO  in  Ar  (Scott)  and  with 
mixtures  prepared  from  pure  NO  and  Ar. 

Flow  rates  of  gases  were  controlled  by  using  orifices,  metering  valves,  and 
0-150  psig  pressure  gauges.  The  gauges  are  8*5"  diameter  gauges  -from  Wallace 
and  Tiernan  and  are  accurate  to  within  +  0,15  psi.  The  orifices  were  all 
operated  at  critical  flow  and  ranged  in  diameter  from  0.0015"  to  .050".  The 
combination  of  orifices,  pressure  gauges,  and  metering  valves  allowed  the 
accurate  control  of  flows  ranging  from  8x10”^  to  3.5  grams/sec  of  .  The 
stability  of  the  flow  system  was  enhanced  by  using  two  stage  regulators  upstream 
of  the  metering  valves.  Although  pressures  and  temperatures  were  continuously 
monitored, only  occasional  adjustment  of  the  metering  valves  was  required. 

Orifices  less  than  .025”  size  were  all  jeweled  orifices  (Swiss  Jewel  Co.  and 
Brockton  Jewel  Co.),  mounted  in  brass  and  aluminum  discs  and  the  discs  were 
mounted  in  1/2-in.  stainless  steel  fittings.  Larger  orifices  were  machined  and 
placed  in  fittings.  Flow  calibrations  were  made  using  a  3-liter  wet  test  meter 
(CGA/Precision  Scientific)  for  flows  less  than  0.3  grams/sec  (N2)  .  A  calibrated 
volume  was  used  for  larger  flow  rates  and  for  hazardous  gases.  Typical  calibra¬ 
tion  curves  are  shown  in  Fig.  4.  Although  most  of  the  calibrations  were  performed 
using  nitrogen,  several  measurements  were  made  using  each  of  the  other  gases  in 
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SCHEMATIC  OF  GAS  HANDLING  SYSTEM 


TYPICAL  CALIBRATION  CURVES  FOR  CRITICAL  ORIFICES 


these  experiments.  Correction  factors  between  N2  and  other  gases  for  a  given 
orifice  were  within  3  percent  of  those  predicted  by  the  equation  (derived  from 
Shapiro,  1953) , 


where  the  subscripts  x  and  N2  refer  to  the  two  gases  of  interest;  m  is  the  mass 
flow  rate;  y  is  the  ratio  of  specific  heats;  and  MW  is  the  molecular  weight.  This 
equation  was  used  when  experimental  data  for  a  given  gas  and  orifice  were  unavail¬ 
able.  Discharge  coefficients  were  assumed  to  be  similar. 


Pressures  downstream  of  the  orifices  were  continually  monitored  to  ensure  that 


the  orifices,  indeed,  were  choked.  Pressure  ratios  (P 
always  3.0  or  greater  while  approximately  2.0  is  requiri 


/p  )  were 

tream  downstream 
to  maintain  choked  flow. 


Calibration  curves  were  repeatable  to  within  22  (usually  to  within  less  than 
1%)  for  measurements  using  either  the  wet  test  meter  or  the  calibrated  volume. 
Measurements  using  both  techniques  for  the  same  orifice  were  within  3%.  The  un¬ 
certainty  associated  with  blending  a  mixture  was  primarily  dependent  on  the 
uncertantles  in  the  orifice  calibrations  (-3%),  pressure  readings  (-0.5%)  and  the 
vendor  certified  and  cross-checked  Ar/NO  mixture  (<3Z) . 


For  the  nitric  oxide  feed  lines,  the  tubing,  fittings,  and  valves  were  con¬ 
structed  of  stainless  steel  to  minimize  corrosion  of  the  lines  and  loss  of  NO  to 
surface  reactions. 

Once  the  gases  mixed,  their  transit  time  to  the  test  rigs  was  less  than  one- 
half  second.  Oxygen  was  added  just  prior  to  the  flat  flame  burner  to  minimize  the 

total  volume  of  a  combustible  gas  mixture  and  to  reduce  the  effect  of  the  reaction. 

NOfNOfOj— —  2N02  (2) 

The  shroud  and  other  support  facilities  are  depicted  in  Fig.  5.  Its  purpose 
is  to  contain  the  gas  flow  and  to  prevent  outside  air  currents  from  disturbing 
the  flow.  It  is  essentially  a  rectangular  box  constructed  of  stainless  steel,  is 
open  at  both  ends,  and  has  internal  dimensions  of  22  cm  long,  14  cm  wide,  and  26  cm 
high.  Copper  tubing  is  silver  soldered  on  the  outside  for  water  cooling.  A  large 
(12.5  x  5  cm)  port  covered  by  quartz  is  installed  on  one  side  to  view  the  quartz 

bed  and  flames.  A  top  view  is  shown  in  Fig.  6. 
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STAJNLESS  STEEL  SHROUD  AND  ASSEMBLY 
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CALIBRATION  ASSEMBLY  (TOP  VIEW) 


Tubes  (3.4  cm  i.d.)  are  welded  on  opposite  ends  of  the  shroud  to  mount  the 
optical  windows  and  to  protect  the  optical  grade  quartz  or  salts  flats  from  the 
flame  environment.  Purge  ports  were  attached  to  these  tubes  to  flush  out  any 
residual  nitric  oxide  and  to  prevent  gaseous  diffusion  of  NO  into  these  colder 
regions.  The  gas  exit  from  these  tubes  was  1.3  cm  in  diameter. 

A  thermocouple  mount  and  traversing  mechanism  were  attached  to  the  outside 
wall  of  the  shroud  directly  above  one  of  the  tubes.  The  thermocouple  bead  could 
be  placed  at  the  same  height  as  the  optical  beam.  The  traversing  mechanism  allowed 
for  free  movement  (x  +  y  translation)  throughout  a  given  horizantal  plane. 

Accuracy  of  positioning  in  both  the  x  and  y  directions  was  approximately  +  1  mm 
relative  and  +  4  mm  absolute. 

The  probe  mount  and  traversing  mechanism  for  the  gas  sample  probes  were  mounted 
outside  the  shroud.  The  probe  mount  could  be  adjusted  vertically  approximately 
1  cm.  This  movement  was  sufficient  to  position  the  probe  tip  at  the  height  of 
the  optical  axis.  A  traversing  mechanism  was  purchased  which  allowed  for  x-y  move¬ 
ment  throughout  a  given  vertical  plane  and  with  a  relative  accuracy  of  +  0.1  mm. 
Probes  were  positioned  with  +  3  mm  absolute  accuracy. 

II. B. 2.  Flowing  Gas  Heater 

A  cut-away  view  of  the  flowing  gas  heater  is  shown  in  Fig.  7.  As  previously 
mentioned,  the  FGH  was  designed  to  be  a  low  temperature  calibration  device  with 
a  quartz  heat  exchanger.  Ceramic  heaters  (Electro-Applications) ,  which  are  capa¬ 
ble  of  reaching  1480  K,  heat  the  quartz  bed  primarily  by  radiative  heat  trans¬ 
fer.  The  quartz  bed  is  contained  by  a  long,  rectangular  quartz  tube  (17.4  x  9.2 
x  255  cm)  open  at  the  top  and  necked  down  to  a  2.5  cm  diameter  tube  at  the  bottom. 
Inside  are  cylindrical  quartz  tubes  ranging  in  diameter  from  5  to  15  mm  o.d.  with 
the  largest  at  the  bottom  and  with  their  sizes  gradually  decreasing  to  the  top. 

From  the  bottom  of  the  heat  exchanger,  two  chromel-alumel  thermocouples  in  alumina 
sleeves  were  inserted  to  monitor  internal  bed  temperatures.  Several  other  Cr-Al 
thermocouples  were  mounted  on  the  outside  of  the  cuartz  wall  using  a  ceramic  paste. 

The  buffer  zone  was  between  the  quartz  wall  and  the  internal  sides  of  the 
ceramic  heaters.  Its  equivalent  area  was  84  cra^.  A  buffer  zone  was  designed  to 
provide  a  region  that  would  have  similar  gas  temperature  as  the  main  flow  yet  be 
void  of  nitric  oxide.  The  gas  flowing  through  the  buffer  zone  was  always  the 
same  as  the  carrier  gas  in  the  main  stream.  The  buffer  region  was  sealed  on  the 
bottom  by  a  horizontal  quartz  plate.  Four  vertical  ceramic  heaters  were  cemented 
in  place  and  were  surrounded  by  a  stainless  steel  heat  shield  (reflector)  and 
approximately  10  cm  of  insulation. 

Due  to  its  bulk  and  to  material  constraints  the  flowing  gas  heater  could  not 
be  raised  or  lowered  once  in  position.  All  profiles  (temperature  and  concentration) 


I 
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were  obtained  at  a  fixed  location  above  the  gas  exit  of  the  quartz  heat 
exchanger.  This  distance  (6.5  cm)  was  at  the  height  of  the  optical  axis. 

II.B.2.a.  Probe  Measurements 


Several  probes  of  differing  design  were  made  in  order  to  access  various 
portions  of  the  optical  path  above  the  flowing  gas  heater.  These  uncooled  probes 
are  described  briefly  in  Table  II -A;  where  applicable,  figures  that  show  profiles 
of  NO  concentrations  also  indicate  the  probe(s)  used  to  obtain  the  data. 

The  extracted  gas  samples  were  analyzed  either  with  a  mass  spectrometer  or 
a  chemiluminescence  analyzer.  Details  of  these  analyses  systems  and  the  corres¬ 
ponding  sampling  trains  are  given  at  the  end  of  this  section. 

1 1 .  B .  2 .  b .  Temper  at  ur£  Mea  s£  r  emen  t_s__ 

Gas  temperatures  above  the  flowing  gas  heater  were  measured  using  a  chromel- 
alumel  thermocouple.  Wires  were  0.010  inches  in  diameter  while  the  diameter  of  the 
junction  is  approximately  0.015  inches,  head  wires  pass  through  an  alumina  rod 
approximately  30  cm  long  mounted  on  the  traversing  mechanism  at  the  side  wall  of 
the  shroud.  The  thermocouple  bead  was  supported  3.0  cm  below  the  ceramic  rod. 
Millivolt  outputs  were  measured  using  a  Data  Precision  Model  3500  DVM  with  a 
room  temperature  junction.  Radiation  corrections  were  calculated  and  found  to 
be  insignificant. 


1 1 .  B  .  2  .  c  .  Ex£e£ iment a. l_Resul 


Tests  were  made  to  investigate  NO  conservation  at  room  temperature  and  NO  pro¬ 
files  across  the  length  of  the  optical  axis.  A  typical  profile  is  shown  in  Fig.  8. 
The  loss  of  approximately  8  percent  NO  from  the  seed  value,  which  is  calculated 
from  relative  flow  rates,  is  apparently  due  to  a  small  percentage  of  intermixing  of 
the  buffer  and  main  flows.  For  these  profiles  and  for  all  data  presented  in  this 
report,  the  window  purge  was  on  at  a  flow  rate  of  .0109  moles/sec  (.0043  moles/ 
sec/cnr)  of  nitrogen.  For  all  calibration  experiments,  the  molar  flow  rate  of 
the  carrier  and  buffer  zones  were  held  constant  at  ,109  and  .0164  moles/sec 
(6.8  x  10-4  and  1.95  x  10-4  moles/sec/cm  )  respectively.  These  values  are  in¬ 
dependent  of  molecular  weight  of  the  carrier  gas  or  exit  temperature  of  the 
quartz  bed.  Under  these  conditions,  the  gas  residence  time  within  the  quartz 
heat  exchanger  is  approximately  three  seconds. 


Elevated  gas  temperatures  were  achieved  by  applying  power  to  the  ceramic 
heaters.  Periodically,  the  flow  rate  of  the  carrier  would  be  increased  to  the 
desired  flow  conditions  for  comparison  of  the  actual  gas  temperatures  to  the 
desired  test  conditions.  Two  hours  were  generally  required  to  reach  peak 
operating  temperatures  due  to  the  large  mass  and,  therefore,  large  thermal 
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TABLE  I I -A 


Description  of  Sampling  Probes  Used  Over 
The  Flowing  Gas  Heater  and  The  Flat  Flame  Burner 

Internal 


Probe 

Material 

Diameter 

(mm) 

Water  Cooled 

Or  if  ice 

inches 

Diameter 

(microns) 

A 

Quartz 

5 

Yes 

.025 

(635) 

Bd) 

Quartz 

5 

Y  es 

.006 

(150) 

C 

Quartz 

4 

No 

.012 

(305) 

D 

Quartz 

5 

No 

.035 

(890) 

E  <2) 

Stainless 

Steel 

1.2 

No 

.047 

(1200) 

p  (3) 

Stainless 

1.2 

No 

.047 

(1200) 

Steel 


(1)  Photographs  in  Figs.  12  and  13. 

(2)  Probe  inserted  horizontally  into  the  sampling  region  through  the 
optical  ports  which  had  their  windows  replaced  with  aluminum  plates 
with  a  small  hole  for  the  probe. 

(3)  Probe  inserted  vertically  down  the  center  of  the  shroud  but  had  a 
right  angle  bend  at  the  height  of  the  optical  axis. 
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NORMALIZED  CONCENTRATION  PROFILES  ALONG  OPTICAL  AXIS  OVER  FLOWING  GAS  HEATER 


a33S(zN/ON>/SV3w(zN/ON> 


7&-04  —  R4-7 


DISTANCE  FROM  CENTERLINE  (cm) 


inertia  of  the  heaters  and  surrounding  insulation.  Peak  gas  temperatures  that 
were  achieved  were  slightly  above  850  K.  Higher  temperatures  were  prevented 
primarily  by  local  hot  spots  in  the  ceramic  heaters  which  underwent  thermal 
runaway  before  other  portions  of  the  heaters  had  reached  their  peak  temperatures. 

The  cause  of  the  runaway  was  most  likely  due  to  nonuniformities  (bends  or 
connections)  in  the  alloy  heating  wires. 

Normalized  concentration  profiles  of  nitric  oxide  at  low  and  at  elevated 
temperatures  are  compared  in  Fig.  9.  Two  prominent  features  are  apparent.  First 
of  all,  nitric  oxide  is  conserved  within  the  quartz  heat  exchanger.  Although 
these  data  were  obtained  using  the  chemiluminescence  analysis  system,  similar 
results  were  obtained  with  the  mass  spectrometer  and  its  sampling  system.  Secondly, 
the  wings  (8.5  to  11  cm  from  the  centerline)  contain  noticeably  less  nitric  oxide 
at  elevated  temperature  than  at  room  temperature.  This  difference  appears  to  be 
real  since  profiles  obtained  at  many  different  dates  exhibit  this  phenomena.  It 
is  believed  that  the  difference  in  the  wings  is  due  primarily  to  a  reduction  in 
the  available  time  for  mixing  of  the  main  and  buffer  flows  at  elevated  temperatures. 
The  reduction  is  due  to  a  velocity  increase  that  is  proportional  to  temperature. 

When  heated  to  gas  temperatures  above  500  K,  it  was  found  that  power  input  to 
the  quartz  bed  was  insufficient  to  maintain  constant  temperatures.  This  problem 
existed  in  spite  of  the  fact  that  even  at  1/2  power  (2. A  kW) ,  the  ceramic  heaters 
provided  sufficient  power,  in  theory,  to  elevate  the  gas  from  room  temperature  to 
800  K.  Consequently  stable  temperature  profiles  above  500  K  could  only  be  main¬ 
tained  for  10-15  minutes.  Furthermore,  "temperature  waves"  could  be  observed 
within  the  quartz  bed  itself,  depending  on  the  immediate  history  of  heating  and/or 
cooling  cycles.  Consequently,  temperature  profiles  were  obtained  for  each  set  of 
experiments  performed  at  elevated  temperatures.  A  typical  profile  is  shown  in 
Fig.  10. 

II. B. 3.  Flat  Flame  Burner 

A  photograph  of  the  flat  flame  burner  is  shown  in  Fig.  11.  Its  surface  is 

sintered  copper  constructed  by  Thermit,  Inc.  Otherwise  the  burner  is  made  of 

sheet  copper  with  two  separate  zones:  the  main,  NO-seeded  zone  and  the  buffer 

zone.  Imbedded  in  the  5/8"  thick  sintered  surface  are  copper  coils  for  water 

cooling.  The  inside  dimensions  for  the  region  of  the  main  gas  flow  was  17.5  x 
o 

9.2  cm  or  161  cm  .  The  buffer  zone  is  slightly  smaller  than  for  the  flowing 
gas  heater  with  an  effective  area  of  76  cm  .  The  flow  systems  feeding  both  the 
buffer  and  main  zones  at  the  bottom  of  the  burner  were  described  earlier. 

II.B.3.a.  Probe  Measurements 

For  gas  sampling  over  the  flat  flame,  water  cooled  quartz  probes  were  used. 
Similar  designs  were  used  both  for  the  MS  analysis  and  for  detection  via  CA. 
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TYPICAL  TEMPERATURE  PROFILES  OVER  FLOWING  GAS  HEATER  AT  ELEVATED  TEMPERATURE 

AND  ALONG  OPTICAL  AXIS 


FLAT  FLAME  BURNER 


The  primary  difference  was  the  orifice  size.  To  feed  the  mass  spectrometer,  both 
low  sampling  pressure  and  very  little  sample  volume  are  required.  For  mass 
spectrometrlc  analysis,  probe  tips  were  small,  approximately  150  microns  in 
diameter  (several  tips  were  made  due  to  occasional  breakage  and/or  devitrification). 
A  typical  probe  and  tip  (probe  B  in  Table  II-A)  are  shown  in  Figs.  12  and  13, 
respectively.  (Although  Figure  13  presents  an  illusion,  the  water-in  tube  is 
indeed  placed  to  one  side  of  the  gas  sample  tube) .  The  CA  required  much  larger 
flow  rates  and  a  tip  opening  of  635  microns  in  diameter  was  used. 

The  large  diameter  of  the  water  cooled  probes  prevented  sampling  of  the 
flame  closer  than  1.5  cm  from  the  wall  of  the  shroud.  To  sample  the  gas  in  this 
region  close  to  the  wall  an  uncooled  quartz  probe  (C)  of  6  mm  o.d.  was  constructed. 
This  probe  was  only  used  on  a  low  temperature  flame  so  that  catalytic  activity  on 
the  hot  quartz  wall  and  divitrification  of  the  quartz  would  both  be  minimized. 

It  should  be  pointed  out  that  calculations  using  the  UTRC  probe  deck  were 
performed  for  the  quartz  probes.  Although  details  of  the  program  will  be  described 
in  the  Task  II  report,  it  is  worthwhile  mentioning  an  important  prediction.  For 
the  microprobes,  of  approximately  150  microns  in  diameter,  aerodynamic  quenching 
of  the  gas  sample  is  Impossible,  even  at  back  pressures  of  'v.  10  torr  used  in  this 
work.  The  limiting  factor  is  the  growth  of  the  boundary  layer  which  prevents  the 
flow  from  accelerating  supersonically.  [Although  plans  are  being  made  to  test  this 
prediction  on  an  atmospheric  pressure  flame,  limited  measurements  made  in  a  flame 
at  76  torr  (Seery,  unpublished,  1978)  suggest  that  the  mass  flow  rate  varies  with 
back  pressure  even  down  to  a  pressure  ratio  of  10].  This  pressure  ratio  should 
be  compared  with  the  ratio  of  2  to  1,  which  is  normally  believed  to  be  sufficient 
to  choke  the  flow.  These  conclusions  may  have  practical  implications  to  many 
flame  researchers  who  have  believed  previously  that  their  microprobes  were  aero- 
dynamically  quenched  even  with  pressure  ratios  as  low  as  2  or  3  to  1. 

II . B .  3 .  b .  Te^_erature  Measurements_ 

For  the  flat  flame  burner,  temperatures  were  monitored  with  a  buttwelded 
Ir/60%  Ir-40%  Rh  thermocouple  (.003"  wires)  coated  with  a  mixture  of  10%  beryllium 
oxide  and  90%  yttrium  oxide.  Kent  (1970),  has  found  that  a  mixture  of  these 
metal  oxides  provides  a  good  non-catalytic  coating  for  thermocouples  above 
2000  K.  The  bead  with  its  coating  had  a  total  diameter  of  .0035”  -  .0005" 
and  was  approximately  1.2  cm  away  from  either  lead/support  wire  (.010"  diameter). 

Photographs  of  the  thermocouple  assembly  and  a  close  up  of  the  yoke  are  shown 
in  Figs.  14  and  15.  Millivolt  output  was  measured  using  a  DVM  and  converted  to 
temperatures  uncorrected  for  radiation  using  published  tables  (Blackburn  and 
Caldwell,  1962). 

Estimates  of  radiation  corrections  were  calculated  by  equating  the  heat 
transfer  from  the  gas  stream  to  total  radiation 


WATER-COOLED  QUARTZ  MICROPROBE 


TIP  OF  QUARTZ  MICROPROBE 
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and  then  solving  for  the  gas  temperature  (Tg) .  The  wall  temperature  (Tw)  was 
assumed  to  be  zero  and  the  spectral  emissivity  (e)  at  6500  A  was  measured  to  be 
0.62  +  .1  using  a  pyrometer  which  agreed  well  with  Kent's  (private  communication) 
value  of  0.64.  o  is  the  Stephan-Boltzmann  constant  and  the  heat  transfer  co¬ 
efficient  (h)  was  calculated  from 


Nu 


(7) 


where  d  is  the  bead  diameter,  X,  the  thermal  conductivity  of  the  gas, and  where 
the  Nusselt  Number  (Nu)  was  determined  from  the  expression  (Hinze,  1959) 

Nu  -  .42Pr*2  +  .57Pr*33Re0,5  (8) 

and  calculated  values  of  the  Prandtl  (Pr)  and  Reynolds  (Re)  numbers.  Thermo¬ 
dynamic  and  transport  properties  of  each  mixture  of  gases  were  calculated  using 
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for  the  specific  heat,  the  Lindsay-Bromley  Formulation  for  thermal  conductivity. 


X 


m 


n  n 

Z  X±/[  1  +  Z  A  (X/X.)], 
i-1  j-1  J  J 

iH 


and  the  Wilke  Estimation  Method  for  viscosity 


(10) 


y 


m 


n  n 

Z  y./[l  +  Z  <fc  (X  /Xt)]. 
i-1  j-1  J  J 

J*1 


(ID 


The  subscript  m  represents  properties  for  the  mixture  of  gases  and  i  and  j  for 
individual  species,  where  n  represents  the  total  number  of  species.  The  expressions 
reported  in  Reid  and  Sherwood  (1966)  were  used  to  calculate  A.,  and  .  The 
molar  fractions  of  gas  (X)  obtained  from  the  flow  calculatlonsJwere  used  in  the 
above  determinations. 
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Properties  of  the  individual  gases  were  obtained  from  JANAF  tables  (C  ) , 

Hanley  (1973)  (y  and  X  for  Ar)  and  Svehla  (1962)  (y  and  X  for  all  other  gases) . 

Due  to  the  uncertainties  in  the  emissivity  calculation  of  transport  properties 
of  mixtures,  and  bead  diameter,  it  is  estimated  that  the  calculation  of  AT 
(=  T  -  T  )  is  accurate  to  within  20-25  percent. 

O  1  •  v*  • 

1 1 . B .  3 . c .  Exhaust^  Analyz er_ 

The  Scott  Model  119  Exhaust  Analyzer  used  in  this  present  investigation  pro¬ 
vides  for  the  simultaneous  analysis  of  CO,  CO2,  NO  or  N0o,  O2  and  total  hydrocarbons 
(THC) .  The  analyzer  is  an  integrated  system,  with  flow  controls  for  sample,  zero 
and  calibration  gases  conveniently  located  on  the  control  panel.  The  incoming  gas 
sample  passes  through  a  refrigeration  condenser  (^275K),  to  remove  residual 
water  vapor.  As  the  sample  passes  from  the  condenser,  it  is  filtered  to  remove 
particulate  matter.  The  system  is  comprised  of  five  different  analytical  instru¬ 
ments.  Beckman  Model  315B  Non-Dispersive  Infrared  (NDIR)  Analyzers  are  used  to 
measure  the  CO  and  CO2  concentrations  in  the  gas  sample.  Concentration  ranges 
available  on  the  CO  analyzer  were  from  0-200  ppm  to  0-15%  on  several  scales. 
Concentration  ranges  available  on  the  CO2  analyzer  were  0-4%  and  0-16%.  The 
accuracy  of  the  NDIR  analyzers  is  nominally  +  1%  of  full  scale.  A  Scott  Model  125 
Chemiluminescence  Analyzer  is  used  to  measure  the  NO  and  NO2  concentrations  in  the 
gas  sample.  Concentration  ranges  available  with  this  instrument  were  from  0-1  ppm 
to  0-10,000  ppm  on  several  scales,  with  a  nominal  +  1%  of  full  scale  accuracy.  The 
thermal  converter  used  in  the  chemiluminescent  analyzer  was  stainless  steel,  and 
was  operated  at  a  temperature  of  approximately  1000  K.  A  Scott  Model  150  Para¬ 
magnetic  Analyzer  is  used  to  measure  the  O2  concentration  in  the  gas  sample.  Con¬ 
centration  ranges  available  with  this  instrument  were  from  0-1%  to  0-25%  on 
several  scales,  with  a  nominal  accuracy  of  +  1%  of  full  scale.  A  Scott  Model  116 
Total  Hydrocarbon  Analyzer  is  used  to  measure  the  hydrocarbon  concentration  in 
the  gas  sample.  This  analyzer  utilizes  an  unheated  flame  ionization  detection 
system  to  provide  for  measurement  of  hydrocarbons  (as  carbon)  in  concentration 
ranges  from  0-1  ppm  to  0-10%,  with  a  nominal  accuracy  of  +  1%  of  full  scale. 

Output  signals  from  the  various  analyzers  are  displayed  on  chart  recorders  and 
a  digital  display. 

For  those  flows  in  which  Ar  was  the  bulk  or  carrier  gas,  the  chemiluminescent 
analyzer  was  calibrated  with  a  gravimetr ically  prepared  NO  in  Ar  standard.. 

This  instrument  can  also  be  calibrated  with  NO  in  N2  standards;  however  a  correc¬ 
tion  factor  of  1.20  must  be  applied  to  the  indicated  reading.  This  difference 
in  calibration  is  primarily  due  to  viscous  effects  in  the  instrument  sample 
capillary  and  varying  quenching  efficiencies  between  Ar  and  N2.  These  effects 
have  been  recently  analyzed  (Dodge,  et  al.  1979). 


A 


Gas  samples  were  transferred  to  the  SCOTT  analysis  instruments  through  a 
four  meter  (13  feet)  sampling  line  purchased  from  Technical  Heaters,  Inc. 

This  sampling  line  was  constructed  with  an  electrical  heater  including  thermo¬ 
couple  and  had  a  TFE  teflon  core  of  .48  cm  (3/16  in.)  internal  diameter.  For 
measurements  over  the  flowing  gas  heater,  samples  were  extracted  using  an  un¬ 
cooled  quartz  probe  (D)  with  an  orifice  diameter  of  890  microns  (0.035  in.). 
Typical  sampling  line  pressures  and  temperatures  were  500  torr  (^2/3  atm) 
and  300  K  respectively.  Under  these  conditions,  the  calculated  residence  time 
in  the  sampling  line  is  approximately  3/4  seconds  (mass  flow  is  assumed  to  be 
50  percent  of  choked  flow) .  For  sampling  from  the  flame,  probe  A  was  used  with 
an  orifice  of  635  microns  (0.025  in.).  Using  lower  pressures  and  elevated 
temperatures  of  380  torr  (1/2  atm)  and  110°C  in  the  sampling  line,  respectively, 
the  calculated  residence  time  is  similar,  approximately  one  second.  Although 
these  times  are  well  within  the  required  residence  times  required  by  the 
Federal  Register  (1976),  it  should  be  pointed  out  that  a  gas  sample  typically 
undergoes  much  longer  transit  times  as  it  passes  through  the  remainder  of  the 
sampling  train  prior  to  analysis.  For  example,  after  the  heated  sample 
line,  the  sample  travels  through  refrigerators,  metal  bellows  pumps,  filters, 
sample  bypass,  flow  valves,  all  connecting  lines  and  finally  the  chemilumenescence 
detector.  Altogether  these  components  increase  the  sample  transit  time  (from 
probe  tip  to  detector)  by  a  factor  of  15  to  20  from  what  is  calculated  in  the 
sample  line  alone.  The  increase  is  due  not  only  to  an  increase  in  total  line 
length  but  also  to  a  pressure  rise  (via  the  pumps)  and  the  sample  bypass,  both 
of  which  act  to  drastically  reduce  the  gas  velocity.  The  overall  effect  of  the 
sample  bypass  is  favorable,  since  it  provides  a  technique  to  continually  flush 
the  sample  lines  and  to  increase  the  flow  rate  and  velocity  through  the  first 
part  of  the  sampling  train.  Ideally,  the  bypass  should  be  located  immediately 
upstream  of  the  analytical  instruments.  These  comments  are  particularly  per¬ 
tinent  to  analysis  of  gas  samples  containing  high  concentrations  of  nitric  oxide 
due  to  the  reaction 


NO  +  NO  +  02  ■+  2N02  (2) 

which  is  second  order  in  NO  and  third  order  in  total  pressure.  In  the  case  of 
the  gas  correlation  measurements,  NO  concentrations  of  nearly  10,000  ppm  were 
required.  With  an  oxygen  rich  environment,  measurements  in  this  laboratory 
demonstrated  that  as  much  as  7  percent  of  the  NO  was  converted  to  N02  via 
Reaction  2.  Indeed,  kinetic  measurements  performed  in  this  sampling  system  pro¬ 
duced  an  evaluation  of  the  rate  constant  k2  which  was  within  25  percent  of  the 
accepted  value  (Hampsonand  Garvin,  1978)  and  were  consistent  with  the  fractional 
conversion  to  NO2  observed  during  the  IR  measurements.  Therefore,  any  NO2 
(N0X  minus  NO)  was  assumed  to  be  formed  in  the  sampling  line  (and  are  reported 
as  NO)  for  these  series  of  tests  where  nitric  oxiv.  concentrations  were 
exceptionally  high. 
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The  calibration  was  conducted  with  mixtures  of  NO  in  Ar  for  the  lower 
temperatures  (T  <  900  K)  and  by  seeding  NO  into  lean  H2/02/Ar  flame  at  the  higher 
temperatures  (900  K  <  T  <  2000  K) .  Standard  instrumentation  used  for  emission 
measurements  is  not  suitable  for  determining  H2  and  Ar.  In  addition,  since  the 
calibration  of  the  chemiluminescent  analyzer  is  dependent  on  the  bulk  or  carrier 
gas,  interpreting  NO  measurements  in  the  region  of  the  buffer  and  shroud  is  not 
s t  ra ight  f  orward . 

Mass  spectrometry  is  a  method  that  with  proper  sampling  technique  can  readily 
measure  homonuclear  molecules  (H2 ,  O2 ,  N2),  inert  gases,  and  nitric  oxide.  More¬ 
over,  proper  sampling  is  more  easily  achieved  since  the  mass  spectrometer  inlet 
and  sampling  line  can  be  operated  at  pressures  of  10  torr  (1.3  kPa) .  This  ensures 
a  rapid  reduction  of  the  sampled  gas  temperature  and  pressure  and  minimizes  sample 
transfer  time.  Since  only  small  mass  flows  are  required  to  make  a  measurement, 
the  probe  orifice  can  be  made  smaller  than  that  used  if  a  standard  analytical 
instrument  train  is  employed. 

The  instrument  used  in  these  measurements  is  a  one-meter,  time-of -flight  mass 
spectrometer  operated  at  a  source  pressure  of  5  x  10“^  torr  (6.6  x  10-i*  Pa). 
Residual  gas  pressure  is  typically  less  than  1  x  10“®  torr  (1.33  x  10-6  Pa).  The 
master  clock  for  the  instrument  is  crystal  controlled  at  10  kHz;  hence,  10^*  spectra 
per  second  are  obtained.  Through  the  use  of  ion  sampling,  the  high  frequency  real¬ 
time  output  is  converted  into  a  lower  frequency  suitable  for  a  data  acquisition 
system.  This  acquisition  system  was  a  Northern  Scientific  NS  575  signal  averager 
and  NS  408  F  tape  interface  with  a  Wang  Mod  7  digital  tape  transport.  For  the 
data  reported  here,  8  to  32  low  frequency  spectra  were  averaged  to  improve  the 
precision  of  the  data.  The  mass  spectrometer  inlet  was  maintained  at  8  torr 
(1.06  kPa) .  The  sample  line  was  heated  3/8"  stainless  steel  which  was  held  at 
a  temperature  no  greater  than  390  K.  The  terminal  pressure  of  the  sample  line 
when  capped  was  1.5  x  10-^  torr  (2.0  Pa).  Sample  transfer  time  was  less  than  0.5 
sec  for  He. 

The  procedure  used  to  reduce  the  data  depended  on  the  calibration  device. 

For  the  flowing  gas  heater,  the  procedure  consisted  of  measuring  the  Intensities 
at  masses  28  (N2) ,  30  (NO),  36  (Ar) ,  40  (Ar) .  ^(>Ar  isotope  was  used  as  a  monitor 

on  instrument  performance.  The  sensitivity  factors  for  these  constituents  were 
empirically  determined  from  gas  standards  and  from  mixes  of  gases  prepared  with 
the  critical  flow  gas  system.  These  sensitivity  factors  are  used  to  take  into 
account  viscous  effects  in  the  mass  spectrometer  gas  inlet,  ion  source  pumping 
speed,  Ionization  cross  section,  ion  gate  transmission,  and  ion  detector  mass 
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discrimination.  The  sensitivity  factors  used  relative  to  Ar,  i.e.,  S 


were  S: 


1.16  and  SAr 
NO 


Mole  fraction  of  NO  was  obtained  from 


XN0  -  S^.  1(30) /{S^.  1(28)  +  1(30)  +  1(40))  (12) 

For  the  flat  flame  burner,  the  procedure  is  made  complicated  by  the  difficulty 
in  making  an  accurate  H20  measurement.  The  standard  method  employed  in  emission 
measurements  is  to  dry  the  sample  in  a  low  temperature  trap.  Such  a  drying  pro¬ 
cedure  is  not  required  for  these  mass  spectral  measurements  since  the  total  molar 
flow  rate  of  Ar  is  known.  A  further  complication  is  introduced  by  the  mixing  of 
gas  from  the  three  zones  of  the  apparatus.  This  latter  problem  can  be  treated  by 
using  the  knowledge  that  the  N2  originates  only  from  the  window  purge  and  the  C0^ 
comes  only  from  the  buffer  zone  flame. 

The  composition  of  the  flow  at  any  point  can  be  described  by 


*M  +  +  %  =  1 

where  X  is  the  fraction  from  the  main  burner  flow,  Xg  is  the  fraction  from  the 
buffer  flow,  and  XN^  is  the  fraction  from  the  window  purge.  Using  N2  and  C02  as 
tracers,  the  local  intensity  ratios  of  N2  and  C02  to  Ar  are 


1(28) 

1(40) 


*N2 

^  MAr  +  S  BAr 


1(44)  =  1 

1(40)  sAl 


x  Brn 
B  co2 


XM  MAr  +  *b  BAr 


where  is  the  mole  fraction  of  Ar  in  the  main  burner  flow;  BAr  is  the  mole 
fraction  in  the  buffer  flow;  1(44)  is  the  CO2  intensity;  SA^  is  the  CO2  sensi¬ 
tivity;  and  BC09  is  the  mole  fraction  of  C02  produced  in  the^buffer  flame. 
and  are  known  from  the  calibrated  flow  system.  Bqq2  is  known  from  the 
equilibrium  value  of  C02  for  the  buffer  flame.  Solving  for  Xjj2  and  Xg  yields 


*N2  “ 


1(44) 


1(28) 


1(40)  \  +  1 
1(28) 


11-32 


1 


(17) 


sco2 

S5J  «28) 


The  results  of  equations  (16)  and  (17)  can  be  inserted  into  equation  (13)  which 
then  can  be  solved  for  X^. 

Given  these  expressions,  the  local  water  to  argon  expression  can  be 
calculated  from 


1(18) 

1(40) 


v  *b  +  V  x„ 

BAr  *B  +  MAr  XM 


(18) 


where  and  Q  are  known  from  the  equilibrium  value  of  H2O  for  the  buffer  and 
main  flakes,  respectively.  The  mole  fraction  at  any  point  for  a  given  ith  molecule 
is 


I(i)/I(40) 


(19) 


where  \ 

T  =  1  +  1X70)  |  S*  .1(44)  +  S*  1(28)  +  1(30)  +  S^.  (32)  +  I(18)j  (20) 

All  terms  are  measured  values  with  the  exception  of  the  I(18)/I(40)  ratio  which  is 
calculated.  Also,  for  the  buffer  region  and  a  small  region  on  either  side  of 
it,  corrections  are  made  for  the  contributions  of  C0(30)  isotope  to  the  nitric 
oxide  intensity  and  for  the  contributions  of  CO(28)  from  the  flame  and  ion 
fragmentation  of  CO2  to  the  nitrogen  intensity. 

I I . B . 3 . e .  Ex£e£iment£l_Results_ 

Initial  tests  of  several  flame  conditions  indicated  a  flow  instability 
several  centimeters  above  the  burner  surface.  To  dampen  these  flow  fluctuations 
a  three  level  tier  of  screens  was  constructed  and  hung  from  the  side  walls  of  the 
shroud.  The  lowest  screen  was  7  cm  above  the  optical  axis.  Slots  were  cut  at 
appropriate  positions  for  movement  of  the  probes  along  the  optical  axis.  Although 
not  mentioned  previously,  this  set  of  screens  was  used  in  all  measurements  above 
the  flat  flame  burner  and  flowing  gas  heater. 
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TABLE  II-B 


FLOW  CONDITIONS  FOR  OPTICAL  MEASUREMENTS 
OVER  FLAT  FLAME  BURNER 


Uncorrected 

Thermocouple 

Temperature 

(K) 

Mole 

Fraction 

B2 

°2 

Ar 

950 

0.073 

0.102 

0.825 

1220 

0.116 

0.114 

0.770 

1400 

0.169 

0.116 

0.715 

1600 

0.330 

0.170 

0.491 

Equivalence 
Rat  io 

Distance  of 
optical  axis 
above  burner 
(cm) 

Total 

Molar 

flow  rate 
(moles/sec) 

0.36 

2.0 

0.128 

0.51 

1.5 

0.134 

0.73 

1.5 

0.145 

0.92 

1.0 

0.194 
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After  selecting  several  flame  conditions  (listed  in  Table  II-B),  initial 
gas  sampling  was  made  with  probe  B  (Table  II-A)  and  analysis  was  made  with  the 
mass  spectrometer.  Typical  profiles  of  nitric  oxide  along  the  optical  axis  are 
shown  in  Fig.  16.  Flow  conditions  for  the  l^/C^/Ar  flames  are  given  in 
Table  II-B.  The  observed  scatter  is  characteristic  of  a  statistical  analysis 
for  the  operating  conditions  of  the  mass  spectrometer.  This  scatter  could  be 
reduced  by  additional  averaging  of  spectra.  It  is  apparent  that  nitric  oxide 
is  conserved  through  the  flame  front  at  least  for  the  low  temperature  flames. 
Similar  data  demonstrates  conservation  at  higher  temperatures. 

A  typical  thermocouple  profile  with  radiation  corrections  is  reproduced  in 
Fig.  17.  The  profile  shows  near  symmetry  around  the  centerline.  The  slight 
temperature  rise  near  the  edges  (around  8  cm  from  the  center)  is  probably  due  to 
a  separation  by  the  sintered  copper  and  the  copper  plate  between  the  main  and  the 
buffer  flows.  Resultant  gas  velocities  are  higher  which  push  the  flame  away  from 
the  burner  and  decrease  the  local  heat  loss  to  the  burner. 

At  approximately  8.5  centimeters  from  the  center  of  the  burner,  the  tempera¬ 
ture  drops  drastically.  This  fall-off  seems  surprising  considering  the  existence 
of  the  buffer  flame  out  to  10.2  centimeters  and  is  caused  by  the  strong  nitrogen 
purge  exiting  from  the  window  tubes.  The  effect  of  the  purge  at  the  selected 
flow  rates  is  quite  noticeable  on  both  the  temperature  and  NO  profiles.  Although 
this  perturbation  is  undesirable,  its  presence  does  not  affect  the  ability  to 
make  the  optical  measurements  since  probe  measurements  provide  sufficient  data  to 
reduce  results  using  multiple  zones.  Probe  measurements  indicated  that  lower  flow 
rates  were  ineffective  in  flushing  out  cold  NO  from  the  optical  arms.  Additional 
experimental  data  on  both  the  flowing  gas  heater  and  flat  flame  burner  are 
presented  in  the  following  chapter. 
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FIG.  16 


NO  HORIZONTAL  PROFILES  OVER  FLAT  FLAME  BURNER 
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Ill .  ULTRAVIOLET  ABSORPTION 


A.  Apparatus 

III.A.l.  Gas  Mixing  System  for  Static  Cell  Optical  Measurement 

The  majority  of  static  cell  calibration  data  were  obtained  by  mixing  gases 
using  the  system  shown  in  Fig.  18.  The  system  was  evacuated  below  25  x  10“^  torr 
(3.3  Pa.)  and  then  filled  with  a  calibration  gas  to  the  desired  pressure.  If  it 
was  desired  to  raise  the  pressure  of  the  diluent  gas  further  without  adding  more 
NO,  the  valve  to  the  static  cell  was  closed,  the  system  evacuated,  and  then  filled 
with  the  diluent  gas  by  slowly  raising  the  regulator  pressure  until  the  final 
desired  pressure  was  reached.  The  valve  to  the  static  cell  was  opened  so  that 
the  pressure  in  the  rest  of  the  system  dropped  and  then  quickly  recovered  via  the 
regulator  to  the  desired  final  pressure  and  the  valve  was  closed. 

Pressures  were  monitored  with  three  gauges.  A  thermocouple  gauge  was  used 
to  verify  the  vacuum  integrity  of  the  system.  A  Barocell  Datametrics  capacitance 
manometer  (model  570  A  sensor  and  model  1173  readout)  with  full  scale  ranges  from 
0-0.1  torr  (0  -  13.3  Pa.)  to  0  -  1000  torr  (0  -  133.3  kPa.)  was  used  along  with 
a  Wallace  and  Tiernan  0  -  150  psia  (0  -  1.034  MPa.)  pressure  gauge  to  determine 
pressures. 

For  the  broadening  and  oscillator  strength  study,  cylinders  of  NO  diluted  to 
about  2000  ppm  in  N2,  Ar,  CO2,  and  CH^  were  obtained  from  Scott  Environmental 
Technology  along  with  the  highest  purity  diluent  gases  commercially  available. 

The  specifications  for  these  gases  are  given  in  Table  III-A.  The  NO/  N2  con¬ 
centration  was  independently  verified;  but,  the  other  cylinders  were  accepted 
as  labeled  based  on  gravimetric  blending  by  the  vendor  and  vendor  analysis. 

III. A. 2.  Light  Sources 

Two  distinct  types  of  light  sources  were  used.  One  lamp  was  a  low  pressure 
hollow  cathode  lamp  which  produced  discrete  emission  lines  in  the  y  (0,0),  y  (1,1), 
and  y  (2,2)  bands  of  NO.  The  other  lamp  was  a  high  pressure  Xe  lamp  which  pro¬ 
duced  continuum  radiation  in  the  region  of  interest. 

I I I . A . 2 . a .  Nari^w^Line  Lamp 

The  hollow  cathode  lamp  was  operated  with  a  dc  discharge  in  flowing  air  and 
produced  emission  lines  from  principally  NO  molecules,  N£  molecules  and  ions,  and 

Ar  atoms.  The  spectral  lines  used  in  this  study  were  in  the  y  (0,0),  y  (1,1),  and 
y  (2,2)  bands  (A2I+  -  X2it)  of  NO. 
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TABLE  II I -A 


Specifications  for  Gases 


Gas  Blends  of  NO 

Diluent  Gas 

Ar 

C02 

cha 

N„ 


Pure  Gases 
Diluent  Gas 

n2 

C02 

CO 

ch4 

Ar 


Concentration 

Vendor  Analysis 

of  NO 

Accuracy 

(ppm) 

1980 

±1% 

2030 

±1% 

2000 

±1% 

2080 

±2% 

Purity 

Maximum  Oj 

99.998% 

3  ppm 

99.99% 

5  ppm 

99.99% 

5  ppm 

99.99% 

5  ppm 

99.999% 

2  ppm 
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The  design  of  the  water-cooled  hollow  cathode  lamp  is  shown  in  Fig.  19  and 
followed  that  of  Meinel  (1975),  who  generously  provided  a  drawing  of  his  lamp. 

The  operating  characteristics  are  similar  to  those  described  by  Meinel  (1975) . 

The  lamp  was  operated  on  flowing  air  at  a  pressure  of  2  torr  and  a  current  of 
25mA  (0.1%  electronic  stability  over  eight  hours).  The  ballast  resistors  were 
combined  for  a  net  resistance  of  1.67K  and  the  operating  voltage  was  about  600V. 

At  constant  pressure,  the  intensity  stability  of  the  lamp  was  excellent  (<.5% 
drift)  over  several  hours. 

The  possible  steps  leading  to  the  production  of  NO  in  an  excited  state  from 
air  in  the  discharge  are  discussed  by  Meinel  (1975).  It  is  clear  from  the  emission 
spectra  that  the  excited  NO  molecules  are  not  in  rotational  or  vibrational  thermal 
equilibrium.  The  population  distribution  of  the  2£(v'=0)  level  is  shown  in  Fig. 

20,  where  the  points  represent  resolved  spectral  lines.  The  lines  drawn  through 
the  points  were  used  as  an  input  to  the  computer  program  to  determine  the 
strengths  of  unresolved  lines,  as  explained  in  the  description  of  the  computer 
program. 

There  is  a  small  amount  of  light  in  the  NO  y  (0,0)  band  which  appears  to  be 
from  some  molecule(s)  other  than  NO.  The  NO  lines  may  be  readily  located,  but  there 
are  some  weaker  spectral  features  between  the  lines  which  have  not  been  identified, 
as  shown  in  Fig.  21.  These  extra  lines  are  not  absorbed  like  the  other  NO  lines 
as  shown  in  Fig.  22.  The  effect  that  this  "extra-light"  has  on  the  measured  f- 
values  is  discussed  in  the  section  describing  the  oscillator  strength  measurements. 

I II . A . 2 . b  Oantinuurn  Lamp 

The  lamp  used  to  obtain  a  continuum  output  was  a  1000  watt-  high  pressure  Xe 
arc  lamp  (Canrad-Hanovia  976C-0010)  mounted  in  an  Oriel  housing  and  using  an  Oriel 
power  supply.  The  lamp  exhibited  a  moderate  amount  of  intensity  drift  (3%  over 
10  min) ,  and  the  lamp  housing  had  to  be  vented  to  the  outside  because  of  the 
significant  quantity  of  ozone  produced. 

1 1 1 .  A .  2 .  c  Spectrometer  and_  As^o  c_i  at  ed_E  lec_t  £on  i£S 

A  1.5-m  focal  length  J-Y  spectrometer  in  a  temperature  controlled  box  with  a 
2400  g/mm  holographic  grating  (110  x  110  ram),  aperture  of  f/12,  and  Fastie  curved 
slits  was  employed  for  all  measurements.  Typical  slit  function,  full  width  at 
half  maximum  (FWHM),  was  observed  to  be  .0018  nm  with  7  pm  slit  settings  for  the 
226  nm  NO  lines  observed  in  the  2nd  order  of  the  grating.  Most  of  the  spectra 
were  recorded  with  a  Hamamatsu  R166  solar  blind  (Cs-Te  photocathode)  photo¬ 
multiplier  tube  cooled  to  -30°C  in  a  Products  for  Research  TE-177  thermoelectr ical ly 
cooled  housing.  Some  measurements  were  also  made  with  an  EMI  9659QB  photo¬ 
multiplier  (extended  S-20  photocathode)  cooled  to  -78°C  with  dry  ice. 
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INTENSITY/HONL  LONDON  FACTOR  (ARBITRARY  UNITS) 


NARROW-LINE  LAMP  EMISSION 


1.  Pn  (6.5)  +  Pn  (11.5) 

2.  P, ,  (7.5)  +  P1 ,  (10.5)  +  P22  (22.5)  +  Q,  2  122.5) 

3.  Pt1  (8.5)  +  Pn  (9.5) 

4.  Q22  (14.5)  +  R,2  114.5)  +  R22  0.5)  +  P,2  (31.51 

5.  P22  (21.51  +  Q12  (21.5) 
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FIG.  22 


EMISSION  FROM  FIG.  21  AFTER  ABSORPTION  BY  5  TORR  10%  NO/AR 
OVER  PATH  LENGTH  OF  18.6  cm 

1.  P,,  (6.5!  +  P,,  (11.5! 

2.  P, ,  (7.5! +P,,  (10.5)  +  P22  (22.5)  +  Q,2  (22.5! 

3.  P,,  (8.5!  +  Pn  (9.5) 

4.  Qj2  (14.51  +  R12  (14.51  +  Rjj  (9.5) 

5.  P,,  (21.5)  +  Q  (21.5) 


4 
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The  signal  from  the  photomultiplier  was  amplified  with  an  Analog  Devices 
AD310K  used  in  the  electrometer  mode  with  feedback  components  of  R,  =  100  M 

/  L 

and  Cf  =  5000pf.  The  scan  rate  was  3.95  •  10”q  nm/sec.  The  resulting  spectra 
were  recorded  with  a  Hewlett-Packard  7100B  strip  chart  recorder.  In  the  case 
of  the  continuum  lamp,  some  of  the  spectra  were  corrected  for  lamp  drift  by  a 
ratiometric  technique.  This  involved  a  reference  measurement  of  the  source  lamp 
prior  to  any  absorption,  and  was  accomplished  by  placing  a  flat  mirror  slightly 
off  the  beam  directed  through  the  calibration  apparatus  which  reflected  light 
through  a  226  nm  filter  and  onto  an  EMI  9601B  photomultiplier  tube.  The  signal 
from  the  spectrometer  was  divided  by  this  reference  signal  in  an  Ithaco  model 
3512  ratiometer  and  this  resultant  ratio  was  recorded  on  the  strip  chart  recorder. 
This  ratiometric  technique  reduced  but  did  not  eliminate  baseline  drift  in  the 
recorded  spectra  while  using  the  continuum  lamp. 

The  ultraviolet  radiation  from  the  source  lamp  was  collimated  and  directed 
through  the  12.7  mm  diameter  apertures  in  the  FGH/FFB  shroud,  across  the  FGH/FFB, 
and  then  imaged  on  the  spectrometer  slit  using  fused  silica  lenses. 

The  spectrometer  was  operated  in  high  resolution  and  low  resolution.  The 
slits  were  set  at  about  5  to  10  pm  for  the  high  resolution  work,  corresponding 
to  slit  functions  with  FWHM  of  about  .0015  nm  to  .0025  nm.  The  slits  were  set 
at  1380  pm  for  the  low  resolution  studies,  with  a  FWHM  value  of  0.146  nm. 

III. A. 3  Static  Cell 

The  static  cell  is  a  room  temperature  absorption  cell  made  of  stainless  steel 
with  valves  toward  each  end  and  UV  grade  fused  silica  windows.  The  cell  was 
leaked  checked  at  10-^  torr  (1.33-10~^  Pa.)  and  115  psig  (791  kPa,).  The  cell 
i.d.  is  22  mm  and  the  optical  path  is  18.6  cm. 
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B.  Theoretical  Development  of  Ultraviolet  Absorption 
III.B.l.  Necessity  for  Theoretical  Model 

A  detailed  theoretical  model  has  been  developed  to  describe  the  absorption 
of  ultraviolet  radiation  by  nitric  oxide  in  the  (0,0)  band.  It  would  not  be 
unreasonable  to  question  the  necessity  for  such  a  first  principles  model,  when 
the  possibility  exists  for  an  empirical  calibration  technique.  But  how  well  can 
known  concentrations  of  NO  be  generated  over  a  temperature  range  of  300  K  to 
2000  K  and  a  pressure  range  of  0.5  to  1.5  atm?  The  answer  is  that  it  is  diffi¬ 
cult  and  expensive  to  generate  a  complete  range  of  conditions  necessary  for 
repetitive  empirical  calibration.  The  approach  adopted  in  this  work  was  to  start 
with  a  theoretical  model  and  then  verify  that  model  by  testing  it  over  a  range  of 
conditions,  similar  to  the  approach  used  by  Davis  et  al.  (1976a)  but  with  addi¬ 
tional  information  from  high  resolution  continuum  absorption  data. 

An  additional  problem  with  empirical  models  is  that  they  are  often  instru¬ 
ment  dependent  and  quite  difficult  to  transfer  from  a  given  instrument  at  one 
laboratory  to  instruments  in  general.  An  example  for  this  particular  experiment 
would  be  that  a  strictly  empirical  calibration  would  depend  on  the  characteristics 
of  the  narrow-line  source  (rotational  distribution)  and  the  spectrometer  (PM  tube 
response  and  slit  function)  which  would  make  duplication  of  results  at  various 
laboratories  difficult.  The  theoretical  model  does  include  some  empirical 
"adjustments,"  as  will  be  discussed,  and  several  experimentally  determined  input 
parameters,  but  the  absorption  process  is  modeled  from  basic  principles. 

III.B.2.  Transmission  Formulas,  Narrow-Line  Source 

The  general  problem  of  relating  measured  absorptions  in  the  y(0,0)  band  to 
molecular  densities  of  NO  may  be  summarized  as  follows.  In  order  to  predict  low 
resolution  absorption  profiles  in  the  y(0,0)  band,  it  is  necessary  to  model  on 
the  order  of  500  lines  whose  relative  intensities  change  with  temperature,  and 
whose  shape  changes  with  temperature,  pressure,  and  gas  composition.  These 
lines  are  so  closely  spaced  that  even  when  recorded  with  a  high  resolution 
instrument,  many  of  the  lines  overlap. 

The  theory  for  the  transmission  of  a  group  of  Doppler  broadened  NO  source 
lines  through  an  absorbing  gas  with  both  Doppler  and  collision  broadened  lines 
has  been  developed  by  Davis  et  al.  (1976).  We  will  use  the  results  of  Eqs.  (1) 
through  (14)  of  Davis  et  al.  (1976),  which  can  be  summarized  as  follows.  The 
total  transmitted  intensity  for  the  jth  emission  line  (T^)  due  to  absorption  by 
many  lines  is  related  to  the  line  center  source  intensity  of  the  jth  line  (I°v?), 
the  Doppler 1  width  of  the  source  (AgV^)^,  the  line  center  Doppler  absorption 
coefficient  of  the  gas  (kv  ) ,  the  broadening  parameter  (a'),  and  the  path  length 
( £ )  by, 


III-10 


oo 

' 

Hi 

l— u 

11 

I°o 

Vj 

/  exp 

o 

— 

^-V°j-lv/lV2 

(As  Oi 


'•  exp 


po  2 

-  —  T'  k  °  /  a1  e~y  dv 

n  ~  v i  /  a '2  +  'ui.  _ 

i  1  -'-oo  v  i 


(20) 


y)‘ 


d  v 


where  a'  is  related  to  the  Lorentz  full  width  at  half  maximum  (FWHM)  in  the 
absorbing  gas  (Aav^)^  and  the  corresponding  FWHM  of  the  Doppler  component 

(Aavi}D 


<AaVL 

<AaVD 


yin  2 


(21) 


and  where 


2  (vj  ~  vj) 

(AaVi>D 


yin  2 


(22) 


The  Dopnler  widths  (FWHM)  of  the  source  and  absorbing  line  are  given  by 


(A  Vj)D 


°  2  /  2  (In  2)  k  T  _  °  _L  /2(ln  i)-Rl 

Vj  c\j  m  Vi  C  V  M 


(23) 


O 

where  is  the  wavenumber  of  line  center,  c  is  the  speed  of  light,  k  is 
Boltzmann's  constant,  T  is  the  appropriate  temperature  (either  source  or 
absorber),  m  is  the  molecular  mass  in  grams  of  the  emitting  cr  absorbing 
molecule,  R  is  the  universal  gas  constant,  and  M  is  the  molecular  weight  in 
grams  per  mole.  The  Lorentz  FWHM  due  to  collisions  is  given  by. 
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where  is  the  frequency  of  collisions  between  an  absorbing  molecule  and  the 
perturbing  molecules,  and  t  is  1/Z^  or  tbe  mean  tl®6  between  collisions. 

It  should  be  noted  that  this  development  assumes  that  the  collision  process 
and  line  broadening  can  be  accurately  explained  by  the  Lorentz  collision  theory, 
an  assumption  which  is  discussed  in  detail  later  in  this  report.  The  use  of  the 
Lorentz  theory  to  examine  experimental  results  is  certainly  the  usual  assumption, 
but  deviations  from  Lorentz-theory  behavior  have  been  observed  in  numerous 
instances  (see  Section  III.B. 5  on  broadening  theory). 
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It  is  worthwhile  considering  the  factors  in  Eq.  (20)  which  must  be  known 

or  determined  in  order  to  predict  the  transmission.  The  relative  intensities 

of  the  source  lines  at  line  center  must  be  measured  for  the  narrow-line  lamp, 

and  intensities  of  lines  that  cannot  be  measured  directly  must  be  interpolated 

from  a  Boltzmann  plot  such  as  shown  in  Fig.  20.  The  Doppler  FWHM  ((AVj)^)  may 

easily  be  calculated  from  Eq.  (23).  The  path  length  (?)  is  easily  determined. 

The  broadening  parameter  (a1)  must  be  determined  experimentally  because 

theoretical  computation  of  the  Lorentz  FWHM  (  ^avi^L^'  which  is  necessary  to 

determine  a'  through  Eq.  (21)  is  not  feasible.  The  experimental  technique 

used  to  determine  a'  is  described  in  a  later  section.  The  final  parameter 

which  must  be  determined  to  solve  Eq.  (20)  for  the  transmission  is  the  absorption 

•  0 

coefficient  at  line  center  for  a  Doppler  broadened  line  (hv.).  The  relation 
between  kv.j  and  the  density  of  molecules  in  the  initial  state  and  the  oscillator 
strength  is  described  below. 


II. B. 3  Computation  of  Doppler  Absorption  Coefficient 


The  line  center  absorption  coefficient  for  a  Doppler  broadened  line 
(kvp  is  given  by  Mitchell  and  Zemansky  (197i) , 


where  e  is  the  electron  charge,  mp  is  the  electron  mass,  c  is  the  speed  of 
light,  is  t*ie  oscillator  strength  of  the  line,  and  N(n"v"Z"J"p") 

is  the  population  of  the  initial  energy  state.  The  single  and  double  primes 
refer  to  upper  and  lower  electronic  state  parameters,  respectively. 


The  Doppler  FWHM  ((AavpD)  can  be  calculated  from  Eq .  (23).  The 
oscillator  strength  for  a  line  (fv, v"j ' K»p»)  is  re^ate(^  to  the  oscillator 
strength  for  a  band  (fv'v»)  by  (Thorne  (1974)), 
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v"J'J"K'K"p" 


Sj' J"K'K"p" 


2J"  +  1 


f  yj  *  y*1 


(26) 


where  Sji j"k'K"p"  t*le  Honl-London  factor  and  J"  is  the  rotational  quantum  number 
in  the  lower  electronic  state.  The  band  oscillator  strength  (f v » v”)  for  the 
NO  y(0,0)  band  has  been  measured  previously  and  another  determination  is 
reported  in  this  work  (see  experimental  results  section).  The  Honl-London 
factors  may  be  computed  from  equations  given  by  Earls  (1935),  but  with  an 
adjustment  factor  for  proper  normalization  (multiply  Earls  values  by  4).  It 
is  important  to  note  that  the  summation  convention  for  the  Honl-London  factors 
suggested  by  Tatum  (1967)  has  been  adopted  here,  that  is: 
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(27) 


£ 

sub-levels 


(2-50>a„)  (2S  +  1)  (2J"  +  1) 


where  the  summation  over  sub-levels  includes  both  spin  splitting  levels 
(the  2S  +  1  factor)  and  A  doubling  levels  (the  2 -SQ  A  factor) .  The  value 
50,^is  unity  for  A  =  0,  but  zero  otherwise.  For  the  NO  y(0,0)  band,  the 
summation  convention  for  the  Honl-London  factor  becomes 

£  £  SJfJ)I  =  A  ( 2J"  +  1)  (28) 

sub-levels  J' 


from  the  fact  that  (2S  +  1)  =  2,  and  (2-6o,A„)  =  2  (see  later  section  on  the 
spectroscopy  of  NO  for  an  explanation  ot  quantum  numbers). 


The  remaining  factor  in  Eq.  (25)  which  must  be  determined  in  order  to 
solve  for  is  represented  by  the  somewhat  awkward  but  definitive  notation 

N(n"v"  E"J"p")  ,  which  is  the  population  of  the  energy  level  of  interest  for  a 
molecule  obeying  Hund's  case  (a).  The  more  often  used  symbol  Nj,,  is  sometimes 
ambiguous  in  that  it  may  or  may  not  include  the  factor  of  1/2  for  the  A  doubled 
levels.  Thus  N(n”v”E"j"p")  represents  the  number  of  molecules  in  the  electronic 
state  n" ,  vibrational  state  v",  spin  state  I",  rotational  level  J",  and  with 
specified  parity  p"  (+  or  -) .  The  fraction  of  the  total  number  of  molecules 
which  are  in  a  given  level  is  represented  by  N(n"v"E"J"p") /N,  and  can  be 
determined  for  the  NO  '2-^-yj'i  and  2  sublevels  of  the  ground  state  as  follows. 


The  ground  state  of  the  NO  molecule  is  intermediate  between  Hund's  cases 
(a)  and  (b) ,  but  to  a  close  approximation  follows  Hund's  case  (a).  Following 
Tatum's  (1967)  convention  and  notation  for  Hund's  case  (a)  for  the  state 

(all  terms  are  for  the  lower  electronic  state  and,  hence,  should  have  double 
primes) , 


N(nvEJp)  _ 
N 
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2(2S  +  1)  exp  ((-hc/kT)  T0)  J  j'exp  ((-hc/kT)  G(v)) 

E  2(2S+1)  exp (( -hc/kT) T0)|  [Ey  exp  ((-hc/kT)  G(v)) 

all  states 


exp  ((-hc/kT)  (-A/2) ) _ 

exp  ((-hc/kT)  (-A/2))  +  exp  ((-hc/kT)  (A/2)) 


(29) 


(2J  +  1)  exp  ((-hc/kT)  F(J)) 

Ej  (2j  +  1)  exp  ((-hc/kT)  F(J)) 
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where  S  is  the  electron  spin  quantum  number,  h  is  Planck's  constant,  Tq  is  the 
electronic  term  value  when  spin  is  neglected,  G(v)  is  the  vibrational  term 
value,  A  is  the  spin  splitting  constant,  and  F(J)  is  the  rotational  term  value. 
The  equivalent  expression  for  the  2tt^^ ^  state  is  identical  except  that  the 
numerator  in  the  third  term  in  brackets  becomes  exp  ( (-hc/kT) (A/2)) . 

The  first  term  in  brackets  in  Eq.  (29)  is  the  fractional  population  of  a 
given  electronic  level,  N(n)/N,  and  differs  from  that  given  by  Tatum  in  that 
we  have  used  the  "average"  electronic  term  value  Tq  (or  electronic  term  value 
when  spin  is  ignored)  rather  than  Tatum's  Te  where  Tg  =  Tn  +  AAI  (for  2tt1  /  2  ’ 

Te  =  Tq- (A/2) ;  for  2^3/2’  Te  =  To  +  (A/2)).  The  use  of  Tg  in  place  of  Tc 
would  produce  (for  the  N0y(0,0)  band),  a  term  essentially*  identical  to  the 
third  term  in  brackets  in  Eq.  (29),  but  following  Tatum's  guidelines,  the 
fractional  populations  of  the  spin-split  levels  are  to  be  calculated  separately 
and  explicitly  as  N(nvZ) /N (nv) ,  the  third  term  in  brackets.  Thus,  it  would 
appear  that  Tatum  has  accounted  for  the  fractional  populations  of  the  spin- 
split  levels  twice,  which  we  have  avoided  by  using  TQ  in  Eq.  (29).  For  the 
electronic  term  value,  as  expressed  in  Eq.  (29),  the  first  term  in  brackets 
is  unity. 

The  second  term  in  brackets  in  Eq.  (29)  is  the  fractional  population  of  a 
given  vibrational  level,  N(nv)/N(n),  which  for  the  v"  =  0  level  is  close  to 
unity  at  300  K  but  decreases  to  about  0.74  at  2000  K.  The  third  term  in 
brackets  in  Eq.  (29)  determines  the  fractional  population  N(nv£)/N(nv)  of  the 
spin-split  2tti/2  anc*  2lT3 / 2  2evels<  The  fourth  term  in  brackets  represents  the 
fractional  population  of  a  given  rotational  level  N(nv£J)/N(nv£) .  By 
employing  the  usual  approximations  (Herzberg  (195o))»  the  sum  in  the  denom¬ 
inator  of  the  fourth  term  in  brackets  may  be  written  as  kT/hcB.  •  The  fifth  term, 
symbolically  written  N(nvIJp)/N(nvEJ) ,  accounts  for  the  A-type  doubling  of 
the  rotational  levels  of  the  state  into  two  sublevels  of  parity  +  and 
Because  the  upper  electronic  state  being  considered  is  a  state,  only  one  of 
the  two  sublevels  is  active  for  a  transition  to  a  specified  upper  state 
rotational  level  (Herzberg  (1950)). 


*  This  assumes  that  exp  ((-hc/kT)Tg)  <<  1  for  all  electronic  states  higher 
than  the  ground  state  (2it),  which  is  true  for  the  NO  molecule  even  at  high 
temperature. 
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Thus,  Eq. 


(29)  for  the 


state  may  be  rewritten: 
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N 


exp( (~hc/kT)G(v) ) 
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(30) 


o 

The  equivalent  expression  for  the  71 3/ 2  state  is: 


N(nvLJp)  ^ 

exp ( (-hc/kT) G(v) ) 
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exp ((-hc/kT)  (A/ 2)) 
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Iv  exp( (-hc/kT)G(v) ) 

exp((-hc/kT)  (-A/2))  +  exp ( (-hc/kT) (A/2) ) 

hcB  (2J  +  1)  exp  ((-hc/kT)  F  (J))  (31) 

2kT 


These  expressions  may  be  modified  intermediate  to  Hund's  case  (a)  and  (b)  by 
absorbing  the  spin-splitting  term  into  the  expression  for  F(J)  (Hill  and  Van 
Vleck  (1928))  and  dropping  the  third  term  in  brackets  on  the  right-hand  side  of 
Eq.  (29).  When  this  is  done,  however,  the  rotational  partition  function  may  no 
longer  (in  general)  be  approximated  as  kT/hcB.  Assuming  that  the  spin-splitting 
energy  is  absorbed  into  F(J)  and  the  resultant  is  labeled  Fi(J),  with  i  =  1  for 

the  ,  state  and  i  =  2  for  the  2  state,  then  instead  of  Eqs.  (30)  and  (31), 
1/2  77  3/2 

one  gets, 


N(nv(£J)p) 
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exp  ((-hc/kT)  G(v)) 

Sv  exp( (-hc/kT)  G(v)) 


x 


(2J  +  1)  exp  ((-hc/kT)  F1  (J) _ 

Si  =  1,  2  Ej  (2J+1>  ((-hc/kT)  F1  (J))  J 


(32) 


For  the  NO  molecule,  we  have  determined  that  the  rotational  partition  function 
can  be  approximated  with  high  accuracy  (*  0.3%  for  300  K  <  T  <  2000  K)  as 

almost  Hund's  case  (a), 
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(33) 


£i  -  1  2  ST  (2J+i)«*p((-hc/kT)Fi  (J))  2(kT/hcB)  [ exp( (-hc/kT) (-A+2B) /2) 

1  I 

+  exp(-hc/kT) (A-2B)/2) I 

for  the  energy  equations  and  constants  as  given  by  Engleman,  et  al  (1970).  This 
approximation  was  subsequently  verified  by  D.  Keefer.  Thus,  Eq.  (32)  for  the 
intermediate  Hund's  case  between  (a)  and  (b)  becomes, 

N(nv(EJ)p)  ^ 

N 

x 


X 


If  Eqs.  23,  26,  and  34  are  substituted  into  Eq.  (25),  the  result  is 
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exp((-hc/kTa)(-A+2B)/2)+exp((-hc/kTa)(A-2B)/2) 


■) 


or,  substituting  for  the  constants, 
kvi  = 


=  1.755-10"10  f  ,  „  s 

v'v"  bj.j 


r3/2 
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exp(86.18/Ta)+exp(-86.18/Ta 
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(36) 


With  Eq.  (36)  for  kv®,  all  terms  necessary  to  calculate  the  transmission 
(Tj)  using  Eq.  (20)  have  been  determined  except  for  the  broadening  parameter  a' 
and  a  confirmation  of  the  band  oscillator  strength  fy' 

I1I.B.4,  Transmission  Formulas,  Continuum  Source 


The  development  of  the  transmission  formulas  for  continuum  radiation  is  very 
similar  to  that  used  for  the  narrow-line  case  given  in  Eq.  (20).  The  source  func¬ 
tion  is  much  simplified  to  just  I(v),  which  Is  a  very  slowly  varying  and  easily 
measured  function.  The  transmitted  intensities  must  be  recorded  continuously  as 


a  function  of  wavenumber  rather  than  as  discrete  intensities  at  the  line  center 
of  emission  lines,  as  is  done  for  the  narrow-line  source.  The  transmission  at  any 
wavenumber  v  is. 


T(v) 


I(v)exp 


I  k  °  "a’e-y2dy _ \ 

i  j  a'2+(ui-y)2/ 
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where  the  symbols  are  the  same  as  explained  previously  for  Eq.  (20). 


III.B.5  Broadening  Theory 

There  is  no  lack  of  theoretical  models  available  to  treat  pressure  broaden¬ 
ing.  The  theories  of  pressure  broadening  are  discussed  in  some  detail  by 
Breene  (1957),  Hindmarsh  and  Farr  (1972),  and  Mitchell  and  Zemansky  (1971). 
Thorne  (1974)  discusses  the  subject  with  less  depth  but  good  clarity.  An 
often  quoted  work  for  microwave  and  infra-red  regions  is  the  paper  by 
Anderson  (1949). 


In  spite  of  the  number  of  theoretical  models  available,  experimentalists 
tend  to  use  the  standard  Lorentz  theory  and  label  discrepancies  as  deviations 
from  the  Lorentz  theory,  rather  than  attempting  to  use  other  models.  In  this 
report,  results  from  the  Lorentz  theory  and  Weisskopf  theory  are  used,  but 
the  reader  is  referred  to  the  previous  references  for  the  details  of  these 
theories. 


The  Lorentz  model  assumes  that  the  optical  collision  diameter  is  indepen¬ 
dent  of  the  relative  velocities  of  the  two  colliding  molecules  (or  atoms). 

From  the  Lorentz  theory  the  FWHM  due  to  collisions,  Av^,  is  given  by  Mitchell 
and  Zemansky  (1971)  , 


Avl 


(38) 


where  Z^  is  the  collision  frequency  of  a  single  molecule,  c  is  the  velocity 
of  light,  and  Av^  is  in  wavenumbers.  It  should  be  noted  that  this  is  true 
only  for  collisions  between  molecules  of  different  types.  If  the  molecules 
are  all  of  the  same  type  then  Fq.  (38)  is  multiplied  by  a  factor  of  2.  This 
is  due  to  the  fact  that  each  collision  terminates  one  mean  free  path  if  the 
collision  is  between  the  gas  being  studied  and  a  foreign  gas,  but  it  termi¬ 
nates  two  mean  free  paths  if  the  collision  is  between  two  molecules  of  the 
same  type.  The  interest  here  is  in  collisions  between  a  minor  species  and 
foreign  gases,  so  Eq.  (38)  applies. 

The  total  number  of  collisions  per  second  between  molecules  of  type  1 
and  type  2  with  number  densities  n1  and  n2  and  molecular  weights  (molar)  Mj 
and  M2  is  given  by  Mitchell  and  Zemansky  (1971)  as, 
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where  d  is  the  hard  sphere  collision  diameter  (d  *  ri+r2»  where  r  is  the  radius),  R 
is  the  universal  gas  constant,  and  T  the  absolute  temperature.  The  collision  fre¬ 
quency  per  molecule  of  a  given  type  is. 


ZL  ‘ 


^  '  2"2dl-2  V  ^  + 


Combining  Eqs.  (38),  (40),  (21) 


2  r~ 

V>  7T 


2  r  mi 

n2dl-2  *  1  + 


(40) 


(41) 


where  v  is  the  wavenumber.  Expressing  n2  in  terms  of  the  partial  pressure  (atm) 
of  the  foreign  gas  and  temperature  (K) ,  d  in  cm,  and  v  in  cit~1 
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2  /  M, 

Jl-2  V  1  -  ^ 


h 


(42) 


or 


a  1  =  C0  12 
2  2  f 


(43) 


If  more  than  one  foreign  gas  is  active  in  broadening,  Z L  is  summed  over  all 
the  species  and  a!  for  each  species  is. 


c.  Li 

1  t 


(45) 


where  is  the  partial  pressure  of  the  ith  species.  The  total  a'  is  the 

sum  of  the  a's. 

1 


An  important  point  must  be  made  regarding  the  collision  diameter  and  the 
effective  collision  cross-section.  The  effective  optical  collision  diameter 
is  given  by  d  in  the  equations  above.  References  in  the  field  of  kinetic 
theory,  such  as  Reif  (1965),  and  Hirschfelder  et  al  (1954),  call  ird  the 


effective  cross-section,  while  Mitchell  and  Zemansky,  Davis  et  al  (1976), 

O 

others  drop  the  it  and  refer  to  d  as  the  effective  cross-section.  Thus, 
caution  must  be  exercised  in  comparing  effective  cross-sections  given  by 
different  authors.  We  define  cross-section  as  ird^  in  this  report. 


and 


According  to  the  Lorentz  theory,  at  constant  pressure,  Av^  is  proportional 
to  T~®-5  and  a'  is  proportional  to  T.-*’^.  There  is  a  significant  amount  of 
evidence  to  indicate  that  the  temperature  dependence  is  stronger  than  this. 

Hence,  the  temperature  dependence  of  Av^  should  be  between  T”®’^  and  T-^*®  and 
that  of  a'  between  T-1*®  and  T-^’^.  Among  the  evidence  for  a  stronger  temp¬ 
erature  dependence  for  neutral  atoms  or  molecules  is  work  by  Hansen  (1978), 

Townes  and  Schawlow  (1955),  Engleman  (1969),  Cann  et  al  (1979),  Planet  et  al  (1978), 
and  Breene  (1967) . 


The  Weisskopf  theory  gives  a  stronger  temperature  dependence  in  agree¬ 
ment  with  a  significant  amount  of  experimental  work  (Thorne  (1974)). 

Weisskopf  developed  a  collision  theory  which  includes  the  effect  of  collision 
time  in  computation  of  the  cross-section.  For  collisions  between  neutral 
atoms  or  molecules,  the  Weisskopf  theory,  for  a  constant  pressure,  gives  the 
collision  width  as  Avw  proportional  to  T-<^' 2  and  a'  proportional  to  T.“^‘2 
This  is  in  almost  exact  agreement  with  the  result  reported  by  Hansen  (1978) 
of  Av  proportional  to  T-®*2^  for  CO  broadening  over  a  wide  range  of  temperatures. 
The  evidence  from  the  flame  measurements  reported  here  indicates  that  the 
temperature  dependence  of  the  Weisskopf  theory  seems  to  fit  the  experimental 
data  for  NO  broadening  better  than  the  Lorentz  theory. 

III.B.6.  Spectroscopic  Details  of  the  Nitric  Oxide  Ultraviolet  Y  (0,0)  Band 
1 1 1 .  B .  6 .  a  .  _Band_S^stems_of^  N0_ 

The  electronic  transitions  of  the  NO  molecule  are  observed  in  the  ultraviolet 
and  vacuum  ultraviolet  spectral  regions.  The  lower  energy  transitions  involving 
the  ground  (X2tt)  state  are  labeled  as  follows  (Pearse  and  Gaydon  (1965)). 


Transition  Type 

Label 

A2I+  -  X2ir 

Y 

system 

B2tt  -  X2t r 

B 

system 

C2ir  -  X2tt 

6 

system 

D2E+  -  X2u 

e 

bands 

The  excited  electronic  states  are  labeled  in  order  of  increasing  energy  by  A,  B, 

C,  etc.,  and  therefore  the  lowest  energy  (longest  wavelength)  transitions  are  in 
the  y  system,  which  makes  that  system  the  most  easily  accessible  for  absorption 
measurement  of  NO. 

The  bandheads  for  some  of  the  lower  vibrational  levels  are  shown  in  Table 
III-B.  For  determining  NO  concentrations  from  spectroscopic  absorption  measure¬ 
ments  at  temperatures  below  1000  K,  it  is  only  necessary  to  consider  those  states 
originating  from  a  ground  vibrational  level  (v"  -  0).  The  longest  wavelength 
transition  of  this  type  is  the  Y  (v'»0,  v"=0)  band  at  -226  nm.  A  number  of  exper¬ 
imental  problems  become  more  severe  at  wavelengths  shorter  than  this:  (1)  the 
transmission  of  quartz  decreases,  (2)  hot  C02  and  O2  bands  absorb  more  strongly 
at  flame  temperatures  (3)  the  chance  of  interfering  absorptions  by  other  species 
increases,  (4)  intense  continuum  sources  such  as  high  pressure  Xe  and  Hg-Xe 
lamps  decrease  rapidly  in  Intensity,  and  (5)  scattering  by  soot  particles  gen¬ 
erally  increases.  Thus,  the  Y  (0,0)  band  is  the  most  reasonable  band  selection 
for  the  measurement  of  NO  over  the  temperature  range  from  300  K  to  2000  K. 
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TABLE  III-B 


Bandheads  in  the  Y  -System  of  Nitric  Oxide 
(Pearse  and  Gaydon  (1965)) 

Double  bandheads  are  given  for  each  band 


(run) 

247.87 

247.11 

244.70 

244.00 

237.02 

236.33 

231.63 

230.95 

228.98 

228.41 

226.94 
226.28 

224.54 

223.94 

222.24 

221.63 

219.96 
219.40 

215.49 

214.91 


v '  ,v" 
0,2 

1.3 

0,1 

2.3 

3.4 

0,0 

1,1 

2,2 

3,3 

1,0 
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The  Y  (0,0)  band  represents  transitions  between  the  zeroth  vibrational  levels 
of  the  excited  state  A^I+  and  the  ground  state  X^ti  .  The  ground  electronic  energy 
level  of  the  NO  molecule  is  split  into  two  sub-levels,  the  lower  energy  state 

and  the  higher  energy  ^3/2  state,  which  are  widely  split  in  energy  (-120  cm~^). 
Because  of  the  considerable  difference  in  energy  of  the  two  sublevels,  the 
state  of  NO  is  usually  treated  as  following  Hund's  case  (a),  but  is  actually 
intermediate  to  cases  (a)  and  (b) ,  approximating  case  (a)  for  small  rotational 
velocities,  but  undergoing  some  spin  uncoupling  at  higher  rotational  rates.  The 
upper  A-^+  electronic  state  belongs  strictly  to  case  (b)  coupling,  as  do  all  I 
states . 


1 1 1 . B .  6 . h .  Ener^y^Levels 

The  equations  and  constants  used  here  to  determine  the  energy  for  various 
rotational  levels  in  both  upper  and  lower  electronic  states  are  almost  exclusively 
from  Engleman,  et  al.  (1970).  Energy  levels  of  a  molecule  (exclusive  of  trans¬ 
lation)  may  be  approximated  as  the  sum  of  the  rotational  (Er) ,  vibrational  (Fv), 
and  electronic  energies  (Ee) , 

E  =  Ee  +  Ev  +  Er  (4b) 

It  is  common  practice  to  express  the  molecular  energies  as  term  values,  which 
are  energy  values  divided  by  he,  and  have  units  of  cm“^.  Using  the  term  value 
notation,  Eq.  (46)  becomes. 


T  =  Te  +  G(v)  +  F(J)  (47) 

In  an  electronic  transition  the  wave  number  (cm”l)  of  a  spectral  line  is  given  by 
the  difference  of  the  term  value  in  the  upper  electronic  state  (T')  and  the  lower 
electronic  state  (T")  , 

u  =  T’  -  T"  -  (V  -  Te")  +  (G'(v')  -  G"  (v">)  +  (Fv •  (J ' ) -Fv» (J") )  (48) 

1 1 1 .  B .  6 .  c  Jiner j[y_Levie  1_E^U£  t  ioni  s_f£r_t  jie_x£ir_S  ta  te_o f_  N0_ 

For  the  state  of  NO,  the  electronic  energy  and  rotational  energy  are 
slightly  coupled  so  that  Te  and  F(J)  are  not  independent.  The  term  value  energy 
equation  is  modified  so  that  F^(J)  contains  not  only  the  rotational  energy,  but 
also  the  energy  associated  with  the  spin  splitting  in  the  2-n  state,  and  Tq  con¬ 
tains  the  electronic  energy  exclusive  of  spin-splitting, 

T"  =•  T  "  +  G"(v")  +  F,"(J")  (49) 

o  1 

where  the  vibrational  term  value  G(v)  is. 
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(50) 


G(v)  =  (v+0.5)  -  weXe  (v+0.5)2 

and  the  rotational  term  F^(J)  is, 

Fi(J)  =  Bv<(J+0.5)2-1)-Dv((J+0.5)4-(J+0.5)2+1)+Bv  Ja  (51) 

a  =  Ofv-2)2+  ( (J+0.5)2-l)  [l+2u(2(J+0.5)2-X)+p2  ((2(J+0.5)2-l)2-ljj  (52) 

Yv  =  [a  +  C(J-0.5)2j  /Bv  (53) 

V  =  Dv/Bv  (54) 

2 

The  minus  (-)  sign  is  for  the  state  in  which  J  *  K  +  0.5  and  i  *  1,  while  the 
plus  (+)  sign  is  for  the  ^3/2  state  in  which  J  =  K  -  0.5  and  i  *  2.  It  should  be 
noted  that  our  equation  for  Fi(J)  includes  a  term  M-Bv"  not  found  in  Engleman, 
et  al.  (1970),  but  which  should  be  included  following  the  derivations  from  the 
original  papers  cited  by  Engleman.  We  have  adjusted  the  upper  state  electronic 
energy  given  by  Engleman,  et  al.  accordingly.  We  have  likewise  added  to  the 
equation  for  a  a  term  "u2,”  but  this  is  numerically  inconsequential.  There  also 
appears  to  be  a  typographical  error  in  the  sign  convention  given  by  Engleman, 
et  al. 


In  addition  to  the  spin-orbit  splitting  of  the  2ir  state,  each  of  the  levels 
is  subject  to  a  slight  perturbation  by  A-type  doubling  due  to  the  interaction 
between  the  rotation  of  the  nuclei  and  the  electronic  orbital  angular  momentum 
L  (Herzberg  (1950)).  The  magnitude  of  this  shift  is  given  by  <(>,  where. 


<t>  =  0.5  (J+0.5) 


(irlx  +  i)(-+  Q  )+  Q_((J+0-5)2-i) 

2/a  ^  t'ct 


(55) 


and  the  sign  convention  is  the  same  as  above.  Because  of  A-type  doubling,  each 
rotational  level  is  split  into  two  components,  but  for  a  given  transition  only 
one  of  the  components  is  active.  The  reason  for  this  is  that  the  two  components 
are  of  opposite  parity  (+  and  -) .  Because  of  the  selection  rule  +  —  -,  only  one 
of  the  components  may  interact  with  an  upper  state  with  specified  parity  (the 
upper  state  is  not  A-type  doubled) .  The  A-type  doubling  has  an  influence  on  the 
spectral  line  locations  as  follows. 


ic 


id 


Fi  ' 
Fi  + 


(56) 


where  Fic  is  the  energy  of 
transitions  P^,  P22’  ^12’ 
transitions  P12,  P21,  ^n* 


2 

the  v  state  (exclusive  of  vibrational  energy)  for  the 
Q21 ,  R12,  and  R22»  and  Fid  is  the  2ir  state  energy  for 
Q22»  ^^2*  and  ^21* 
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The  zero-point  energy  of  a  molecule  is  arbitrary  in  that  only  energy  differ- 
ances  have  significance.  For  the  constants  given  by  Engleman  et  al.,  if  the 
vibrational  and  rotational  (but  not  spin-splitting)  energies  are  set  to  zero,  the 
level  of  zero  electronic  energy  is  midway  between  the  energies  of  the  two  spin 
split  states  (  and  ^3/2^  '  addition,  the  spin  splitting  term  A  approached 

zero,  the  two  spin  split  states  would  collapse  into  one  energy  level  defined  as 
the  zero  electronic  energy  for  this  system.  However,  we  have  shifted  this  zero 
point  energy  by  A/2  to  correspond  approximately  to  the  zero  level  energy  of  the 
lower  2r^  state.  Thus  Tq  =  A/2,  and  Eq .  (49)  becomes 


T"  =  —  +  G" (v")  +  1C  (57) 

2  I  F"  (J") 

id 

and  the  upper  electronic  state  energy  given  by  Engleman  et  al.  is  increased  by 
A/2. 

The  constants  which  describe  the  lower  electronic  state  are  given  in  Table 
III-C.  The  A  doubling  constant  P  has  a  different  sign  in  Engleman  et  al . ’ s 
Tables  VII  and  XI.  We  have  chosen  the  negative  sign  because  it  gives  better 
agreement  between  the  experimental  and  theoretical  line  locations. 

2  + 

III . B. 6. d  ^ner^y_Level_E£U£ti.ons_f£r_the_A_I_  States  of_NO 

2  *4- 

The  energy  levels  of  the  upper  electronic  state  (A  I  )  are  described  by 
the  general  Eq .  (47) 

T’  =  T;  +  G(v’)  +  F’(J’)  (58) 


where  Tg  is  a  constant,  G(v)  is  given  by  Eq.  (50),  and  F(J)  is  given  by  Engleman 
et  al.  as, 

FX(J)  -  Bv(J+0.5)(J-0.5)-Dv(J+0.5)2(J-0.5)2+0.5y(J-0.5)  (59a) 

F2(J)  =  Bv(J+0. 5) (J+l . 5)-Dv(J+0. 5) 2(J+1 . 5) 2-0. 5y( J+l . 5)  (59b) 

where  FX(J)  is  for  states  with  J  =  K  +  1/2  and  F2(J)  is  for  states  with  J  =  K  -  1/2. 
Of  course,  all  values  of  J,  K,  and  v  are  for  the  upper  electronic  state  and  would 
have  single  primes  if  any  were  indicated.  The  y  is  a  small  spin  splitting  constant 
and  the  value  used  was  determined  by  Bergeman  and  Zare  (1972)  from  an  rf  resonance 
study.  The  constants  used  to  describe  the  upper  state  are  given  in  Table  III-D. 
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TABLE  III-C 


X  n  State  Constants  for  NO 
(from  Engleman,  et  al.  (1970)) 


T  =  61.595  cm-1* (=A/2) 

o 

-1 


123.19  cm 
-5.8.10-4  cm-1 


( i).  =  1904.405  cm“l 

e 


cjexe=  14.1870 

cm-^ 

B  =  1.69568 

V 

cm~^ 

Dv  =  4.5-10-6 

cm~^ 

P  =  -1.17-10' 

-2 

cm 

Q  =  7.8-10-5 

-1 

cm 

>  text 
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TABLE  III-D 


? 


i 


A^E  State  Constants  for  NO 
(from  Engleman,  et  al.  (1970)  except  as  noted) 


Tg  =  43966 . 2643  cm~1+ 
we  =  2374.307  cm-1 

aigX  e  =  16.106  cm“l 

Bv  =  1.98576  cm-1 

Dy  =  4.6 -10-6  cm-1 

V  =  .00276  cm“l* 


+adjusted  for  different  zero  references 
from  Engleman,  et  al. 

*from  Bergeman  and  Zare  (1972). 
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Ill . B . 6 . e.  MO  y  (0,0)  Band  Structure 


For  cases  where  the  2ir^  and  ^3/2  ener@y  levels  are  noticeably  separated 
(Hund's  case  (a)),  the  band  structure  of  a  2E  -  2tt  transition  is  usually  considered 
as  two  sub-bands,  2Z  -  2tti^  and  2Z  -  ^ tt 3 / 2 «  which  are  separated  from  each  other  by 
the  amount  of  doublet  splitting  of  the  2tt  state.  There  are  six  branches  possible 
for  each  of  two  sub-bands,  making  a  total  of  twelve  branches.  Because  of  the 
small  spin  splitting  in  the  2Z  state  of  NO,  only  eight  distinct  branches  are  re¬ 
solved,  but  in  this  study  the  spin  splitting  in  the  2Z  state  has  been  included 
based  on  the  value  given  by  Bergeman  and  Zare  (1972). 


The  selection  rules  which  determine  which  upper  and  lower  electronic  states 
are  involved  in  transitions  are  explained  in  detail  by  Herzberg  (1950),  and  may 
be  summarized  as  follows.  For  changes  in  vibrational  energy,  Av  =  v'  -v  =  0,  -  1, 

±  2,  -  3...,  with  the  relative  strengths  determined  by  the  overlap  of  the  potential 
energy  curves.  These  relative  strengths  are  given  by  the  Franck-Condon  factors. 

Concerning  rotational  levels,  only  states  of  opposite  parity  may  interact,  +~  - 
+  ^.+  This  has  an  effect  on  which  of  the  A-type  doubled  levels  in  the  tt 

state  may  interact  with  given  upper  state  levels.  The  total  angular  momentum  J 
is  restricted  such  chat  AJ  =  J’  -  J"  =  0,  ±  1.  Since  the  2*  state  is  intermediate 
between  Hund's  cases  (a)  and  (b) ,  the  selection  rule  AN  =  N'-N”  -  0,  i  1  (ol<^ 
notation  AK  =  K'-K")  applies,  but  AN  =  ±  2  can  appear  subject  to  AJ  =  ±1,  although 
the  intensity  is  much  reduced.  Branches  for  which  AJ  *  AN  will  genera  y  e  wea  er. 


These  selection  rules  imply  that  there  are  twelve  branches  and  the  labeling 
of  those  branches  is  shown  in  Table  III-E.  The  notation  is  the  same  as  that  used 
by  McGregor  et  3I.  (1973),  and  there  are  additional  diagrams  in  that  report  ex¬ 
plaining  the  notation.  The  following  line  pairs  are  overlapped  due  to  the  small 
spin  splitting  in  the  upper  state:  ?22  and  Ql2;  Qji  and  ?21’  Q 22  and  Rl2’  and  *11 
and  Q23* 


TABLE  III-E 


Notation  for  Transitions  of  NO 


Transition 

Rotational 
Upper  State 

Rotational 

Lower  State 

&  j». 

Pl2<J"> 

* 

Fj(J"-l)  - 

F”  (J") 

2d 

-i 

PU(J") 

m 

F](J"-1)  - 

F"  (J") 
lc 

-i 

p22(j") 

m 

F^(J"-1)  - 

F”c(JM) 

-i 

P21(J") 

= 

F2(J"-1)  - 

F"  <J") 

1C 

-i 

Ql2(J”) 

- 

FJ(J") 

F2C  (J") 

0 

QU(J") 

- 

F'(J") 

F"  (J") 

Id 

0 

q22(j") 

X 

F*(J") 

F*2d ( J") 

0 

Q21(J,,) 

m 

f2  ( J") 

F’.’  (J") 
lc 

0 

r12(j") 

- 

F[(J"+1) 

F"  (J") 

2d 

+1 

R11(J") 

a 

F^(J"+1)  - 

F"  (J") 
lc 

+1 

r22(j") 

- 

F^(J"+1) 

F"  (J") 

2c 

+1 

r2i(J") 

- 

F^(J"+1)  - 

F"  (J") 

Id 

+1 

+the  newer  notation  replaces  K  with  N. 


C  Experimental  Results  -  Spectroscopic  Measurements 
II1.C.1.  Determination  of  Broadening  Parameters 

Ill.C.l.a.  Bac_k£round_ 

In  the  discussion  of  the  development  of  transmission  formulas,  it  was  shovm 
that  all  quantities  necessary  to  calculate ‘ the  transmission  have  been  determined 
except  for  the  broadening  parameter  a'  and  the  oscillator  strength  *v'v"*  Reas¬ 
onable  literature  values  exist  for  the  oscillator  strength,  but  there  is  no 
concensus  for  the  broadening  parameters. 

There  are  several  techniques  for  studying  pressure  broadening.  When  the 
rotational  lines  are  closely  spaced,  as  is  the  case  for  the  NO  y  (0,0)  band, 
direct  line  shape  measurements  are  difficult.  A  curve  of  growth  technique  may 
be  used  at  various  pressures  and,  with  a  knowledge  of  the  oscillator  strength, 
the  broadening  may  be  determined.  This  technique  was  used  by  Thorson  and  Badger 
(1957)  to  determine  a  collision  diameter  for  NO  broadening  in  the  y  (0,0)  band 
by  N2.  A  similar  analysis  was  performed  more  recently  by  Tajime  et  al.  (1978), 
but  serious  deficiencies  in  their  theoretical  model  were  pointed  out  by  Dodge  and 
Dusek  (1979).  Davis  et  al.  (1976)  utilized  a  narrow-line  absorption  technique 
to  study  NO  broadening  by  N2  but,  again,  serious  problems  in  their  theoretical 
model  were  pointed  out  by  Dodge  and  Dusek  (1978)  .  Hadeishi  et  al.  (1976)  used 
a  Zeeman  tuned  Cd  lamp  to  determine  the  profile  of  a  single  NO  line  in  the  y 
(1,0)  band.  (They  report  measurements  of  the  0-1  vibrational  band  of  NO  at 
214. 438  nm,  but  the  y  (0,1)  band  is  at  about  236  ram,  so  we  assume  the  measure¬ 
ments  were  for  the  y  (1,0)  band.)  The  results  of  Tajime  et  al.  (1978)  and 
Davis  et  al.  (1976)  are  in  question  because  of  the  theoretical  problems,  but  the 
measurements  by  Thorson  and  Badger  (1957)  and  Hadeishi  et  al.  (1976)  may  be  in 
reasonable  agreement.  Unfortunately  (or  significantly),  the  results  reported 
here  are  substantially  different  than  in  those  two  papers. 

A  number  of  measurements  of  the  NO  y  (0,0)  band  oscillator  strength  have 
been  made,  and,  in  contrast  with  the  broadening  measurements,  there  is  reasonable 
agreement  for  most  of  the  literature  values.  The  results  reported  here  are 
similar  to  published  values.  The  problems  associated  with  different  techniques 
used  to  measure  oscillator  strengths  are  discussed  by  Thome  (1974) . 

Ill .  C .  1  .b .  J?rocedur£ 

The  technique  used  in  this  study  for  the  oscillator  strength  and  broadening 
measurements  of  the  NO  y  (0,0)  band  was  unique  in  that  the  rotational  structure 
of  the  band  was  resolved  and  actual  individual  line  shapes  and  strengths  were 
determined  photometrically.  This  requires  a  resolution  of  about  130,000  at 
225  nm.  With  the  additional  ability  to  scan  a  spectrum  with  photoelectric  de¬ 
tection,  it  is  possible  to  avoid  the  nonlinear  response  of  photographic  plates. 
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Commercial  instruments  with  these  specifications  have  only  recently  become 
available.  In  addition,  a  sophisticated  computer  model  developed  jointly  at 
Arnold  Engineering  Development  Center  (AEDC)  and  this  laboratory  was  required 
for  reducing  the  data. 

The  individual  line  profiles  were  recorded  mostly  in  the  spectral  region 
between  the  and  P22  bandheads  where  the  lines  were  separated  most  clearly. 

All  of  these  lines  originated  in  the  ^r2j2  8rounc*  state.  A  few  measurements 
were  recorded  in  the  Pjj  bandhead  region  to  examine  lines  originating  in  the 
in  the  ^^1/2  ground  state. 

All  room  temperature  broadening  measurements  were  made  with  the  high  pressure 
Xe  DC  arc  lamp  (continuum  source) ,  the  static  cell  and  associated  gas  mixing 
apparatus,  and  the  high  resolution  spectrometer  which  were  previously  described. 
The  ratiometric  electronics  were  used  for  recording  some  spectra  while  others 
were  recorded  directly  (Fig.  18). 

The  cell  was  first  filled  with  about  5  torr  (6.7  Pa)  of  10%  NO  diluted  in  Ar 
and  a  scan  was  taken  to  determine  the  slit  function  of  the  instrument,  along  with 
a  small  contribution  from  the  Doppler  broadened  gas.  The  slit  function  varied 
slightly  in  shape  and  width  over  several  hours  while  data  were  being  recorded, 
but  the  best  fit  to  the  experimental  values  could  be  obtained  with  a  Gaussian. 

Our  slit  function  was  slightly  wider  near  the  base  than  a  Gaussian,  but  it  did 
not  possess  the  wings  characteristic  of  a  Lorentzian.  A  more  complicated  slit 
function  such  as  given  by  Kusch  et  al.  (1977)  did  not  seem  warranted.  The  effect 
of  convolving  a  slit  function  with  the  actual  spectral  shape  is  discussed  by 
Thorne  (1974),  Kusch  et  al.  (1977),  and  Sulzmann  et  al  (1976).  The  convolution 
technique  used  here  is  described  in  the  computer  program  description  (Appendix  D) 
The  measured  FWHM  for  the  slit  function  was  typically  0.0018  nm  while  the  Doppler 
FWHM  was  0.00051  nm.  This  slit  function  was  monitored  several^times  during 
each  data  set.  For  Gaussians,  the  half  widths  combine  as  (Av)  «  (Av^)^  + 

(Av2)2  +  .  .  .  (Thorne  (1974)),  and  since  both  the  slit  function  and  Doppler 
shape  are  Gaussians,  the  actual  slit  function  FWHM  was  typically  .00173  nm.  In 
comparison,  the  range  of  FWHM  for  N2  for  which  broadening  data  were  reduced  was 
0.00562  to  0.00965  nm,  so  slight  variations  in  the  slit  function  (^0.0002  nm) 
had  negligible  effect  on  the  measured  profiles.  An  example  spectrum  at  low  pres¬ 
sures  is  shown  in  Fig.  23. 

After  the  slit  function  was  determined,  the  gases  were  blended  to  arrive  at  a 
constant  NO  number  density  for  a  given  final  pressure,  independent  of  diluent  gas. 
Spectra  were  recorded  at  total  pressures  of  0.50  atm  (50.7  kPa),  0.75  atm  (76.0 
kPa) ,  1.00  atm  (101  kPa),  1.50  atm  (152  kPa) ,  and  2.00  atm  (203  kPa) .  Experi¬ 
mentally  recorded  spectra  for  NO  broadening  by  N2  over  this  range  of  pressures  are 
shown  in  Figs.  24-28.  Similar  spectra  were  obtained  for  NO  diluted  in  CO2, 

CO,  Ar,  and  CH4.  These  experimental  spectra  were  compared  with  computer  generated 
spectra  such  as  those  shown  in  Figs.  29-31,  which  correspond  to  Nj  broadening  at 
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pressures  of  0.75  atm  (76.0  kPa) ,  1.00  atm  (101  kPa),  and  1.50  atm  (152  kPa) . 

The  spectral  lines  in  Figs.  24-28  are  identified  in  Table  III-F. 

At  these  NO  densities,  there  is  significant  absorption  between  even  the  most 
widely  spaced  lines  over  most  of  the  pressure  range  studied.  Thus  the  zero  ab¬ 
sorption  baseline  must  be  determined.  Even  with  the  dual-beam  ratiometric 
electronics  small  amounts  of  drift  occurred  in  the  zero  absorption  signal.  The 
baseline  was  determined  in  two  ways.  First,  the  zero  absorption  signals  before 
and  after  the  runs  were  averaged  and  called  the  estimated  baseline  (EST.  BASE,  in 
Figs.  24-28).  Second,  for  the  higher  pressure  cases,  the  computer  spectra  were 
used  to  estimate  a  baseline  and  called  computer  baseline  (COMPUTER  BASE,  in  Figs. 
24-28).  These  agreed  very  well  for  the  N2  and  Ar  data.  The  CO2  and  CH^  generally 
showed  somewhat  more  absorption  in  the  wings  than  predicted  by  the  computer  model, 
and  thus,  the  estimated  baseline  was  as  much  as  5  percent  higher  than  the  computer 
baseline  relative  to  zero  transmission  for  the  worst  case  at  2  atm  (203  kPa)  con¬ 
dition.  Such  a  comparison  for  CO  could  not  be  made  easily  because  pure  CO  gas 
absorbed  significantly  in  the  region  of  interest  as  shown  in  Fig.  32.  None  of  the 
other  gases  showed  any  absorption  in  this  region.  All  data  were  reduced  relative 
to  the  computer  baseline. 

Reported  data  are  only  for  the  Q22  (9.5)  +  (9.5)  line  pair  and  the  Q22 

(10.5)  +  (10.5)  line  pair  and  for  the  pressure  range  of  1.00  atm  (101  kPa) 

to  2.00  atm  (203  kPa) .  This  pressure  was  selected  to  reduce  uncertainties  from 
changes  in  the  slit  function.  However,  at  the  lower  pressures  where  half  widths 
could  be  measured  for  all  lines,  the  widths  were  consistent,  and  independent 
of  J  value  within  experimental  precision.  At  all  pressures  there  was  excellent 
agreement  between  the  computer  predicted  spectral  shape  and  the  observed  profiles. 
The  same  pressure  broadening  value  gave  an  excellent  match  between  experimental 
and  computed  spectra  for  transitions  originating  in  the  ^1/2  state  for  the 
pressures  0.75  atm  (76.0  kPa) ,  1.00  atm  (101  kPa) ,  and  1.32  atm  (133  kPa) ,  as 
shown  in  Figs.  33  and  34. 

III.C.l.c  Low  Temperature^  Broaden in g_D  at £  and  Di£cu_ss^i£n_ 

A  listing  of  the  typical  experimental  results  for  the  broadening  of  NO  by  N£ , 
C02,  CO,  CH^,  and  Ar  is  given  in  Table  III-G.  This  table  gives  an  indication  of 
scatter  in  the  data.  All  error  bands  in  Table  III-G  are  lo  of  the  experimental 
data  and  do  not  reflect  systematic  errors.  It  is  estimated  that  the  broadening 
parameters  C  and  K  are  accurate  to  within  15  percent. 

As  discussed  in  the  section  on  broadening  theory,  the  Lorentz  theory  gives  the 
temperature  dependence  of  a'  as  a'  =  CP/T,  while  the  Weisskopf  theory  results  in 
a'  -  KP/T*'^  (Thorne  (1974)).  The  optical  collision  diameter  d  is  also  shown. 

As  discussed  previously,  collision  cross-sections  are  defined  as  either  d or 
Trd2. 
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TABLE  III-F 


SPECTRAL  LINES  USED  FOR  BROADENING  MEASUREMENTS 


Group 

Number 

Line 

Identification 

Wavelength* 

(nm)  or  (X/10) 

Wavenumber 

(cm"*) 

1 

P22(15.5) 

Q12(15.5) 

226.66669 

226.66647 

44117.642 

44117.685 

2 

Q22(8.5) 

R12(8.5) 

P12<26.5) 

226.65247 

226.65234 

226.65156 

44120.410 

44120.436 

44120.587 

3 

R22 (4.5) 

226.64758 

44121.363 

4 

p22(16 • 5) 
Q12(16.5) 

226.63137 

226.63113 

44124.519 

44124.564 

5 

Q22 (9.5) 

R12(9.5) 

226.61623 

226.61608 

44127.466 

44127.496 

6 

R22(5.5) 

226.60194 

44130.248 

7 

P22(17.5) 

Q22(17.5) 

226.59325 

226.59300 

44131.941 

44131.989 

8 

Q22(10.5) 

R12(10.5) 

226.57724 

226.57708 

44135.060 

44135.091 

+TheoreticaIly  determined  (difference  in  upper  and  lower  energies) 
These  are  slightly  different  if  experimental  line  locations  by 
Engleman  et  al.  (1970)  or  Deezsi  (1958)  are  used. 
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TABLE  III-G 

BROADENING  PARAMETERS  FOR  NO  y(0,0) 

(A2I-X2tt)  TRANSITIONS  (ALL  TOLERANCES  ARE  +  lo  OF  MEASURED  DATA) 


Gss 

Pressure 

FWHM 

Slit 

Function 

(atm) 

(nm) 

(nm) 

n2 

1.00 

.00562 

.00176 

n2 

1.50 

.00743 

n2 

2.00 

.00965 

n2 

co2 

AVG 

1.00 

.00565 

.00191 

CO  2 

1.50 

.00720 

CO  2 

2.00 

+.00040 

.00961 

co2 

CO 

AVG. 

1.00 

+.00035 

.00580 

.00148 

CO 

1.50 

+.00013 

.00747 

CO 

2.00 

+.00050 

.01038 

CO 

ch4 

AVG. 

1.00 

+.00010 

.00635 

.00160 

ch4 

1.50 

+.00033 

•C0836 

ch4 

2.00 

+.00013 

.01041 

ch4 

Ar 

AVG 

1.00 

.00490 

.00181 

Ar 

1.50 

+.00029 

.00638 

Ar 

2.00 

+.00016 

.00803 

Ar 

AVG 

+.00040 

a' 

C 

K 

d 

atm"l  K 

atm-1  K1,2 

(nm) 

6.04 

1800 

8.50 

1689 

11.31 

1685 

1725 

5383 

1.13 

5.75 

1705 

8.12 

1605 

11.17 

1656 

1655 

5165 

1.17 

6.02 

1782 

8.47 

1671 

12.10 

1791 

1748 

5455 

1.15 

6.90 

2046 

9.63 

1904 

10.41 

1803 

1918 

5986 

1.10 

4.25 

1258 

6.95 

1371 

9.22 

1365 

1331 

4154 

1.04 

II 1-44 


There  are  two  items  of  note  in  Table  III-G.  First,  the  collision  para¬ 
meters  for  the  gases  shown  are  very  nearly  the  same  with  the  exception  of 
argon  which  is  only  a  minor  constituent  of  air  fed  combustion  except  in  some 
research- type  flames.  This  implies  that  the  calibration  for  NO  absorption 
will  be  relatively  unaffected  for  various  proportions  of  the  gases  shown 
assuming  the  gases  have  the  same  relative  collision  efficiency  at  elevated 
temperatures.  Absent  from  Table  III-G  are  two  major  species,  H2O  and  02- 
Water  vapor  measurements  could  not  be  made  with  a  room  temperature  cell, 
although  data  were  extracted  from  broadening  measurements  in  H2/02/Ar  flames 
as  will  be  discussed.  An  effort  was  made  to  measure  broadening  by  O2,  but 
the  oxidation  of  NO  to  NO-?  in  the  presence  of  strong  ultraviolet  light  was 
too  rapid  to  make  precise  measurements.  However,  it  was  estimated  from  the 
recorded  spectra  that  the  broadening  parameter  for  0^  was  not  significantly 
different  from  the  other  gases  tested  and  that  its  value  was  between  that 
of  Ar  and  N2.  Because  of  the  similarity  in  molecular  size,  calculations  for 
flames  were  made  assuming  that  N2  and  O2  have  the  same  collision  parameter. 

The  convenience  of  a  calibration  almost  independent  of  the  gas  composition  should 
be  contrasted  with  the  extreme  dependence  shown  for  NO  measurements  by  fluore¬ 
scence  (Schwarz,  1975).  it  should  be  noted  that  the  data  in  Table  III-G  are 
for  a  constant  NO  density  independent  of  pressure,  so  anomalously  large  self¬ 
broadening  by  MO  cannot  explain  the  relative  independence  of  molecular  type  on 
the  broadening  parameter. 

The  magnitude  of  the  collision  cross-sections  shown  in  Table  III-G  are  much 
larger  than  previously  reported,  and  may  revive  an  argument  that  has  existed  in 
literature  for  some  time.  It  was  reported  by  Lambrey  (1929,  1930)  and  by  Naude 
(1930)  that  the  y  system  of  NO  showed  abnormally  large  pressure  broadening.  An 
explanation  offered  by  Moore,  Wulf,  and  Badger  (1953)  was  that. this  might  be  due 
to  induced  predisssociation .  The  suggestion  of  abnormally  large  broadening  was 
opposed  by  Gaydon  and  Fairbairn  (1954)  and  Thorson  and  Badger  (1957) .  All  of 
these  measurements  were  made  with  photographic  plates,  which  present  a  much  more 
difficult  problem  for  reducing  line  widths  than  the  photometric  data  recorded  in 
this  study.  It  is  not  clear  what  constitutes  abnormally  large  broadening, 
but  our  results  are  in  conflict  with  those  of  the  latter  two  references.  For 
a  partial  pressure  of  NO  (PN0)  of  2  torr  (267  Pa)  and  a  path  length  ££)  of  15  cm, 
and  a  total  pressure  with  N2  of  1  atm  (101  kPa) ,  Gaydon  and  Fairbairn  (1954) 
estimate  the  true  half-breadth  or  -width  of  a  line  is  less  than  .025  A  (.0025 
nm) .  In  this  study,  the  corresponding  value  is  .00562  nm  including  the  slit 
function  and  ~.0050  nm  after  deconvolving  the  slit  function.  Here,  the  optical 
depth  was  less,  i.e.,  PN0  =  .85  torr  (113  Pa)  and  *■  =  18.6  cm.  Since  the  total 
pressure  was  the  same,  narrower  lines  rather  than  broader  lines  should  have  been 
observed.  Similarly,  collision  cross-sections  of  this  study  are  significantly 
greater  than  those  of  Thorson  and  Badger  (1957)  for  y  (0,0)  lines  and  those  of 
Hadeishi  et  al.  (1976)  for  y  (1,0)  lines  as  shown  in  Table  III-H.  Included  in 
the  table  are  collision  diameters  from  viscosity  data,  although  it  is  not 
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unusual  for  optical  collision  diameters  to  be  larger  than  those  from  viscosity 
data  (Engleman  (1969),  Townes  and  Schawlow  (1955),  Mitchell  and  Zemansky  (1971)). 

1 1 1 .  C .  1 .  d  ,  Hi£h_Tem£era_tu_re_B_roadejrin£  j)at_a_an.d_Dis£ussi_on_ 

The  same  procedure  was  used  to  measure  the  actual  line  profiles  in  the 
^2^2/Ar  flames.  However,  the  relative  precision  of  these  measurements  was  less 
than  that  for  the  room  temperature  static  cell  data  because  the  line  widths  were 
much  narrower.  A  FWHM  of  about  0.0031  nm  was  typical  with  a  slit  function  FWHM 
of  about  0.0018  nm.  Thus,  the  actual  line  width  without  the  instrument  function 
was  about  0.0024  nm  which  is  not  significantly  larger  than  the  slit  function. 

The  measured  a'  values  are  shown  in  Fig.  35  for  four  H2/02/Ar  flames.  Data  are 
also  shown  for  three  CH4/O2/N2  flames,  which  are  quite  similar.  The  solid  line 
is  for  an  a'  given  by  the  Weisskopf  theory  and  the  dotted  line  for  the  Lorentz 
theory.  However,  values  of  the  broadening  parameters  C  and  K  were  estimated  for 
both  H2O  and  O2.  Oxygen  was  assumed  to  be  as  efficient  as  N2  as  a  broadener, 
which  is  probably  reasonable.  Water  vapor  is  not  similar  to  any  of  the  molecules 
tested,  but  was  assumed  to  be  an  efficient  broadener  slightly  stronger  than  CH4, 
which  has  a  similar  molecular  weight.  For  these  assumptions,  the  Weiskopf  model 
fit  the  measured  data  within  experimental  precision.  Perhaps  more  significantly, 
actual  line  profiles  were  determined  at  flame  temperatures  independent  of  any 
extrapolations  from  room  temperature,  albeit  with  some  uncertainty. 

III.C.2.  Determination  of  Oscillator  Strength 

III. C. 2. a.  Continuum  Lamp,  Procedure  and  Results 


The  oscillator  strength  was  determined  from  the  same  experimental  data  used 
for  the  broadening  measurements.  The  procedure  was  to  measure  the  peak  absorp¬ 
tions  on  the  same  line  pairs  over  the  same  pressure  range  and  to  compare  that 
data  with  results  from  the  computer  model.  Peak  heights  rather  than  areas  were 
used  to  reduce  the  errors  due  to  baseline  uncertainty,  and  to  simplify  data 
reduction.  For  the  region  between  the  P^  and  P22  bandheads  used  for  the 
broadening  measurements,  all  absorption  lines  are  due  to  transitions  connected  with 
the  ^3/2  level.  Since  Spindler  et  al.  (1970)  suggested  that  the  strengths  of  the 
transitions  connected  with  the  ^3/2  state  are  about  6  percent  larger  than  those  for 
the  ^3/2  state,  separate  measurement  of  the  oscillator  strength  was  made  near 
the  P;q  bandhead  region  for  transitions  originating  in  the  ^3/2  state,  as  shown 
in  Fig  33.  The  results  are  given  in  Table  III-I.  Data  are  shown  for  NO 
mixtures  in  various  diluent  gases,  but,  of  course,  the  oscillator  strength  must 
be  independent  of  the  diluent  gas.  Since  the  NO/N2  mixture  was  confirmed  indep¬ 
endently  of  the  gas  vendor,  the  oscillator  strength  values  for  that  mixture  were 
used.  The  uncertainty  in  these  values  is  estimated  to  be  +  10  percent. 
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MEASURED  NO  BROADENING  PARAMETER  IN  FLAMES 


FIG.  35 
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TABLE  III-I 


OSCILLATOR  STRENGTHS  FOR  THE  NO  y(0,0)  (A2I+-X2ti)  BAND 
FOR  DIFFERENT  GASES* 


NO  Diluted  In 

A2£+-X2tt3/2 

»2V+  v2„ 

A  I  -X 

n2 

3.57  +  0.06+ 

3.72  +  .23 

co2 

3.54  +  0.004 

ch4 

3.51  +  0.02 

Ar 

3.27  +  0.06 

*The  oscillator  strength  must  be  independent  of  diluent  gas,  so  the 
variation  in  values  is  indicative  of  the  experimental  error. 
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These  results  are  -.red  with  some  previous  values  in  Table  1I1-J.  The 

value  for  the  ~"\/Z  state  was  chosen  for  comparison  as  it  is  the  more  heavily 
populated  state  (bv  a  factor  of  1.8)  at  room  temperature.  The  oscillator 
strength  determined  in  this  study  agrees  well  with  the  previous  results,  but  appears 
|  1°  he  slight  iv  level  limn  tin-  most  reliable  values.  However,  this  studv  is 

|  unique  in  that  the  oscillator  strength  lias  been  measured  directly  from  individual 

lines  in  high  resolution.  We  are  not  familiar  with  any  previous  work  which 
:  provides  such  a  direct  measure  of  individual  line  strengths. 

IIl.C.2.b.  Narrow-line  Lamp,  Procedure  and  Results^ 

The  broadening  data  determined  by  the  continuum  lamp  were  used  in  the  model 
to  predict  transmission  by  the  narrow-line  lamp.  The  oscillator  strengths  were 
adjusted  from  those  determined  by  the  continuum  lamp  in  order  to  match  experimental 
data  obtained  with  Che  narrow-line  lamp.  The  oscillator  strength  should  be  same 
for  both  models,  but  several  sources  of  error  can  cause  small  differences. 

The  narrow-lino  lamp  contains  a  small  amount  of  radiation  not  accounted 
for  in  the  model.  This  can  be  seen  by  comparing  Figs.  21  and  22  which  display 
output  from  the  narrow-line  lamp  before  and  after  absorption  by  a  low  pressure 
NO  sample.  The  major  peaks  denoted  1-5  show  absorption  as  predicted  by  the 
model,  but  the  low  level  light  between  the  peaks  is  not  absorbed  by  the  NO.  This 
tends  to  decrease  the  f  value  required  to  fit  experimental  data.  An  estimate  of 
the  effect  of  the  "extra-light"  on  f  value  is  shown  in  Table  ITI-K  along  with  the 
values  determined  for  the  narrow-line  lamp  before  correction.  There  may  be  other 
factors  associated  with  this  difference  Lr:  f  values,  such  as  a  shift  in  frequency 
for  the  emission  lines  due  to  the  strong  electric  fields  in  the  lamp,  hut  these 
have  not  been  investigated. 

III.C.3  Summary  of  Hot  Calibration  Device  Data 


1 1 1.  C.  3.  a  Fl£w_in£  fas  Hett  t  e  r  j _ Summary  _of_  Concent  rat_ion  anti  Temp  <2  r  <a  t  ure  Dat_a 

The  carrier  flow  rate  through  the  flowing  gas  heater  was  held  constant 
at  0.109  moles/sec  for  all  temperatures  and  gases.  Typical  concentration  and 
temperature  profiles  are  depicted  in  Figs.  8-10  and  36.  Various  concentrations 
of  nitric  oxide  were  added  to  the  main  gas  flow  at  each  temperature,  and  for 
each  condition  gas  samples  were  extracted  from  the  center  of  the  optical  axis 
and  analyzed  using  the  chemiluminescence  detector.  Although  measured  values 
were  nearly  always  lower  than  the  calculated  seed  level,  the  two  values  were 
usually  within  7%.  Predicted  nitric  oxide  concentrations  from  the  optical 
measurements  are  compared  to  the  average  of  the  seid  and  measured  values. 
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TABLE  III-K 


OSCILLATOR  STRENGTHS  (fQ  Q)  FOR  NO  y(0,0)  BAND: 
A  COMPARISON 


fn  n  x  10* 

Transition 

Continuum 

Lamp 

Narrow-Line 

Lamp 

Narrow-Line  Lamp 
With  "Extra-Light" 
Correction 

2y  2n 

Z-  tt1/2 

3.72 

3.53 

3.76 

?  2 

E_  *3/2 

3.57 

3.20 

3.43 

HI-52 
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FIG.  36 
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1 1 1 .  C .  3  •  b  Flai t_Flame_Burne£:  _  Summary  £f  _Con£e£t  ra £i£n_and_Tem£e ra tu  re_Dat£ 

For  the  optical  measurements  of  NO,  four  flame  conditions  were  selected. 

The  thermocouple  temperatures  (uncorrected  for  radiation)  at  the  centerline 
were  950,  1220,  1400  ana  1600  K.  The  flow  conditions  for  these  flames  are 
listed  in  Table  II-B.  Profiles  of  thermocouple  temperatures  and  radiation 
corrections  are  shown  in  Figs.  17  and  37-39.  Various  NO  concentrations  were 
introduced  into  the  unburned  gas  for  all  four  flame  conditions.  Concentration 
profiles  (normalized)  are  shown  in  Fig.  16.  Except  for  a  few  high  seed  levels 
of  NO  during  the  gas  correlation  experiments,  all  nitric  oxide  was  introduced 
using  the  10%  NO  and  90%  Ar  mixtures.  The  excess  argon  associated  with  the 
added  NO  is  not  included  in  the  tables  and  typically  amounts  to  1  to  2%  of  the 
total  gas  flow.  For  several  of  the  gas  correlation  experiments,  the  25%  NO, 

75%  Ar  mixtures  were  used  to  reduce  the  effect  of  the  added  argon.  Once  con¬ 
servation  of  nitric  oxide  had  been  confirmed  using  mass  spectrometric  analysis, 
it  was  not  necessary  to  repeat  these  measurements  for  all  tests.  During  the 
UTRC  optical  measurements,  therefore,  only  cold  flow  measurements  of  NO  in  Ar 
were  made  before  and  after  each  set  of  experiments.  These  measurements  were 
made  to  verify  orifice  calibrations  and  to  ensure  that  leaks  or  other  problems 
did  not  develop  in  the  gas  handling  system.  Concentrations  of  nitric  oxide 
during  flame  conditions  were  determined  using  calculated  flow  rates  and  assuming 
complete  combustion  to  water.  Since  the  total  concentration  of  other  species 
such  as  H2,  H2O2  and  radicals  is  rather  small,  this  later  assumption  provides 
a  good  and  simple  approximation.  According  to  techniques  outlined  in 
Section  II.B.3.d,  the  mass  spectrometric  data  were  reduced  to  provide  the 
profiles  of  other  species  within  the  flame.  These  data  are  reproduced  in 
Figs.  40  and  41  for  two  flame  conditions. 

Profiles  of  nitric  oxide  and  optical  measurements  over  the  burner  assembly 
were  also  made  at  room  temperature  for  comparison  with  the  FGH  and  the  static 
cell.  The  profile  data  is  reproduced  in  Fig.  42. 

III.C.4  Narrow-Line  Calibration  Study 

As  previously  mentioned,  one  of  the  principal  goals  of  the  first  phase  of 
this  contract  was  the  development  of  a  model  for  predicting  the  transmission 
of  narrow-line  NO  radiation  through  a  gas  containing  NO  molecules  at  temperature 
range  from  300K  to  2000K.  The  validity  of  the  model  would  be  certainly  estab¬ 
lished  if  an  agreement  of  approximately  *  10%  in  NO  concentration  between  theo¬ 
retical  prediction  and  measurement  was  achieved,  The  bulk  of  the  data  presented 
in  this  section  are  within  the  specified  tolerance.  Those  data  that  are  not 
within  this  tolerance  reflect  some  of  the  experimental  difficulties  associated 
with  NO  chemistry  and  measurement. 
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1 1 1 .  C .  4 .  a  St  a_t  ic_Ce_l  1_  Opt_ica  1  M  easur emends 


A  number  of  tests  were  performed  in  an  18.6  cm  long  static  cell  at  room 
temperature.  Tests  were  performed  in  high  resolution  on  single  line  pairs 
and  in  low  resolution  (.146  nm  FWHM)  using  gases  blended  with  the  gas  mixing 
system  and  gas  mixtures  obtained  from  vendors.  Data  were  obtained  for  a  pres¬ 
sure  range  of  .0627  atm  (6.35  kPa)  to  1.32  atm  (134  kPa) .  Some  of  these  data 
for  low  resolution  scans  with  gases  (NO  diluted  in  N?)  checked  against  National 
Bureau  of  Standards  reference  gases  and  used  without  dilution  are  shown  in 
Table  III-L.  The  table  lists  the  cylinder  concentration  (NO  in  ^  by  volume), 
total  pressure,  optical  depth,  and  the  measured  transmission  for  the  two  band- 
heads  which  are  displayed  for  a  spectra  at  this  resolution  in  Fig.  D-8.  The 
first  bandhead  is  located  at  about  226.8  nm  and  the  second  is  located  at 
about  226.15  nm.  Also  given  in  Table  III-L  is  the  transmission  as  predicted 
by  the  computer  model,  and  an  estimate  of  the  error  that  would  be  incurred  in 
predicting  the  NO  concentration  in  the  static  cell  from  the  measured  trans¬ 
mission  and  the  computer  model.  Since  the  NO  concentration  is  approximately 
proportional  to  In  T,  where  T  is  the  transmission,  the  error  E  is  defined  by 
F.  = 1  -  (In  Tpggp/ln  TMEAS  ),  which  is  multiplied  by  100  to  convert  a  percent 
error.  Thus,  if  the  known  NO  concentration  is  multiplied  by  the  factor  (1  +  E) , 
the  predicted  and  measured  transmission  values  will  agree.  This  same  definition 
is  used  for  the  flowing  gas  heater  and  flat  flame  burner,  except  than  the  factor 
(1  +  E)  is  multiplied  by  the  probe  determined  NO  profile  for  agreement  between 
predicted  and  measured  transmissions. 

These  tests  were  used  to  obtain  an  estimate  of  the  apparent  line  width  of 
the  narrow-line  source,  and  the  effective  oscillator  strength.  The  line  width  in 
the  source  could  not  be  measured  directly,  but  could  be  inferred,  from  measurements 
in  the  static  cell  over  a  range  of  pressures.  At  the  higher  pressures,  the  ab¬ 
sorption  is  almost  independent  of  the  source  line  width,  and  the  absorption  data 
at  these  pressures  were  used  to  determine  the  effective  oscillator  strength  when 
using  the  narrow-line  lamp.  The  oscillator  strength  determined  in  this  manner 
resulted  in  a  slightly  smaller  value  than  that  determined  with  the  continuum  source 
as  shown  in  Table  III-K,  with  the  difference  attributable  (at  least  partly)  to  the 
"extra  light"  in  the  narrow-line  lamp,  as  discussed  previously.  At  the  lower 
pressures,  the  emission  line  width  from  the  lamp  becomes  comparable  with  the  ab¬ 
sorption  line  width,  and  the  transmittance  at  these  pressures  depends  on  the 
value  assumed  for  the  source  line  width.  The  best  fit  to  the  static  cell  data 
were  obtained  for  a  source  Doppler  width  with  a  temperature  of  600  K.  This 
compares  favorably  with  the  value  of  561  K  for  the  lower  temperature  part  of  the 
Boltzmann  plot  for  the  line  strengths  in  the  narrow-line  lamp  (Fig.  20).  Meinel 
(1975)  has  suggested  that  this  low  temperature  part  of  the  Boltzmann  plot  should 
give  a  measure  of  the  translational  temperature  of  the  lamp,  which,  in  turn, 
determines  the  Doppler  width  of  the  emitted  lines.  This  agreement  may  be 
fortuitous,  because  the  "extra  light"  phenomena  should  also  affect  the  apparent 
lamp  line  width. 


III-61 


TABLE  III-L 


Static  Cell  Calibration  Data  (NO  in  N2) 


Cylinder 

Concen¬ 

Total 

Optical 

tration 

Pressure 

Depth 

(ppm) 

(atm) 

(cm-2) 

93.3 

1.315 

5.679-1016 

93.3 

1.000 

4.317-1016 

93.  3 

.751 

3. 2«0'  10 

93.  3 

.500 

2.159-1016 

93.3 

.474 

2.044-1016 

93.3 

.250 

1.079-1016 

93.3 

,125 

5.405-1015 

470 

1.001 

2.176-1017 

470 

.751 

1.633' 1017 

470 

.503 

1.093' 1017 

470 

.249 

5.426-1016 

470 

.125 

2.721-1016 

470 

.0627 

1.363-1016 

1st 

Is  t 

Error 

Bandhead 

Bandhead 

in  Pred 

Trans . 

Trans . 

NO  Cone 

Meas. 

Pred . 

(%)  + 

.896 

.890 

-6.5 

.900 

.900 

-0.2 

.912 

.909 

-3.2 

.917 

.920 

+3.1 

.918 

.921 

+4.4 

.934 

.939 

+9.2 

.956 

.959 

+5.7 

.627 

.625 

-0.7 

.646 

.633 

-4.5 

.674 

.  666 

-2.9 

.738 

.734 

-1.5 

.816 

.812 

-2,4 

.887 

.882 

-4.2 

2nd 

2nd 

Error 

Bandhead 

Bandhead 

in  Pred 

Trans. 

Trans. 

NO  Cone 

Meas. 

Pred. 

(%)  + 

.824 

.821 

-1.9 

.843 

.842 

-0.8 

.864 

.859 

-3.5 

.879 

.879 

-0.1 

.880 

.881 

+0.9 

.906 

.909 

+3.1 

.938 

.938 

+1.3 

.465 

.469 

+1.2 

.513 

.498 

-4 . 6 

.568 

.555 

-4.2 

.668 

.651 

-6.2 

.762 

.748 

-6.8 

.846 

.835 

-7.5 

+  These  errors  are  based  on  more  significant  figures  than  shown  in  this  table.  See  discus¬ 
sion  in  Section  III. C. 4. a  for  definition  of  error. 
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1 1 1 . C . 4 . b  Flowing  Gas_  Heater_(FGH) _0£t  ic a 1  JMea su r emen t£ 

Calibration  tests  were  performed  with  the  flowing  gas  heater  for  NO  diluted 
in  pure  N2  and  pure  Ar  at  a  total  pressure  of  1  atm  (101  kPa)  and  over  a  tempera¬ 
ture  range  from  295K  to  825K.  All  of  these  absorption  spectra  were  obtained  in 
low  resolution  (FWHM  -  .146  nm  »  1.46  8),  and  transmission  data  were  reduced  for 
the  first  bandhead,  located  at  about  226.8  nm  and  the  second  bandhead,  located 
at  about  226.15  nm.  A  comparison  of  the  computer  model  prediction  versus  measured 
transmission  is  shown  in  Table  III-M.  The  values  from  the  computer  model  were 
based  on  an  NO  concentration  which  was  the  average  of  the  measured  value  (chemi¬ 
luminescent  analyzer)  and  the  value  computed  from  flow  calculations  for  the 
critical  orifices.  The  ratio  of  the  measured  NO  concentration  to  that  preducted 
by  the  flow  system  was  'v  0.96,  which  was  less  than  unity  because  of  diffusion  of 
purge  gas  into  the  main  stream  and,  possibly,  calibration  errors.  The  actual 
temperature  and  NO  concentration  profiles  were  approximated  as  two  or  three  zones 
of  constant  temperature  for  computer  predicted  transmissions.  The  data  were 
reduced  assuming  Lorentz  broadening,  that  is  a1  =  CP/T.  For  comparison  with  the 
data  reduction  procedure  used  for  the  flat  flame  burner,  the  Weisskopf  dependence 
(a’  =  KP/T^*^)  could  have  been  employed  for  the  FGH.  At  room  temperature  the 
Weisskopf  theory  would  give  the  same  result  as  the  Lorentz  theory,  but,  at  the 
highest  temperatures  in  the  FGH,  the  Weisskopf  theory  would  systematically  reduce 
the  predicted  transmissions  as  shown  for  Run  43  in  Table  III-M.  At  these  tem¬ 
peratures  the  predicted  absorptance  (-In  T)  is  increased  about  11  percent,  so 
that  the  numbers  in  column  "Error  in  Pred.  NO  Cone."  would  be  reduced  by  about 
11  percent.  The  overall  errors  would  probably  be  slightly  reduced  by  this  pro¬ 
cedure. 


III.C.4.C  Flat  Flame  Burner  (H2/02/Ar/NO)  Optical  Measurements 


Optical  calibration  for  NO  absorption  was  performed  for  H2/02/Ar/NO  flames  at 
1  atm  (101  kPa)  over  a  temperature  range  from  1000K  to  1860K.  Data  were  obtained 
independently  for  the  continuum  lamp  and  the  narrow-line  lamp. 

The  continuum  lamp  data  were  used  principally  to  determine  the  line  widths  in 
the  flames,  but  the  transmission  information  was  also  available,  and  since  this 
high  resolution  technique  was  used  on  some  of  the  CH^/N^/02/NO  flames  of  Task  III, 
data  were  reduced  for  the  H2/02/Ar/NO  flames  to  determine  its  applicability.  Trans¬ 
mission  data  for  three  flames  at  1000K,  1360K  and  1610K  are  compared  with  computer 
model  predicted  transmissions,  in  Table  III-N.  These  were  reduced  based  on  NO 
centerline  levels  computed  from  flow  calculations  for  the  critical  orifices, 
which  were  checked  in  cold  flow,  and  profile  data  as  shown  in  Fig.  16.  Two  or 
three  zones  of  constant  temperature  were  assumed. 

The  narrow-line  lamp  data  were  obtained  for  four  different  flame  temperatures 
with  at  least  two  NO  concentrations  at  each  temperature.  The  computed  and  actual 
transmissions  for  the  two  bandheads  observed  in  the  low  resolution  spectra  are 
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TABLE  III-M 


Flowing  Gas  Heater  (FGH)  Optical  Measurements 


1st 

1st 

Error  in 

2nd 

2nd 

Error  in 

Run 

Temp 

NO  Cone. 

Diluent 

Bandhead 

Bandhead 

Pred.  NO 

Bandhead 

Bandhead 

Pred.  NO 

Number 

(Centerline) 

(Centerline) 

Gas 

Trans . 

Trans . 

Cone . 

Trans . 

Trans . 

Cone  . 

X 

(K) 

(ppmV) 

Meas . 

Pred . 

(%)  + 

Meas. 

Pred. 

(%) 

31 

295.8 

287 

Ar 

.735 

.752 

+7.9 

.609 

.643 

+12.3 

32 

295.8 

112 

Ar 

.882 

.895 

+13.3 

.818 

.841 

+16.3 

33 

295.8 

108 

n2 

.904 

.912 

+9.1 

.848 

.861 

+9.9 

34 

295.8 

280 

*2 

.779 

.790 

+6.2 

.661 

.684 

+9.2 

35 

4°3. 

112 

Ar 

.926 

.930 

+6.4 

.899 

.903 

+4.7 

36 

453. 

293 

Ar 

.793 

.815 

+13.3 

.726 

.750 

+11.4 

37 

524. 

111 

n2 

.945 

.944 

-2.6 

.920 

.920 

+0.4 

38 

452. 

280 

n2 

.823 

.843 

+14.4 

.764 

.783 

+10.1 

39 

825. 

110 

Ar 

.962 

.961 

-2.6 

.947 

.946 

-0.6 

40 

788. 

287 

Ar 

.891 

.896 

+5.0 

.854 

.859 

+3.5 

41 

768. 

284 

n2 

.907 

.908 

+1.9 

.869 

.873 

+3.1 

42 

793. 

109 

n2 

.962 

.965 

+11.8 

.949 

.951 

+4.9 

43 

756. 

430 

n2 

.856 

.862 

+5.2 

.805 

.811 

+3 . 6 

44 

738. 

521 

Ar 

.791 

.809 

+10.7 

.733 

.747 

+6.2 

43' 

756 

4  30 

n2 

.856 

.847 

-  6.8 

.805 

.794 

-6.1 

-  These  errors  are  based  on  more  significant  figures  than  shown  in  this  table.  See  discus¬ 
sion  in  Section  III. C. 4. a  for  definition  of  error. 


§Test  case  for  a' determined  from  Weisskopf  theory.  See  text. 
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shown  in  Table  III-O.  The  computed  transmissions  were  based  on  a  centerline  value 
which  was  the  average  of  calculated  levels  for  the  critical  orifices  and  chemi¬ 
luminescent  measurements  (corrected  for  Ar  response),  and  profile  data  as  shown 
in  Fig.  16.  Two  or  three  zones  of  constant  temperature  were  assumed.  The  a'  values 
used  to  generate  the  computer  spectra  were  taken  from  the  solid  line  in  Fig.  35. 

An  estimate  of  the  error  in  predicted  transmission  due  to  the  uncertainty  in  the 
a'  values  is  also  shown  in  Table  III-O  for  the  highest  temperature  flame.  The 
uncertainty  in  a'  of  0.6  +  0.2  results  in  a  variation  of  predicted  NO  of  +9.7 
percent . 

The  overall  agreement  between  predicted  and  measured  values  was  within  -  16 
percent  for  the  flame  data,  which  was  slightly  greater  than  the  original  goal  of 
i  10%,  but  was  reasonable  given  the  experimental  uncertainties  and  relative  to  the 
previously  reported  discrepancies.  It  should  be  noted  that  the  overall  errors 
result  from  the  summation  of  errors  in  the  model,  scatter  in  the  experimental 
spectra,  uncertainties  in  the  broadening  parameter  a',  uncertainties  in  the  flame 
profiles  for  temperature  and  NO,  and  errors  in  the  gas  blending  procedure. 
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IV.  DISCUSSION 


A.  Introduction 

The  computer  model  that  was  developed  by  W.  K.  McGregor,  M.  Davis,  J.  D.  Few  and 
their  colleagues  at  Arnold  Research  Organization  (ARO)  was  one  of  the  starting  points 
of  this  study.  Initially,  it  was  thought  that  the  model  rested  on  a  sound  theoretical 
basis  and  was  supported  by  strong  experimental  evidence.  However,  experimental  data 
obtained  with  the  static  cell,  i.e.,  under  the  most  ideal  conditions,  were  not  in 
agreement  with  the  predictions  of  the  model.  A  thorough  examination  of  the  program 
for  errors  in  translation  onto  the  UTRC  UNIVAC  1110  computer  revealed  no  such  errors. 
Additional  verification  of  proper  translation  was  obtained  when  the  results  of  test 
computations  made  at  both  UTRC  and  ARO  were  in  agreement.  An  intense  review  of  gas 
standards,  mixing  procedure,  optical  and  electronic  instrumentation  did  not  yield  the 
cause  of  the  disagreement.  Consequently,  the  theoretical  basis  of  the  model  was  studied 
and  found  deficient.  The  errors  that  were  discovered  are  reported  in  Appendix  B. 

These  errors  appear  in  reports  published  since  1973;  however,  a  precise  knowledge  of 
the  various  forms  of  the  model  from  1973  through  1977  is  not  known.  Suggested  refine¬ 
ments  to  the  model  delivered  to  UTRC  are  also  noted  in  Appendix  B.  A  new  model  was 
developed  on  the  theory  presented  in  Section  III  of  this  report.  A  description  of 
the  new  model  is  given  in  Appendix  D.  A  major  subroutine  entitled  N0*Spect  (Appendix 
D)  was  also  provided  to  ARO  in  July  1978.  This  subroutine  can  be  used  to  interpret 
only  resonant  lamp  data. 

B.  Comments  on  ARO  AEDC  TMR-79-P7 

As  part  of  this  study,  measurements  were  made  on  the  UTRC  calibration  devices  by 
ARO  personnel  using  a  capillary  discharge  lamp.  The  results  of  those  measurements  are 
contained  in  an  internal  ARO  Report  (AEDC  TMR-79-P7)  by  Few  ,  McGregor,  and  Keefer 
(1979).  This  report  will  be  published  externally  at  a  later  date. 

In  order  to  prevent  further  confusion  in  the  literature  stemming  from  the  errors 
contained  in  the  reports  listed  in  Appendix  B  and  problems  contained  in  AEDC  TMR-79-P7, 
it  seemed  necessary  and  was  requested  by  the  FAA  that  the  content  of  this  latest  ‘t 
report  be  discussed  here.  That  report  has  four  major  sections  each  of  which  will  be 
considered  below. 

IV.B.l.  TMR  Introduction 

The  last  part  of  this  section  should  be  expanded  to  include  a  statement  that  the 
agreement  of  the  optical  measurements  with  the  calculated  input  of  MO  concentration  in 
conjunction  with  the  probe  profiles  are  generally  within  the  uncertainties  of  the  probe 
and  optical  measurements.  As  written  it  could  be  inferred  that  the  uncertainties  rest 
only  in  the  probe  measurements.  Such  an  interpretation  would  not  be  accurate  as  will 
be  Indicated  later  in  this  review. 
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IV. B. 2  TMR  Review  of  Model 


In  subsection  A  which  is  on  the  spectroscopic  theory,  it  is  stated  that  the 
Honl-London  factors  are  those  of  Earls  as  previously  defined  in  a  report  by  Davis, 
McGregor,  and  Few  (AEDC-TR- 76-12,  February  1976).  Unfortunately  some  of  the  Honl- 
London  factors  included  in  that  report  are  in  error.  In  subsection  B  on  model 
parameters,  the  authors  state  that  the  oscillator  strength  value  of  4.09  ±  0.1  x 
10-4  Hasson  et  al)  appears  to  be  the  most  reliable.  The  author's  judgment  is  neither 
in  agreement  with  Huber  and  Herzberg  (1979)  who  prefer  3.8  x  10"^  nor  with  3.7  x  10“^ 
determined  by  UTRC.  With  regard  to  source  temperature,  it  should  be  remembered 
that  950  ±  25  K  is  a  factor  of  3  higher  than  that  previously  reported  (e.g.  McGregor, 
Few.  Litton  (1973),  Few  Bryson  and  McGregor  (1976)).  Moreover,  the  broadening 
parameters  reported  are  a  function  of  both  the  oscillator  strength  that  was  selected 
from  the  literature  and  the  source  temperature,  and,  b«nce,  are  inferred  and  not 
directly  determined. 


IV. B. 3  TMR  Optical  Results 

In  the  description  of  the  calibration  facility,  it  is  incorrectly  stated  that  the 
majority  of  the  sampling  measurements  on  the  flowing  gas  heater  weremadewith  the  mass 
spectrometer.  In  fact,  the  bulk  of  these  measurements  was  made  with  the  chemilumine¬ 
scent  analyzer.  It  is  also  incorrectly  stated  that  the  buffer  flame  of  the  flat  flame 
burner  was  lean  burning  ^/C^/Ar.  The  buffer  flame  was  actually  a  lean  CH^/C^/Ar 
flame.  Figures  5  and  6  give  respectively  concentration  and  temperature  profiles  over 
the  flowing  gas  heater.  In  these  figures,  the  uncertainty  is  listed  as  unknown. 

These  uncertainties  are  given  in  Figures  9  and  10  in  the  main  body  of  the  Task  I  Report 
Similarly,  Figures  7  and  8  of  the  ARO  report  give  concentration  profiles  with  "unknown" 
uncertainties.  See  Figures  42  and  16  of  the  Task  I  Report.  Also,  in  Figure  11  of  the 
ARO  report,  the  40%  uncertainty  refers  only  to  the  radiation  correction  factor  and  is 
not  the  uncertainty  of  the  temperature  measurement  itself.  It  is  not  clear  why,  as 
in  the  Introduction,  the  implication  is  given  that  the  uncertainties  in  the  experimen¬ 
tal  data  reside  only  in  probe,  thermocouple,  and  calibration  device  performance.  In¬ 
deed,  this  is  an  incomplete  assessment  of  the  uncertainties.  A  complete  assessment 
of  the  errors  should  also  include  variations  in  lamp  intensity,  drift  in  system  elec¬ 
tronics,  trigger  and  grating  scan  uncertainties  and  their  influence  on  the  signal 
averaging  process  used  by  ARO.  These  experimental  uncertainties  were  readily  identi¬ 
fiable  during  the  ARO  measurements  made  at  UTRC. 

With  regard  to  data  treatment,  the  new  model  developed  by  ARO  is  briefly  described 
This  new  model  does  not  employ  a  Beer's  Law  zonal  treatment  of  temperature  and  con¬ 
centration  ^homogeneities  encountered  in  combustion  systems.  Instead,  the  new  model 
keeps  a  detailed  knowledge  of  the  emission  line  shape  from  zone  to  zone.  Figure  17 
illustrates  in  an  exaggerated  manner  the  importance  of  a  problem  which  can  occur  when 
the  absorber  line  is  narrower  than  the  source  line.  Such  a  situation,  it  is  stated, 
would  be  encountered  for  a  jet  engine  exhaust  under  altitude  conditions.  This  state¬ 
ment  is  extremely  vague  and  in  the  context  of  this  study  is  misleading.  No  information 
is  given  on  what  altitude  such  a  condition  may  occur.  Moreover,  no  data  were  obtained 
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in  this  study  that  required  consideration  of  this  phenomena.  In  addition,  if  the  ICAO 
standard  atmosphere  is  used  to  provide  temperature  and  pressure  as  a  function  of  alti¬ 
tude  and  if  600  K  is  assumed  to  be  a  typical  static  temperature  of  a  jet  exhaust,  then 
the  absorbing  lines  at  35,000  feet  are  still  considerably  wider  than  the  emission  lines 
Even  at  65,000  feet,  significant  transmitted  line  center  "burnout"  is  not  expected. 

Equation  (17)  [Ac  =|(fn  Tmeas/^n  Tcal)“l  }•  100 1  that  was  used  to  compute  the  per¬ 
cent  error  in  concentration  given  in  Tables  3,  4,  and  5  is  not  the  same  expression 

used  in  Task  I  Renort  Tables  III-L,  III-M,  and  III-O.  The  expression  used  in  Task  1 
Report  was 


A  c 


Tcal 
£n  t 


meas 


•  100 


Using  equation  (17),  the  root-mean  square  average  departures  of  the  calculated 
NO  concentration  is  reported  in  the  TMR  as  ± 9%  for  the  absorption  cell,  - 18%  for  the 
flowing  gas  heater,  and  -21%  for  the  flat  flame  burner.  This  method  of  summarizing 
the  errors  obscures  the  fact  that  the  errors  indicated  in  Tables  4  and  5,  which  res¬ 
pectively  contain  data  from  the  flowing  gas  heater  and  flat  flame  burner  are  not  really 
random  but  systematic.  In  fact,  for  the  bulk  of  the  data  reported  in  those  Tables, 
the  error  is  positive  and,  generally,  increases  with  temperature.  The  average  error 
for  all  the  ARO  tests  from  room  temperature  to  the  highest  temperatures  was  equiv¬ 
alent  to  an  overprediction  of  NO  using  the  optical  technique  by  +  15%,  but  at  the 
highest  temperature  the  overprediction  was  +  29%.  This  systematic  error  is  most 
likely  due  to  deviations  from  Lorentz  broadening  theory.  This  deviation  is  dis¬ 
cussed  in  Section  III  of  Task  I  Report.  Moreover,  the  errors  (-15%  to  +4.8)  listed 
in  Table  III  for  the  absorption  cell  give  an  indication  of  the  errors1 associated 
with  the  optical  measurement  system  used  by  ARO. 

IV. B. 4  TMR  Summary 

With  regard  to  the  broadening  study  results,  it  is  suggested  that  the  broadening 
due  to  H2O  is  the  same  as  that  of  N2,  Ar,  and  C02-  The  basis  for  this  suggestion  is 
that  the  broadening  cross  sections  for  N2,  Ar,  and  CO2  are  similar.  This  suggestion 
cannot  be  accepted  because  of  the  evidence  presented  in  Figure  of  Task  I  Report. 
Furthermore,  H2O  is  a  polar  molecule  while  N2,  Ar,  and  CO2  are  not. 

It  is  restated  in  this  summary  that  the  ARO  procedure  avoids  the  incorrect  assump¬ 
tion  of  Beer's  law  which  is  reasonable  for  this  data  where  the  absorber  lines  are  con¬ 
siderably  wider  than  the  source  lines.  Furthermore,  serious  difficulty  with  Beer's  law 
would  be  expected  during  measurement  on  turbine  engine  exhausts  at  simulated  altitudes. 
However,  there  is  no  evidence  presented  in  this  ARO  report  that  supports  these  state¬ 
ments.  In  fact,  for  all  the  measurements  made  in  this  study,  the  Beer's  law  assump¬ 
tion  is  valid.  Moreover,  it  is  also  a  valid  assumption  for  those  altitudes  in  which 
the  majority  of  the  jet  aircraft  fleets  operate. 


Finally,  the  ARO  report  concludes  with  a  statement  that  a  20%  projected  uncet - 
tainty  can  be  expected  by  the  use  of  the  ARO  UV  resonance  line  absorption  technique 
based  on  the  results  obtained  with  uncertain  temperature  and  concentration  profiles 


in  the  calibration  devices.  It  must  be  reiterated  that  uncertainties  also  arise  from 
the  performance  of  the  optical  systems  as  indicated  by  errors  in  the  static  cell  measure¬ 
ments  which  were  obtained  under  the  most  ideal  conditions. 


C.  Comparison  of  Hollow  Cathode  and  Capillary  Discharge  Performance 

Although  these  lamps  are  very  different  in  design,  the  emission  characteristics 
are  quite  similar  as  opposed  to  the  statement  made  in  TMR-79-P7.  Overall,  the  trans¬ 
mission  characteristics  through  media  containing  NO  at  various  temperatures  can  be 
seen  to  be  comparable  from  the  data  in  Table  IV-A. 

D.  Comparison  of  the  UTRC  and  ARC  Computer  Programs 

Except  for  input  data,  e.g.,  lamp  temperature, broadening  coefficients,  etc.,  the 
main  body  of  the  present  ARO  computer  program  appears  to  be  substantially  similar  to 
the  version  developed  by  UTRC  from  the  1977  ARO  model  (Davis  et  al  (1976a)).  At  room 
temperature,  the  predictions  of  the  ARO  model  agree  with  those  of  the  UTRC  model. 
However,  the  predictions  diverge  significantly  at  elevated  temperature  due  to  the 
different  temperature  dependence  for  the  broadening  parameter  a'.  The  ARO  model  em¬ 
ploys  Lorentz  theory  while  the  UTRC  model  employs  Weisskopf  theory.  The  ARO  model 
will  systematically  predict  more  nitric  oxide  than  the  UTRC  model  for  a  given  trans¬ 
mission,  pressure,  and  elevated  temperature  distribution.  As  stated  in  Section  IV. B. 3, 
the  different  broadening  dependence  is  the  probable  cause  for  the  systematic  overpre¬ 
diction  of  NO  by  the  ARO  model  at  the  high  temperatures  in  the  flat  flame  burner  tests 
at  UTRC. 

Moreover,  it  should  be  noted  that  the  data  reported  in  Table  III-N  on  the  con¬ 
tinuum  lamp  transmissions  for  the  flat  flames  was  reduced  using  a  -model  which  stored 
the  line  shape  for  each  zone.  For  a  continuum  lamp  measurement,  this  procedure  is 
necessary.  However,  for  a  resonant  lamp,  such  a  procedure  was  not  necessary  for  this 
study.  At  typical  commercial  jet  exhaust  temperatures  and  operating  altitudes  (Hi 
40,000  ft)  this  procedure  is  not  necessary.  Moreover,  even  at  65,000  ft,  the  resonant 
line  sources  will  have  lines  which  are  narrower  than  the  exhaust  absorber  lines. 

E.  Comment  on  Gas  Correlation  Mzasurements 

As  stated  in  Appendix  A,  this  instrument  was  designed  to  measure  several  pollu¬ 
tant  molecules  emitted  from  smokestacks,  NO  being  one  of  those  molecules.  For  that 
application,  this  instrument  has  considerable  merit.  However,  for  measurements  in 
this  study,  this  instrument  in  its  present  configuration  was  not  well  suited  for 
several  reasons.  The  first  was  related  to  the  failure  of  the  detector  that  was  used 
in  the  Wright-Patterson  measurements.  The  active  area  of  the  original  detector  was 
circular  (1.3  mm  dia)  while  the  replacement  has  a  rectangular  area  (4  mm  x  0.4  mm). 

This  geometric  mismatch  is  considerable.  The  second  was  that  the  intrinsic  noise 
of  the  replacement  detector  was  higher  than  that  of  the  original  detector.  These 
two  items  alone  were  most  likely  responsible  for  a  factor  of  four  increase  in  the 
noise  equivalent  NO.  A  third  reason  was  due  to  the  lines  chosen  for  correlation. 
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For  temperatures  up  to  900K,  this  selection  of  lines  was  such  to  produce  a  temperature 
independent  calibration  curve.  Above  900K,  however,  the  selection  is  no  longer 
optimum;  hence,  the  sensitivity  decreases  by  almost  a  factor  of  three  at  the  highest 
temperatures  of  this  study. 

The  center-line  concentrations  necessary  for  adequate  signal-noise  ratios  for 
the  infrared  measurements  were  typically  factors  of  2  to  10  greater  than  those  used 
in  the  ultraviolet  measurements.  Because  of  the  loss  in  sensitivity,  a  side-by-side 
comparison  data  is  not  possible. 


I 


TABLE  IV-A 


COMPARISON  OF  TRANSMISSION  DATA 
ARO  LAMP  VERSUS  UTRC  LAMP 


TRANSMISSION 

1st  BH  2nd  BH 


ARO 

UTRC 

ARO 

UTRC 

STATIC  CELL 

93.3  ppm 

.892 

.901 

.905 

.899 

.835 

.849 

.850 

.843 

470  ppm 

.621 

.628 

.627 

.482 

.504 

.497 

FGH 

280, ppm  452  K 

.823 

.764 

288ppm,  448  K 

.817 

.753 

FFB 

999  ppm,  1360  K 

.822 

(0=1.96-10-4) 

.777 

(0=2.53- 

1066  ppm,  1360  K 

.811 

(a=1.97.10-A) 

.747 

(0=2.74-10-4) 
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V.  SUMMARY  AND  CONCLUSIONS 


Methods  for  providing  known  amounts  of  NO  from  room  temperature  to  2000  K  were 
developed.  At  room  temperature,  a  stainless  steel,  static  cell  was  used.  This  cell, 
which  was  18.6  cm  long  and  2.2  cm  in  diameter,  was  leak  checked  at  1  x  10-&  torr 
(1.33  x  10-4  Pa)  and  could  be  pressurized  to  115  psia  (792  kPa) .  Its  windows  were 
ultraviolet  grade  fused  silica.  The  cell  was  attached  to  a  precision  mixing  apparatus 
and  used  in  obtaining  the  majority  of  the  broadening  and  oscillator  strength  informa¬ 
tion  presented  in  this  report.  For  temperatures  up  to  850  K,  a  quartz-bed  heat  ex¬ 
changer,  through  which  mixtures  of  NO  in  N2  and  Ar  were  flowed,  served  as  a  calibra¬ 
tion  device.  NO  decomposition  in  the  bed  was  not  significant  (<10%).  NO  concentra¬ 
tions,  measured  with  both  uncooled  quartz  and  metallic  probes  at  and  near  the  center- 
line,  were  in  close  agreement  with  NO  concentrations  calculated  for  a  gas  mixing 
apparatus  which  employed  critical  flow  orifices.  The  NO  determinations  were  made 
with  both  chemiluminescent  analysis  and  mass  spectrometry.  For  those  cases  where  Ar 
was  the  bulk  gas,  proper  account  of  viscosity  and  quenching  phenomena  was  taken  in 
the  calibration  of  the  chemiluminescent  analyzer.  Kinetic  analysis  indicated  that 
substantial  decomposition  would  not  occur  in  this  exchanger  up  to  1000  K.  Thermal 
runaway  of  the  electrical  heating  elements,  however,  prevented  its  use  at  that  tem¬ 
perature.  For  temperatures  of  1000  K  to  2000  K,  a  lean  flame  seeded 

with  NO  was  used.  The  burner  surface  was  constructed  of  water-cooled,  sintered 
copper  with  two  zones:  a  main,  NO-seeded  zone  and  an  unseeded  buffer  zone.  Measure¬ 
ments  made  with  water-cooled  quartz  probes  indicated  downstream  conservation  of  NO 
in  flames  whose  stoichiometries  were  varied  from  0.36  to  0.92,  Ar  was  used  as  a 
bulk  gas  instead  of  N2  to  preclude  the  formation  of  thermal  NO.  The  data  confirm 
that  water-cooled  quartz  probes  can  be  used  in  measuring  NO  in  lean  H2/02/Ar  flames. 
Detailed  concentration  and  temperature  distributions  were  obtained  along  the  optical 
path  for  both  high  temperature  devices.  These  distributions  had  to  be  available  for 
the  reduction  of  the  optical  data  and  their  comparison  with  probe  results.  The 
temperature  measurements  were  made  with  thermocouples.  For  the  quartz-bed  heat  ex¬ 
changer  measurements,  uncoated  chromel-alumel  was  used.  In  the  flame  measurements, 
Ir/60%Ir-40%Rh  wires  (76um)  coated  with  10%  beryllium  oxide  and  90%  yttrium  oxide 

were  employed.  The  purpose  for  the  coating  was  to  minimize  catalysis  on  the  wires. 

Corrections  for  radiation  losses  were  applied  to  the  data  obtained  on  the  flames. 

Two  distinct  ultraviolet  sources  were  used  in  this  study.  The  first  was  a 
hollow  cathode  lamp.  A  dc  discharge  in  air  at  low  pressure  (2  Torr;  266  Pa)  pro¬ 
duced  emission  lines  mainly  from  NO  molecules,  N2  molecules  and  ions,  and  Ar  atoms. 

The  spectral  lines  used  in  this  study  were  in  the  Y(0,0),  Y(l,l)  and  Y(2,2)  bands 
(A2£+  -  X2tt)  of  NO.  The  second  source  was  a  high  pressure  Xe  lamp  which  produced 
radiation  of  the  continuum  type  through  the  above  spectral  region. 


JL. 


V-l 


r 


A  detailed  review  of  the  previously  developed  spectroscopic  theory  and  computer 
model,  which  was  supplied  to  these  authors  at  the  beginning  of  this  study  by  ARO, 

Inc.,  revealed  several  significant  errors  (Appendix  B) .  The  accuracy  of  the  optical 
results  published  by  the  originators  of  that  theory  and  model  is,  hence,  in  doubt. 

Due  to  the  complex  history  of  the  computer  model,  inconsistencies  in  calibration  data 
and  the  unavailability  of  the  original  raw  data,  it  is  not  possible  to  comment  with 
certainty  on  the  accuracy  of  the  optical  data  existing  in  the  literature.  Where 
possible,  some  reanalysis  of  the  original  data  will  be  included  in  TASK  III  Report. 

Because  of  these  difficulties,  the  spectral  theory  used  in  the  previously 
published  model  was  corrected,  and  the  computer  model  was  expanded  so  that  not  only 
data  from  resonant  line  sources  can  be  analyzed  but  also  that  from  continuum  sources 
(Appendix  D) .  Experimental  data  were  obtained  at  low  and  elevated  temperatures 
using  these  significantly  different  spectral  sources.  A  comparison  of  these  data 
showed  excellent  agreement;  hence,  the  validity  of  the  theory  and  model  were  estab¬ 
lished  . 

In  order  to  use  this  model,  a  knowledge  of  the  broadening  of  NO  spectral  lines 
in  foreign  gases  was  necessary.  Data  on  broadening  in  Ar ,  N2,  CO2,  CO  and  CH.^  were 
obtained  by  the  direct  observation  in  high  resolution  of  isolated  lines  using  the  con¬ 
tinuum  source.  If  Lorentz  theory  is  assumed,  i.e.,  a'=CP/T,  the  following  values 
of  C  were  obtained:  C(Ar)  »  1331  atm~l  K;  C(N2>  =  1725  atm-1  K;  C(C02)  =  1655  atm-1  K 
C(CO)  *  1748  atm~l  K;  and  C(CH-)  -  1918  atm~^  K.  If  Weisskopf  theory  is  assumed, 
i.e.,  a'  =  K'P/T-^’2,  then  the  same  experimental  data  yield:  K'  (Ar)  =  4154  atm-1K1,2; 
K'  (N2)  -  5383  atm-1*1-2;  K'  (C02)  =  5165  atm"3*1-2;  K’(CO)  =  5455  atm"1*1-2;  and 
K  (CH4)  =  5986  atm“^K'*”2.  These  data  were  also  reduced  to  determine  the  following 
collision  diameters:  d(Ar)  =  1.04  nm;  d(N2>  *  1.13  nm;  d(C02)  =  1.17  nm;  d(C0)  - 
1.15  nm;  and  d(CH4)  =  1.10  nm.  Because  of  the  oxidation  of  NO  to  NO2  by  02,  precise 
measurement  of  the  broadening  parameter  for  02  was  not  possible.  However,  for  the 
data  obtained,  it  was  estimated  that  the  broadening  parameter  for  02  has  a  value 
which  lies  between  that  of  Ar  and  N2*  Data  taken  on  the  flames  indicated  that  at 
elevated  temperature  Weisskopf  theory  seemed  to  provide  a  better  fit  to  the  data. 

If  that  theory  is  used  and  if  it  is  assumed  that  O2  is  as  efficient  a  broadener  as 
N2,  then  the  broadening  parameter  for  H20  is  6260  atm-3-  K3,2. 

Oscillator  strengths  were  determined  from  data  obtained  with  both  radiation 
sources.  From  the  continuum  lamp  data,  the  oscillator  strengths  for  the  NO  Y(0,0) 

A2E+  -  X2nj/2  band  when  NO  is  diluted  in  N2»  C02,  CH^,  and  Ar  were  3.57  ±  0.06  x  10“^, 
3.54  ±  0.004  x  10”^,  3.51  ±  0.02  x  10“^,  3,27  ±  0.06  x  10”^,  respectively.  A  measure¬ 
ment  of  the  A2!^  -  X?r[H2  band  with  N2  as  the  diluent  yielded  an  oscillator  strength 
of  3.72  ±  0.23  x  10~4.  This  value  was  chosen  for  comparison  with  values  in ’the  liter- 
ature  because  the  tt^ y 2  state  is  more  heavily  populated  at  room  temperature  than  the 
2tt3/2  state.  The  agreement  is  good  but  slightly  lower  that  those  considered  to  be 
the  most  reliable.  However,  it  should  be  noted  that  this  study  is  unique  since  the 
oscillator  strengths  were  measured  directly  from  individual  lines  in  high  resolution. 
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With  the  use  of  the  broadening  data  obtained  with  the  continuum  lamp,  the  oscillator 
strengths  in  the  model  were  adjusted  so  that  transmissions  predicted  by  the  model 
agreed  with  the  experimental  data  obtained  with  the  hollow  cathode  lamp.  The  values 
for  the  ^nx/2  an< d  "3/2  were  3.53  x  10“^  and  3.20  x  10~^  respectively.  These  values 
were  low  relative  to  the  continuum  values  because  of  the  presence  of  radiation  in 
the  hollow  cathode  lamp  due  to  species  other  than  NO.  If  an  estimate  is  made  of  the 
effect  of  this  excess  radiation,  the  oscillator  strengths  for  these  transitions  are 
3.76  x  10"^  (^1/2)  an<*  3.43  x  10“^  (^3/2).  The  source  temperature  used  in  this 
and  all  calculations  involving  the  hollow  cathode  lamp  was  600  K. 

Sufficient  measurements  were  conducted  and  compared  with  model  predictions  of 
transmission  that  the  optical  system,  based  on  a  hollow  cathode  lamp  and  associated 
spectral  model,  can  be  considered  calibrated  to  measure  NO.  If  accurate  temperature 
and  pressure  data  are  available  and  if  the  NO  absorption  is  large  relative  to  the 
noise,  then  measurements  with  accuracies  of  at  least  -  20%  are  possible. 

Similarly,  if  a  proper  model  is  used,  the  capillary  discharge  lamp  (see  Section 
IV. D.)  can  also  be  considered  calibrated. 

In  order  to  reach  this  conclusion  for  the  hollow  cathode  lamp,  data  were  ob¬ 
tained  and  processed  in  the  following  manner.  Two  concentrations  of  NO  in  N^  (93.3 
ppm  and  470  ppm,  independently  certified)  were  introduced  into  the  static  cell  at 
pressures  ranging  from  47  Torr  (6.3  kPa)  to  1000  Torr  (133  kPa) .  The  first  and 
second  bandhead  transmissions  of  the  Y(0,0)  transitions  were  measured  and  compared 
with  model  predictions.  The  averaged  errors  in  concentration  betw>  n  predicted  and 
measured  were  -0.3%  and  -2.23%  for  the  first  and  second  bandhead  data  sets,  respect¬ 
ively.  For  the  quartz-bed  heat  exchanger  and  the  flat  flame,  the  temperature  and 
concentration  profiles  were  divided  into  zones,  and  the  transmissions  through  each 
zone  were  calculated  using  the  average  temperature  and  concentration  in  each  zone. 

For  the  quartz-bed  heat  exchanger,  the  averaged  errors  in  concentration  for  the  total 
data  set  based  on  first  and  second  bandhead  transmissions  were  +7.14%  and  +6.79%, 
respectively.  Similarly,  for  the  flat  flame  data,  the  averaged  errors  for  the  data 
sets  were  +7.04%  and  -2.19%,  respectively.  The  maximum  error  encountered  for  any 
single  measurements  was  22.3%.  Finally,  for  the  high  resolution  continuum  measure¬ 
ments,  the  averaged  error  of  the  set  was  +3.6%  in  concentration.  It  must  be  noted 
that  these  errors  reflect  the  errors  associated  with  the  optical  measurements,  model 
predictions,  and  the  temperature  and  probe  concentration  measurements. 

Finally,  an  empirical  calibration  of  the  infrared  gas  correlation  spectrometer 
was  performed  (Appendix  A).  Because  this  instrument  was  originally  designed  for  stack 
monitoring,  i.e.,  low  temperatures  and  high  densities,  it  was  not  well-suited  for  the 
measurements  of  interest  here.  The  low  temperature  data  indicated  that  the  instru¬ 
ment  was  20%  more  sensitive  relative  to  the  calibration  previously  used  in  jet  com¬ 
bustor  measurements.  This  variation  is  attributed  to  changes  in  grating  alignment. 

A  dependence  of  the  calibration  on  broadening  gas  was  observed  and  determined.  For 
temperatures  up  to  900  K,  the  calibration,  within  the  scatter  of  the  instrument  output, 
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remained  constant.  Above  900  K,  a  significant  decrease  in  sensitivity  was  observed 
This  dependence  is  most  likely  due  to  significant  changes  in  the  populations  of  the 
lines  selected  by  the  grating  assembly.  However,  sufficient  data  were  obtained  to 
allow  measurements  to  be  made  at  high  temperatures  if  high  NO  seeding  of  the  media 
is  used. 
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INTRODUCTION  AND  SUMMARY 

As  a  result  of  tests  perforated  at  Wright  Patterson  Air  Force  Base^  the 
EPA  Smokestack  Instrument^^  was  requested  to  be  used  on  additional  jet  exhaust 
tests  under  controlled  conditions.  Tests  were  performed  using  a  flat  flame 
burner,  flowing  gas  heater  and  in  the  future,  will  be  made  on  a  jet  burner 
and  a  modified  combustor  can  from  a  Pratt  and  Whitney  FT 12  combustor.  This 
phase  of  the  tests  deals  with  the  calibration  procedures  used  for  the  flat 
flame  burner  and  flowing  gas  heater  conducted  at  United  Technology  Research 
Center  (UTRC) . 

The  EPA  Instrument  is  a  laboratory  instrument  designed  to  detect  SO2,  CO, 
HCL,  HF  or  NO  contaminants  being  emitted  from  a  smokestack.  The  details  of  the 
Instruments  are  well  documented^  ’  ^  and  are  not  repeated  here.  The  instrument 
uses  the  principle  of  gas  correlation  to  detect  small  amounts  of  contaminant 
gases  and  reduce  or  eliminate  the  effect  of  other  gases  such  as  H20,  C02  etc. 

For  this  series  of  tests  it  is  being  used  to  detect  NO  in  the  exhaust  aircraft 
jet  engines  for  temperatures  up  to  approximately  2000K. 

Three  different  types  of  tests  were  performed.  Static  cell  tests  were 
performed  using  a  sample  cell  to  contain  the  gas  at  room  temperature.  At  UTRC 
a  glass  cell  10  cm  long  with  NaCl  windows  was  filled  with  NO  premixed  with  either 
Ar  or  N2.  At  Ford  two  cells  were  used,  each  with  a  20  cm  path  length.  One  cell 
with  A1203  window,  the  other  a  glass  cell  with  A1203  windows.  The  samples  were 
NO  premixed  with  N2.  The  Flowing  Gas  Heater  tests  at  UTRC  consisted  of  flowing  a 
carrier  gas  of  either  N2  or  Ar  through  a  heater  unit  into  a  sample  area  that  is 
17  by  9.2  cm.  Known  amounts  of  NO  were  seeded  into  the  carrier  gas  at  tempera¬ 
tures  up  to  800R.  Buffer  gases  surrounded  and  contained  the  smnple  gases  and 
flooded  the  remainder  of  the  sample  area.  NaCl  windows  allowed  the  light  beam 
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to  traverse  the  sample  area  while  confining  the  gases.  The  Flat  Flame  Burner 
tests  at  UTRC  consisted  of  flowing  an  argon  carrier  gas  containing  H2  and 
O2.  The  H2  +  O2  was  burned  at  the  proper  ratio  and  flow  rates  to  create 
the  desired  temperatures  from  1000K  up  to  1600K.  The  argon  was  seeded  with 
known  amounts  of  NO.  The  same  buffering  system  and  external  housing  used 
for  the  flowing  gas  heater  were  used  for  the  flat  flame  burner. 

Because  of  the  versatility  designed  into  the  EPA  instrument,  it  is  not 
optimized  for  any  particular  gas  contaminant,  physical  pathlength,  particular 
environment  or  temperature.  It  was  designed  to  be  relatively  insensitive  to 
NO  sample  temperature  changes  up  to  900K.  This  was  achieved  by  selecting  a 
spectral  bandpass  that  contains  NO  absorption  lines  that  will  increase  in 
strength  and  some  that  will  decrease  as  the  sample  temperature  increases.  This 
results  in  the  sum  of  the  strengths  of  the  lines  remaining  relatively  constant 
for  the  temperature  range  from  300  to  900K.  Above  900K  the  response  of  the 
instrument  decreases  with  increasing  temperature.  The  present  tests  were 
performed  to  calibrate  the  instrument  and  to  determine  if  an  empirical  relation- 
ship  could  be  formulated  to  account  for  the  temperature  effect.  A  curve  has 
been  determined  to  relate  the  original  calibration  to  the  higher  temperatures. 

The  instrument  was  found  to  be  more  responsive  by  15  to  207.  to  samples  pressurized 
with  N2  under  static  conditions  at  room  temperature  than  when  used  at  Wright 
Patterson  Air  Force  Base.  This  is  due  primarily  to  a  shift  in  the  bandpasses 
of  the  grating  box.  The  instrument  was  also  found  to  be  more  responsive  to 
samples  diluted  with  Ar  than  to  samples  pressurized  with  N2. 

Differences  of  as  much  as  287.  occurs  for  room  temperature  samples  between 
static  cell  tests  and  flowing  gas  heater  tests.  These  differences  are  most 
] ikely  due  to  improper7  purge- flow. 
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INSTRUMENT  REFURBISHING 

The  instrument  was  received  from  EPA  with  no  apparent  damage  in  shipment. 

A  few  mirrors  that  were  broken,  previous  to  shipment,  were  replaced  with 
mirrors  of  similar  optical  quality.  The  instrument  was  originally  designed 
to  be  used  on  a  smokestack,  with  an  across-stack  path  using  a  retro* reflector 
to  doublepass  the  sample  area.  With  the  type  of  test  conducted  at  United 
Technology  Research  Center  we  could  not  doublepass  the  beam  because  of  the 
physical  dimensions  of  the  apparatus  involved.  Consequently,  a  source  section 
was  built  to  be  placed  on  the  opposite  side  of  the  sample  section  to  direct 
the  beam  through  the  sample  area  to  the  main  instrument.  This  is  similar  to 
the  technique  used  at  Wright  Patterson  Air  Force  Base  when  the  same  instrument 
was  used  for  NO  measurements  on  a  simulated  jet  aircraft  engine.  The  source 
optics  consisted  of  some  focusing  mirrors  and  an  electronic  reference  that 
was  actuated  by  the  435  Hz  chopper  that  also  chopped  the  light  beam  from  a 
Nernst  source.  The  main  instrument  was  checked  for  optical  alignment.  The 
transfer  optics  had  to  be  adjusted  when  the  grating  box  was  placed  into  the 
instrument.  The  mask  on  the  transfer  optics  was  removed  in  order  to  increase 
the  throughput  of  the  instrument. 

The  grating  box  is  the  same  one  U6ed  for  the  Air  Force  Wright  Patterson 
experiment.  The  grating  box  passes  three  spectral  intervals  where  NO  absorption 
lines  occur.  Figure  1  shows  three  spectral  curves  over  the  spectral  Interval 
of  interest.  The  upper  panel  shows  the  water  vapor  absorption  by  laboratory 
air  at  room  temperature.  The  middle  panel  shows  the  spectral  bandpasses  of  the 
grating  box,  and  the  third  panel  shows  the  absorption  by  NO  in  those  bandpasses. 
As  can  be  seen  by  the  middle  panel,  the  spectral  bandpasses  mask  out  the  strongly 
absorbing  R2O  lines.  In  the  spectral  interval  between  1925  and  1940  cm'*  the 
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TRANSMITTANCE 


Figure  1.  Spectral  transmittance  curves  over  the  spectral  interval  passed 

by  the  NO  grid  assembly.  The  upper  panel  shows  the  H2O  absorption 
by  laboratory  air  at  room  temperature.  The  middle  panel  shows  the 
spectral  bandpasses  of  the  grid.  The  lower  panel  shows  the  absorp- 
by  NO  in  the  bandpasses  of  the  grid.  The  dotted  lines  Indicate  the 
bandpass  shift. 


Che  NO  absorption  lines  increase  in  strength  as  the  temperature  Increases. 

In  the  spectral  interval  between  1895  and  1920  cm”^  the  NO  absorption  lines 
decrease  in  strength  as  the  temperature  increases.  The  total  absorption 
tends  to  remain  constant  as  the  temperature  changes  in  the  range  from  300  to 
900K.  The  gracing  box  was  placed  into  the  instrument  and  optically  checked. 

The  InSb  detector  was  installed,  however,  the  dewar  had  leaked  causing 
the  window  of  the  detector  to  frost  with  water  condensation.  Our  detector 
department  attempted  to  reevacuate  the  dewar.  The  pins  had  epoxy  deposited 
around  them  to  stop  the  leak  that  was  created  by  the  cracks  in  the  glass  near 
the  pins.  The  dewar  was  tubulated,  evacuated  and  pumped  for  about  four  days, 
then  the  tubulation  was  sealed.  The  dewar  subsequently  was  tested  and  found 
to  hold  liquid  nitrogen  for  four  hours.  However,  because  the  dewar  was  filled 
with  llguid  N2  to  the  top  the  epoxy  cracked  and  allowed  the  dewar  to 
leak  again. 

A  replacement  detector  that  was  not  optically  suited  for  the  Instrument 
because  of  its  element  size,  was  installed  in  the  instrument.  The  area  of  the 
element  of  this  detector  is  4  sso  x  0.4  mo  ,  whereas  the  area  of  the  original 
detector  element  was  1.3  mm  in  diameter.  The  original  preamps  of  the  Instrument 
were  found  to  be  defective  in  that  two  of  the  operational  amplifiers  had  to  be 
replaced.  A  new  preamp  was  used  with  the  replacement  detector.  The  rest  of  the 
electronics  were  found  to  be  functioning  very  well. 

Because  of  the  late  arrival  of  the  Instrument  and  the  length  of  time  that 
It  took  to  resolve  the  electronic  and  detector  problems,  limited  time  was 


available  Co  check  ouC  Che  instrument  and  make  Che  needed  repairs.  The  gracing 
box  was  inscalled  in  che  inscrumenc  and  Che  cransmiccance  of  Che  correlacion 


cell  was  checked.  This  was  done  by  comparing  che  carrier  signal  when  che 
correlacion  cell  was  noc  in  che  beam  Co  che  carrier  signal  when  che  correla- 
C ion  cell  was  in  che  beam.  We  found  chaC  Che  cransmiccance  of  Che  correlacion 
cell  was  577.,  which  was  che  same  as  when  used  on  che  Wrlghc  Paccerson  Air  Force 
Cescs.  This  cesc  seemed  Co  be  a  good  indicacor  ChaC  Che  correlacion  cell  and 
Che  gracing  box  were  in  good  specCral  alignment  and  ChaC  Che  correlacion  cell 
had  probably  noc  changed  lcs  specCral  characCerlscics. 

A  calibracion  cesc  was  Chen  performed  Co  check  che  previous  calibracion. 
The  calibracion  curve  used  ac  WrighC  Paccerson  is  shown  in  Figure  2.  A  sample 
of  1000  ppm  in  a  20  cm  sample  cell  produced  a  V'  of  68  x  10”^.  From  Che  pre¬ 
vious  calibracion  chis  V'  would  represenC  an  absorber  chickness  of  2.1  x  10“^ 
acm  cmS7p  .  Whereas  che  sample  in  che  sample  cell  was  1.8  x  10~2  acm  cm  . 
The  inscrumenc  was  indicacing  approximately  157.  higher  chan  che  sample  in  the 
sample  cell.  The  short  term  noise-equivalent  V*  was  4  x  10  ,  which  gave  a 

noise-equivalent  u  of  1.1  x  10"  ^  atm 


Figure  2.  Calibration  curve  for  NO  used  at  Wright  Patterson 
Air  Force  Base,  instrument  response  V'  vs  absorber 
thickness  u. 


TEST  AT  TTTRC 


The  instrument  arrived  ac  UTRC  in  Connecticut  in  very  good  shape,  with 
no  damage  during  shipment.  The  optics  were  set  up,  and  mirrors  had  to  be 
cut  to  accomodate  some  of  the  thermocouple  connectors  that  were  sticking  out 
of  the  chimney. 

Figure  3.  shows,  an  optical  diagram  of  the  equipment.  Sodium  chloride 
windows  were  used  on  the  flat  flame  burner,  and  an  image  of  the  Nernst  was 
formed  in  the  center  of  the  flat  flsme  burner.  Mirror  MS  then  focused  the 
image  on  to  mirror  XI  of  the  EFA  Instrument.  A  glass  cell  with  sodium 
chloride  windows  could  be  inserted  on  the  source  side  of  the  flat  flame 
burner  in  order  to  test  the  instrument  under  static  sample  conditions.  A 
burner  with  a  sintered  copper  top  produced  a  flat  flame  by  burning  H2  +  (>2 
in  an  Ar  carrier.  The  flat  flame  burner  was  replaced  with  a  flowing  gas 
heater  wherein  hot  gases  as  high  as  900K  could  be  introduced  into  the 
infrared  monitoring  beam. 

Tests  using  the  static  cell  were  done  by  first  evacuating  the  cell  in 
order  to  record  a  zero  and  then  filling  it  with  a  sample.  Two  types  of 
samples  were  used,  one  with  argon  gas  as  a  broadener  and  the  other  with 
nitrogen  gas  as  a  broadener.  After  the  sample  was  introduced  into  the  cell 
and  the  signal  recorded  on  the  recorder,  the  cell  was  evacuated  to  record 
another  zero.  The  two  zeros  on  the  chart  were  used  to  determine  the  zero 
position  for  the  instrument  when  the  gas  was  in  the  sample  cell.  If  any 
drift  had  occurred  it  could  be  compensated  for  by  drawing  a  line  through  the 
two  zeros.  Typical  shifts  were  less  than  1  *  10"4  for  V'  and  0. 3  x  10“  3  atm  cm 
for  u. 
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the  Plowing  Gas  Heater  was  used. 


A  similar  procedure  was  used  for  the  flowing  gas  heater  and  the  flat 
flame  burner.  The  flowing  gas  heater  was  heated  up  to  the  temperature 
desired  with  either  Ar  or  N2  as  a  carrier  gas;  the  temperature  of  the  center 
of  the  flowing  gas  heater  sample  section  was  measured  and  recorded.  A  zero 
level  was  recorded  on  the  chart,  then  the  sample  was  seeded  with  NO  to  a 
predetermined  amount  and  the  NO  concentration  at  the  center  of  the  sample 
area  was  determined  by  T7TKC  using  quartz  probe  techniques.  The  signal  was 
recorded  on  Che  recorder  chart  paper  and  the  NO  gas  was  turned  off.  The 
Instrument  then  gave  another  zero  level.  The  two  zeros  on  either  side  of  the 
signal  recorded  for  the  sample  were  used  in  order  to  determine  the  zero  when 
the  NO  sample  was  in  the  sample  section.  A  similar  procedure  was  used  for 
the  flat  flame  burner.  The  burner  was  brought  up  to  temperature,  heating  the 
argon  carrier  gas.  Two  zero  levels  were  recorded,  one  before  and  one  after 
the  carrier  gas  was  seeded  with  NO.  In  the  flowing  gas  heater  and  the  static 
cell  tests  both  nitrogen  and  argon  were  used  as  a  carrier  gas.  The  nitrogen 
was  not  used  for  the  flat  flame  burner  for  fear  of  creating  other  NO  or  N0x 
products  which  would  contaminate  the  sample  area. 

RESULTS 

The  techniques  used  to  reduce  the  data  are  very  similar  to  the  one  used 
for  the  tests  conducted  at  Wright  Patterson  Air  Force  Base.  The  details  of 
the  data  reduction  can  be  found  in  their  report  (2).  Only  a  brief  description 
will  be  given  here. 

The  signal,  V' ,  from  the  EPA  instrument  is  proportional  to  the  NO  absorber 
thickness  for  small  samples  at  constant  temperature  and  concentrations.  For 
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Che  samples  encountered  In  these  tests,  u  *  2  x  10" ^  atm  cm  STP,  the  response 
of  the  Instrument  Is  linear  as  can  be  seen  from  the  calibration  curve  in  Figure 
2.  When  the  sample  has  temperature  and  concentration  gradients,  the  sample 
path  can  be  considered  to  be  divided  into  small  increments  such  that  each 
incremental  path  is  at  a  constant  temperature  and  concentration.  The  absorber 
thicknesses  for  each  increment  can  be  determined  and  because  of  the  linear 
response  of  the  instrument  their  sum  will  give  the  total  absorber  thickness 
in  the  sample  path.  In  order  Co  compare  the  results  of  the  instrument  with 
the  actual  sample,  the  temperature  and  concentration  gradients  must  be  known. 

A  computer  program  was  created  to  calculate  the  absorber  thicknessess  of  the 
samples  at  the  elevated  temperatures. 

To  determine  the  absorber  thickness,  values  of  temperature  and  concentra¬ 
tion  were  determined  from  gradient  and  mole  fraction  gradient  curves  supplied 
by  TJTRC.  For  the  flowing  gas  heater  tests,  four  temperature  profiles  were 
supplied  and  shown  in  Figure  4.  Each  curve  is  flat  for  a  distance  out  to 
approximately  4  or  5  cm  from  the  centerline  of  the  burner.  The  actual 
centerline  temperatures  of  the  samples  are  listed  above  the  curves  for  which 
they  were  used.  These  profiles  were  determined  within  an  hour  from  the  tests 
using  Che  EPA  instrument  were  performed.  The  profiles  are  a  result  of  probe 
measurements  made  by  the  T7TRC  staff  on  samples  whose  ten?eratures  and  con¬ 
centrations  were  similar  to  those  used  when  making  the  measurements  with  the 
EPA  instrument.  From  the  curves  in  Figure  4,  the  actual  centerline  temperatures 
of  the  NO  sample  can  be  off  as  much  as  107,  from  the  centerline  temperature 
of  the  curve  that  was  used  to  determine  the  sample  absorber  thickness.  Figure  5 

shows  the  mole  fraction  profile  for  three  different  temperatures.  The  300K 
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Figure  4.  Temperature  Profile*  for  the  Flowing  Gas  Beater. 

The  sample  numbers  for  which  the  curves  were  used 
are  written  below  each  curve.  The  temperature  at 
the  center  of  the  sample  is  listed  above  each  curve. 
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curve  was  used  for  the  room  temperature  samples  and  the  dotted  curve  was  used 
for  the  elevated  temperatures. 


The  absorber  thickness,  u,  for  a  sample  of  constant  temperature  and  con¬ 
centration  can  be  determined  from  the  UTRC  data  from 


where 


u(atm  cmgYp)  ■  p(atm)  L(cm)  ^ 

u  -  Absorber  thickness  In  atm  cm  g^-p 
p  •  Absorber  pressure  in  atm  *  c(ppm)  x  10"^ 
c  ■  Absorber  concentration  in  ppm 
L  ■  Path  length  in  cm 

0  •  Temperature  in  K  (2) 


However,  if  the  concentration  and  temperature  varys  over  the  sample  path,  then 
the  absorber  thickness  for  a  small  incremental  path  length  is  given  by 

Ui  -  pj^Li  273 
8i 

then  the  total  absorber  thickness  for  the  path  can  be  found  by  summing  over  the 
entire  path. 

u  -  z  ■  2Pi  £Li  ^12  (4a) 

1  1 

and  if  the  path  la  divided  into  equal  increments,  aL, 

u  -  aL  273  x  10“6  5  ii 

1  ei 
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The  concentration  can  be  found  from  the  mole  fraction  normalized  concentration 
profile,  M,  and  the  centerline  concentration,  C,  from  the  following  expression. 

ci  "  Mic  (5) 

Therefore 

uc  -  £LC  273  x  10*6  ? 

ei 

where  uc  denotes  the  calculated  absorber  thickness.  The  values  used  in  the 
calculations  were  determined  at  0.5  cm  increments  from  Figures  4  and  5.  In¬ 
tegration  was  done  from  the  center  of  the  sample  through  10  cm.  The  results 
were  doubled  on  the  assumption  that  the  profile  is  symmetrical  about  the 
center  line.  Some  of  the  representative  curves  of  previous  data  indicate  that 
this  might  not  be  an  accurate  assumption.  Asymmetry  might  cause  errors  as  large 
as  871. 

Because  the  calibration  curve  for  the  EPA  instrument  has  been  determined 
for  temperatures  in  the  range  of  300K  to  900K  it  is  not  expected  to  give 
absolute  results  for  samples  pressurized  with  Ar  or  for  samples  vrtiose  tempera¬ 
tures  are  higher  than  90QK.  In  addition  it  has  previously  been  shown  that 
the  EPA  instrument  is  giving  values  that  are  157.  higher  than  the  original 
calibration  curve.  The  signal,  V* ,  from  the  EPA  Instrument  will  be  used  to 
determine  a  measured  equivalent  absorber  thickness,  urn,  from  the  original 
calibration  curve.  These  measured  equivalent  absorber  thicknesses  will  be 
compared  to  the  calculated  absorber  thicknesses  to  determine  an  empirical 
temperature  relationship  between  the  Instrument  response  and  the  calculated 
values  for  samples  diluted  with  Ar. 


Table  1  Uses  the  results  of  the  static  cell  tests.  The  first  column 
lists  the  broadening  gas.  The  second  columns  lists  the  results  of  the  EPA 
instrument.  The  third  column  lists  the  calculated  absorber  thickness  and 
the  fourth  column  lists  the  ratios  of  um/uc.  The  noise  equivalent  u  is 
approximately  2  x  10  The  values  obtained  for  the  fifth  and  sixth  samples 
in  Table  1  are  for  small  samples  diluted  with  Nj.  The  noise  equivalent 
u  represents  a  larger  percentage  error  than  for  the  samples  for  higher 
concentrations.  The  first  three  samples  listed  on  Table  1,  consequently  would 
give  the  more  accurate  result.  This  indicates  the  EPA  instrument  is  giving 
values  of  about  207.  higher  which  agrees  with  the  results  that  we  obtained  at 
Ford  prior  to  the  UTRC  tests.  The  seventh  and  eighth  samples  are  for  the 
calibration  test  performed  at  Ford  prior  to  and  after  the  UTRC  tests,  res¬ 
pectively.  The  last  two  samples  are  for  argon  as  a  broadener  gas  and  the 
values  of  um/uc  are  larger  by  about  357.;  It  has  been  shown  (3)  that  the 
broadening  abilities  of  Ar  and  N?  on  CO  are  different  and  a  similar  effect 
on  NO  might  explain  the  observed  effect  above. 

The  results  from  the  flowing  gas  heater  tests  are  shown  in  Table  2.  The 
upper  portion  of  the  table  shows  the  nitrogen  carrier  data,  the  lower  portion 
the  argon  carrier  data.  The  left  hand  column  list  the  teat  number;  some  of 
the  tests  were  repeated  at  a  later  date  and  consequently  given  the  same  test 
number.  The  third  column  lists  the  temperature  at  the  centerline  of  the 
sample.  The  next  column  lists  the  centerline  temperature  of  the  curve  chat 
was  used  in  the  calculation.  The  next  column  lists  the  centerline  concentra¬ 
tions  supplied  by  UTRC.  The  next  column  lists  the  measured  equivalent  absorber 
thickness,  um,  in  atm  emg^p,  as  determined  with  the  EPA  instrument.  Next  column 
gives  the  calculated  absorber  thickness  based  on  the  probe  curves  given  by  UTRC. 
The  last  column  gives  the  ratio  of  um/uc.  The  lower  portion  of  the  table  lists 
the  results  when  argon  is  used  as  a  carrier  gas. 


TABLE  1 


Static  Cell  Results 


Broadener 

Gas 

un 

(atm  co  STP 
xlO3) 

uc 

(atm  co  stp 
xlO3) 

um/uc 

N2 

24.0 

19.2 

1.25 

N2 

23.5 

19.2 

1.22 

»2 

23.6 

19.2 

1.22 

N2 

4.2 

3.64 

1.15 

n2 

8.5 

7.56 

1.12 

N2 

21.2 

18.5 

1.15 

n2 

21.2 

18.5 

1.15 

Ar 

24.5 

18.3 

1.34 

Ar 

24.5 

18.3 

1.34 
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TABLE  2 


Flowing  Gas  HeaCer  Results 


Test 

Carrier 

e(K> 

9(K) 

C 

um 

uc 

um/uc 

# 

Gas 

atm  cm  stp 

atm  cm  ctp 

Measured 

Curve 

Ppm 

xl03 

xl03 

20T 

n2 

296 

669 

12.0 

9.22 

1.30 

21  T 

n2 

296 

1315 

24.5 

18.1 

1.35 

22  T 

n2 

608 

611 

1025 

7.5 

6.69 

1.12 

23  T 

n2 

573 

611 

1850 

16.6 

13.3 

1.25 

24  T 

n2 

723 

746 

1580 

8.9 

8.53 

1.04 

25  T 

n2 

703 

746 

2220 

14.6 

11.9 

1.23 

20  T 

n2 

296 

650 

11.8 

8.97 

1.32 

20  T 

n2 

296 

668 

12.1 

8.97 

1.35 

22A 

Ar 

678 

667 

1018 

7.4 

6.11 

1.21 

23A 

Ar 

603 

667 

2047 

18.0 

12.3 

1.47 

24A 

Ar 

793 

796 

1585 

8.8 

8.06 

1.09 

25A 

Ar 

800 

796 

2200 

13.9 

11.2 

1.24 

2QA 

Ar 

295 

15.6 

9.07 

1.72* 

20A 

Ar 

295 

15.2 

9.07 

1.68* 

*  Window  purge  was  most  likely  improperly  set. 
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The  last  column  Indicates  a  tendency  for  the  ratio  to  decrease  as  the 
temperature  Increases.  The  ratio  for  the  room  temperature  samples  are 
much  larger  than  the  ratio  for  the  samples  at  the  elevated  temperatures.  The 
calculated  results  for  sample  25T  may  be  reduced  If  the  temperature  curve, 

746K  were  corrected  to  the  temperature  of  the  centerline  of  the  sample,  703K. 

It  would  cause  the  ratio  to  decrease  coming  closer  to  the  results  of  sample 
24T.  The  values  obtained  for  the  nitrogen  samples  tend  to  be  lower  than  the 
samples  with  argon  as  the  carrier  gas  in  agreement  with  the  tests  on  the 
static  cell.  The  ratio  for  the  argon  sample  tests  also  tend  to  get  smaller 
as  the  temperature  is  increased. 

t 

For  the  flat  flame  burner  tests,  similar  type  calculations  were  carried 
out.  However,  the  temperatures  given  at  the  centerline  of  the  sample  are  not 
the  true  temperatures  of  gas  but  the  temperature  of  the  thermocouple.  These 
temperatures  have  to  be  corrected  because  the  thermocouples  are  radiating 
thereby  decreasing  their  temperature  giving  a  result  that  is  lower  then  the 
actual  gas  temperature.  The  temperatures  that  were  used  were  determined  from 
curves  given  to  us  by  UTRC.  They  are  shown  in  Figure  6.  Four  different 
temperature  profiles  are  represented.  Listed  above  each  curve  is  the  corrected 
centerline  temperature. 

These  curves  were  determined  for  similar  sample  conditions  in  the  flat 
flame  burner.  However  they  were  obtained  prior  to  May  1978,  six  months  before 
the  tests  with  the  EPA  instrument.  Tt  was  assumed  that  they  would  give  a  true 
representation  pf  the  temperature  profile  of  the  sample  conditions  of  the  present 
tests.  The  Centerline  temperatures  for  the  present  samples- were  slightly 


different.  However,  corrections  to  these  temperature  curves  were  made  to  get 
representative  temperature  profiles.  We  have  interpolated  between  the  curves 
for  the  different  samples  using  a  linear  interpolation  from  296®K  to  the 
calibration  curve.  The  radiation  correction  was  done  in  the  following  matter: 

The  correction  was  assumed  to  follow  the  centerline  temperature  of  the  sample 
raised  to  some  power,  probably  near  the  third  power.  The  temperature 
corrections  for  the  four  different  curves  shown  were  plotted  against  the  re¬ 
corded  centerline  temperature  on  log-log  paper.  The  slope  of  the  curve  gave 
the  exponent  for  the  following  equation: 

6  +M i(6)n.  (6) 

From  the  slope  of  the  curve,  n  was  found  to  be  2.94  and  by  using  the  values 
that  were  used  to  determine  the  curve,  B  was  found  to  be  1.24  x  10"^.  With 
the  use  of  this  equation  the  intermediate  temperatures  were  corrected.  Then 
with  the  use  of  the  computer  we  determined  the  values  of  the  integrated  u,  for 
the  different  samples.  Figure  7  shows  the  curve  determined  using  probe  techniques 
prior  to  May  1978  that  was  used  for  the  mole  fraction  profile  for  the  flat  flame 
burner  calculations. 

Table  3  lists  the  results  for  the  flat  flame  burner  tests.  The  second 
column  in  Table  3  shows  the  uncorrected  centerline  temperature;  the  third 
column  lists  the  corrected  center  line  temperature.  The  fourth  column  lists 
the  value  of  the  measured  equivalent  absorber  thickness  as  measured  by  the 
EPA  instrument,  the  sixth  list  the  calculated  values  and  the  seventh  column 
gives  the  ratio  of  um/uc.  As  can  be  seen  from  Table  3  the  data  are  grouped 
into  three  different  temperatures  regions.  As  the  temperature  increases, 
the  tendency  for  the  ratio  is  to  decrease,  similar  to  the  flowing  gas  heater 
results. 


A-21 


Distance  from  Centerline  (cm) 

Figure  7,  Mole  fraction  profile  for  the  Fiat  Flame  Burner 


TABLE  3 

Flat  Flame  Burner  Results 


Test 

# 

800 

Oncorrected 

000 

Corrected 

C 

DPtB 

un 

a£jq  cttSTP 
xlO3 

uc 

atm  cmsTP 
xlO3 

um/uc 

2  OF 

940 

1008 

1389 

5.8 

6.41 

0.905 

2  IF 

935 

1002 

1954 

7.8 

9.08 

0.859 

20FA 

934 

1001 

2392 

7.2 

11.1 

0.647 

21FA 

939 

1007 

3775 

12.0 

17.5 

0.688 

20FB 

919 

983 

3248 

13.0 

15.4 

0.846 

21FB 

919 

983 

2293 

8.2 

10.9 

0.756 

21FC 

919 

983 

3606 

16.0 

17.1 

0.937 

22FA 

1189 

1325 

3089 

6.0 

10.7 

0.562 

23FA 

1180 

1313 

4246 

8.5 

14.8 

0.  574 

26FA 

1227 

1377 

4018 

5.8 

13.4 

0.434 

27FA 

1227 

1377 

6074 

8.5 

20.2 

0.420 

27FA 

1189 

1325 

6074 

7.2 

20.2 

0.356 

22FA 

1189 

1325 

3089 

8.4 

10.7 

0.787 

23FA 

1189 

1325 

4246 

10.8 

14.7 

0.736 

24FA 

1393 

1600 

4136 

5.1 

12.0 

0.425 

25FA 

1396 

1615 

6125 

8.5 

17.7 

0.480 

24FA 

1381 

1593 

4136 

7.5 

12.1 

0.618 

25FA 

1381 

1593 

6125 

10.6 

20.2 

0.524 

25FA 

1381 

1593 

6125 

11.8 

20.2 

0.584 

25FB 

1381 

1593 

6125 

11.5 

20.2 

0.569 
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TESTS  AT  FORD 


After  the  tests  &t  UTRC,  additional  tests  were  performed  at  Ford 
that  were  not  done  prior  because  of  time  limitations.  Also  tests  were 
performed  to  help  explain  some  of  the  effects  that  were  observed. 

The  grating  box  was  tested  by  directing  a  beam  of  light  through  it  and 
into  a  grating  spectrometer.  The  spectral  transmittance  curve  of  the  grating 
box  was  displayed  on  a  chart  recorder.  The  spectral  bandpasses  of  the  grating 
box  were  compared  to  the  previous  spectral  bandpasses  and  found  to  have  been 
displaced  approximately  2  em~^  toward  lower  wavenumber.  The  displacement  of 
the  spectral  bandpasses  probably  resulted  in  the  calibration  shift.  It 
c  ould  have  changed  the  response  of  the  lnstrtsnent  to  H2O;  however,  during 
our  tests  at  DTRC,  when  the  H2  +  Oj  burner  was  fired  creating  H20,  no  change 
in  signal  was  observed,  indicating  chat  the  instrument  still  had  good  H20 
rejection.  Additional  calibration  tests  were  performed  at  the  Ford  facility 
confirming  the  previous  shift  of  15H  in  the  calibration  curve  from  the  earlier 
data.  The  spectral  transmittance  curve  of  the  correlation  cell  was  obtained 
and  compared  to  its  previous  transmittance  curve,  and  no  observable  difference 
was  noted. 

There  is  a  noticeable  effect  when  Ar  is  used  as  a  carrier  or  broadening 

f  31 

gas  as  compared  to  .  We  have  observed  a  similar  effect  between  Ar  and 
N2  with  CO  as  an  absorbing  gas.  N2  broadens  the  CO  lines  more  than  Ar.  The 
wings  of  the  broadened  lines  would  have  more  absorption  that  the  wings  of 
less  broadened  lines.  Because  cylinders  of  Ar  were  not  available  to  check 
the  effect  of  the  broadening  ability  of  Ar  on  NO,  a  test  was  performed  to 
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determine  Che  effect  of  broadening  the  NO  lines  with  Nj.  A  sample  of  NO 
at  two  different  pressures  near  1  atm  was  tested.  A  1000  ppm  NO  in  N2 
sample  at  0.7  atm  in  a  20  cm  cell  was  pressurized  with  N2  to  1.0  atm  and 
the  signal,  V',  decreased  from  55  x  10"^  to  51.5  x  10*^.  A  decrease  would 
be  expected  because  the  wings  of  the  absorption  lines  contribute  to  V*  as  a 
negative  correlation,  when  the  lines  are  broadened  the  wing  absorption  is 
increased,  thereby  reducing  the  signal.  Because  Ar  is  not  expected  to 
broaden  as  much  as  N2,  the  wings  of  the  lines  will  absorb  less  for  Ar-broadenlng 
thereby  increasing  the  signal. 


CONCLUSIONS 


There  Is  a  calibration  shift  of  approximately  15  to  20%  higher  from  the 
calibration  used  for  the  Wright  Patterson  Air  Force  tests  in  1975.  This  is 
probably  due  to  the  spectral  shift  of  the  gratin  g  box  bandpasses  of  2  cm"* 
toward  lower  wavenumber.  The  N2 -broadened  NO  lines  are  broadened  more  than  the 
Ar-broadened  NO  lines  and  consequently  absorb  more  in  the  wings  of  the  lines. 
However  more  absorption  in  the  wings  of  the  lines  will  cause  a  decrease  in 
the  signal,  therefore  the  samples  diluted  by  Ar  give  larger  signals  than 
samples  pressurized  by  N2. 

The  ratio  um/uc  for  the  flowing  gas  heater  tests  are  considerably  higher 
than  for  the  static  cell  tests  when  Ar  is  used  as  the  broadening  gas.  This  is 
very  clear  when  the  ratio  um/uc  is  compared  in  Table  1  and  Table  2  for  room 
temperature  samples.  The  ratio  is  1.34  for  the  static  cell  tests  and  1.7 
for  the  flowing  gas  heater.  Personnel  at  UTRC  indicate  that  the  window  purge 
was  improperly  set.  The  difference  is  not  as  great  for  the  samples  pressurized 
with  N2.  The  N2  results  are  within  experimental  error. 

Figure  8  shows  a  plot  of  the  results  for  samples  that  are  pressurized  with 

Ar.  The  ratio  um/uc  is  plotted  against  temperature.  The  large  difference  in  the 

static  cell  tests  and  the  flowing  gas  heater  results  is  obvious.  There  is  a 

large  amount  of  scatter  in  the  data.  The  noise-equivalent  u  was  approximately 

.3 

1  x  10  atm  ctngTp,  5  to  10%  of  the  observed  signal,  and  accounts  for  some 
of  the  scatter.  The  use  of  standardized  temperature  and  concentration  profiles 
rather  than  using  a  profile  for  each  sample  could  add  5  to  10%  to  the  scatter. 

As  can  be  seen  from  Tables  2  and  3  the  centerline  temperatures  varied  by  as  much 


Figure  8.  Plot  of  um/uc  vs  Temperature  for  Ar»broadened  samples. 
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as  10%  from  t  he  centerline  temperatures  of  the  profile  curves  used  In  the 
calculation.  Another  source  of  error  that  could  add  up  to  7  to  10%  is 
the  assumption  that  the  profiles  are  symmetrical.  Previous  data  supplied 
by  UTRC  Indicate  that  the  profiles  may  not  be  syrmnetrical. 

A  smooth  curve  was  drawn  through  the  points  in  Figure  8.  The  curve 
shows  the  decrease  in  um/uc  as  the  temperature  increases.  This  apparent 
change  in  calibration  was  not  observed  for  temperatures  below  900K  in  the 
previous  tests  done  in  1975.  As  can  be  seen  in  Figure  2  for  temperatures 
up  to  900K,  under  static  conditions,  a  single  calibration  curve  could  be 
used.  The  absorption  lines  used  in  the  detection  of  NO  tend  to  account  for 
temperature  changes  because  some  of  the  absorption  lines  increase  in  strength  and 
and  same  decrease  in  strength  with  Increasing  temperature.  The  total  strength 
of  the  sum  of  the  lines  would  then  tend  to  remain  constant  over  the  temperature 
range  of  interest.  However  as  the  temperature  increases  above  1000K  the 
strengths  of  the  lines  all  start  to  decrease  and  the  response  of  the  instru¬ 
ment  falls  off.  For  the  subsequent  testson  this  contract  the  calibration 
curve  should  be  modified  by  a  factor  determined  from  the  curve  from  Figure  8. 

A  similar  plot  for  ^-broadened  samples  is  shown  in  Figure  9.  The 
differences  between  the  static  cell  tests  and  the  flat  flame  burner  are 
not  as  great  as  for  the  Ar-broadened  samples.  Zn  addition  the  difference 
between  the  static  cell  tests  and  the  flowing  gas  heater  tests  not  as 
apparent  for  Nj-broadened  samples. 


Plot  of  um/uc  v«  Temperature  for  ^“broadened  »*®pl*8 


Figure  9. 
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APPENDIX  B 


COMMENTS  ON  THE  PROBLEMS  IN  THE  PREVIOUSLY  REPORTED  SPECTRAL  MODEL 

The  review  of  the  spectroscopic  model  developed  by  ARO/AEDC  revealed  the 

following  major  problems: 

1.  There  is  an  error  in  equation  18  of  Ref.  B3  (the  same  error  is  reflected  in 
Refs.  B3-B8)  relating  fj'  "  to  fv'  ",  which  is  off  by  a  factor  of  4. 

2.  There  is  an  error  in  equation  17  of  Ref.  B3  (and  the  other  reports)  giving 
the  ratio  Nj"/N0,  which  is  off  by  a  factor  of  about  2. 

o  2 

3.  In  Ref.  B1  there  was  no  distinction  in  population  between  the  and  *3/2 
states  (i.e.,  Hund's  case  (b)  was  assumed)  although  this  was  partially  cor¬ 
rected  in  Ref.  B3  and  the  computer  program  delivered  to  UTRC,  but  without  the 
correct  normalization  in  the  denominator. 

4.  An  error  exists  in  the  equation  for  the  Honl -London  factor  for  Q33  lines. 

This  error  can  be  found  in  the  text  of  a  1973  report  (Bl) ,  the  text  of  a 
1976  report  (Ref.  B3) ,  and  the  computer  program  delivered  to  UTRC  in  1978. 

ARO  personnel  have  stated  that  this  error  was  not  always  present,  but  they 
are  uncertain  as  to  when  it  first  appeared.  The  error  results  in  a  mis¬ 
calculation  of  NO  concentration  by  about  15%  when  considered  alone,  and  an 
incorrect  temperature  dependence  for  the  absorption. 

5.  An  error  exists  in  Eq.  24  of  Ref.  B3  which  relates  the  measured  broadening 
parameter  to  the  collision  cross  section.  This  error  does  not  affect  the 
predicted  NO  values  from  the  model,  but  does  result  in  a  factor  of  10  error 
in  the  collision  cross  section  reported  in  Ref.  B7. 

In  addition  to  the  major  problems  listed  above,  a  number  of  minor  problems 

were  discovered  in  the  model  delivered  to  UTRC  which  do  not  have  much  numerical 

significance,  but  which  should  be  corrected. 

6.  The  Qj2  (30.5)  and  R^  (30.5)  lines  were  incorrectly  labeled. 

7.  The  lines  Q32  (3.5),  Rj^  (36.5),  (31.5)  and  P-^  (34.5)  through  P32 

(40.5)  should  be  added. 

8.  There  is  an  error  in  the  equation  for  Bv  in  subroutine  HONNUM  (HOnl-London 
factor  computation). 


Finally  UTRC  has  shown  that  the  use  of  theoretically  generated  line  locations 
offers  an  improvement  in  the  model  when  compared  with  the  experimental  line 
locations  as  given  by  Deezsi  (Ref.  B9) .  This  does  not  represent  an  error  in  the 
earlier  model,  but  does  represent  a  considerable  improvement,  particularly  at 
low  pressures  or  elevated  temperatures. 
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APPENDIX  C 


COMMENTS  ON  THE  EXPERIMENTAL  TECHNIQUE  OF  WISE  AND  FRECH 

Wise  and  Freeh  observed  the  amount  of  NO  decomposition  by  measuring  the 
formation  of  the  assumed  products,  N2  and  02-  Molecular  nitrogen  was  measured 
by  trapping  the  remaining  NO  in  liquid  nitrogen  and  using  'manometric'  tech¬ 
niques.  Oxygen  was  allowed  to  react  according  to 

NO  +  NO  +  02  -*•  2  N02  (lc) 

The  nitrogen  dioxide  was  measured  via  optical  absorption  techniques.  Several 
problems  with  their  paper  and  experimental  procedure  may  be  found.  First  of 
all,  no  comparison  is  made  between  the  N2  and  NO2  data.  Secondly,  no  allowance 
is  made  for  alternative  reactions  or  their  products.  For  example,  consider  the 
overall  reactions 

4N0  -+  2  N2  +  2  N02  (2c) 

or 

2N0  -*•  N20  4  1/2  02  (3c) 

for  which  experimental  data  must  be  interpreted  differently.  Finally,  and  quite 
importantly,  their  pressure  measurement  of  the  product  N2  appears  to  be  at 
best  extremely  difficult.  Since  very  low  levels  of  conversion  were  examined 
(usually  about  0.5  percent),  very  small  pressures  of  N2  are  expected,  typically 
about  one  torr.  In  order  for  their  technique  to  be  useful,  the  nitrogen  pres¬ 
sure  must  be  measured  accurately  in  a  background  of  approximately  one  torr  nitric 
oxide  (vapor  pressure  at  liquid  nitrogen  temperatures),  one  torr  of  oxygen  (if 
a  product),  and  about  1/4  torr  impurity  (likely  to  be  nitrogen).  Certainly,  a 
measurement  of  one  torr  nitrogen  with  these  other  gases  present  would  not  be  an 
easy  task.  If  the  dominant  impurity  in  the  initial  NO  is  N02  rather  than  N2  or 
N20  (stated  purity  is  99.93  percent  NO),  then  a  measurement  of  0.5  percent  de¬ 
composition  would  be  in  error  by  (.07  percent)  x  2/. 5  percent  =  28  percent  when 
the  N02  data  is  interpreted. 
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APPENDIX  D 


UTRC  SPECTRAL  COMPUTER  PROGRAM  DESCRIPTION  AND  LISTING 


D.I.  Program  Description  of  NO  Absorption  of  Continuum  Radiation 
D.I.a.  Introduction 


Several  different  versions  of  the  NO  spectral  model  have  been  developed. 

All  of  these  are  based  on  the  model  developed  by  Davis  et  al  (1976)a  and  sup¬ 
plied  to  us,  with  corrections  and  additions  made  by  us  as  noted  in  Appendix  B. 

The  plot  routines  are  substantially  different.  The  corrections  made  to  the 
original  program  developed  by  Davis  et  al  (1976a)  make  a  significant  differ¬ 
ence  in  the  amount  of  NO  necessary  to  produce  a  given  amount  of  absorption  - 
about  a  factor  of  2  at  room  temperature  and  atmospheric  pressure. 

The  different  versions  of  the  NO  spectral  model  may  be  summarized  as 
follows.  The  absorption  of  continuum  radiation  by  a  homogeneous,  isothermal 
layer  of  NO  diluted  in  a  foreign  gas  is  modeled  by  the  program  (DUSEK*)NO-ABS 
(described  in  D.I.),  in  conjunction  with  the  plot  routine  (DUSEK*)PLOT. 

PLOT 5  (detailed  in  D.IV.)  and  the  data  file  (DUSEK*)DATA-2 .  or  (DUSEK*)DATAENG . 
The  similar  problem  for  absorption  of  continuum  radiation  by  multiple  zones, 
each  of  which  is  homogeneous  and  isothermal,  is  treated  with  two  programs  not 
described  here,  (PAGE*)NO-ABS .  and  (PAGE*)PLOT . ,  along  with  the  data  file  (PAGE*) 
DATA-2 .  The  multiple  zone  includes  the  proper  accounting  of  the  complete  line 
profile  through  the  different  absorbing  layers  before  convolving  the  resultant 
with  the  slit  function.  The  model  describing  the  absorption  of  narrow-line 
radiation  is  called  (DUSEK*)  NO-SPECT.  (described  in  D.II.),  with  the  same  plot 
routine  and  data  files  as  for  (DUSEK*) NO-ABS . 

The  program  NO-ABS.  predicts  the  transmission  as  a  function  of  wavelength 
of  incident  continuum  radiation  by  a  known  optical  depth  of  nitric  oxide  at  a 
given  temperature  and  pressure.  This  program  was  written  specifically  for  pre¬ 
dicting  transmission  in  high  resolution  over  a  narrow  wavelength  region.  An 
additional  program,  PLOT.,  convolves  this  predicted  transmission  with  a  spectro¬ 
meter  slit  function  and  displays  the  results  graphically  for  comparison  with 
experimental  spectra. 

The  program  NO-ABS.  is  composed  of  a  main  program,  four  subroutines,  and 
a  block  data  file.  The  program  elements  for  the  main  body  and  the  block  data 
are  NO-ABS. MAIN  and  NO-ABS. BLOCK-DATA,  respectively.  The  element  names  of 
the  subroutines  are  FUPPER,  FLOWER,  HONNUN,  AND  WFUNC.  Each  of  these  names 
is  preceded  by  the  file  name  NO-ABS,  as  before.  A  definition  of  variables 
for  each  element  of  NO-ABS.  and  a  listing  of  the  elements  are  included  at  the 
end  of  the  program  description  section. 
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D.I.b.  Main 


The  main  program  is  divided  into  four  parts.  A  description  of  the  function 
of  each  part  follows. 

D.I.b.l.  Main,  Part  1^ 


Part  One  reads  parameters  such  as  the  optical  depth,  temperature,  pressure, 
broadening  information,  options  for  line  location,  region  of  interest  and  a 
code  that  describes  the  version  of  the  program  presently  in  use.  The  options 
available  for  line  location  include  either  a  theoretical  calculation  of  line 
center  or  the  use  of  experimentally  determined  line  center  wavenumbers.  The 
experimental  values  for  the  line  center  wavenumbers  were  obtained  from  Engleman 
et  al  (1970).  Theoretical  line  locations  using  spectroscopic  equations  and 
constants  for  NO  (Engleman)  yield  predicted  locations  that  differ  from  experi¬ 
mental  observation  by  less  than  the  room  temperature  Doppler  width  of  NO 
(.005  X  or  .0005  ran)  for  almost  all  lines.  The  wavenumber  region  of  interest 
is  defined  by  the  first  and  last  line  numbers  of  that  region.  The  present 
version  of  the  program  reads  474  lines.  Normally,  1  to  30  of  these  lines 
are  examined  in  a  given  run.  Also  contained  in  Part  One  of  the  main  program  is 
the  assignment  routine  for  the  f-number  of  the  transition.  Since  experimental 
evidence  indicates  that  transitions  originating  in  the  2*jy2  state  have  a  larger 
f-number  than  those  originating  in  the  ^"3/2  state,  provision  has  been  made  for 
a  spin  dependent  f  number. 

i: •  __M a in,_P art  2 

Part  Two  of  the  main  program  reads  and  decodes  the  input  line  designation. 
The  line  designation  consists  of  the  branch  of  the  transition  (P,  Q,  R  for 
J'-J"  =  -1,  0,  +1),  upper  spin  state,  lower  spin  state,  and  the  rotational 
level  in  the  lower  state.  The  upper  electronic  state  is  slightly  split  (spin 
split)  into  two  levels  by  a  weak  interaction  of  the  electron  spin  with  the 
magnetic  field  generated  by  the  spinning  molecule  (and  other  effects,  Ref. 
Herzberg  1939).  Absorption  transitions  terminating  in  the  J  -  K  -  1/2  level 
are  designated  with  a  2,  while  those  terminating  in  the  J  »  K  +  1/2  level  are 
assigned  a  1.  A  similar  system  is  used  to  define  the  level  of  the  lower  elec¬ 
tronic  state  in  which  a  transition  originates.  The  lower  state  (2ir)  is  strongly 
split  by  spin-orbit  coupling.  The  separation  of  the  spin  1/2  and  spin  3/2 
levels  is  M.20  cm“^.  Transitions  originating  in  the  lower  energy  2"i/2  state 
are  assigned  the  value  1,  while  those  originating  in  the  higher  energy 
state  are  assigned  a  2. 

Thus  a  line  with  J  *  10.5  and  K  =  11  in  the  upper  electronic  state  and 
with  11.5  in  the  lower  spin  energy  level  of  the  lower  electronic  state 
(i.e.  nl/2)  would  be  labelled  P21  (11.5).  From  the  above  discussion  it  can 

be  seen  that'  in  general  an  absorption  transition  originating  from  a  given  J 
value  can  terminate  in  any  one  of  twelve  possible  systems,  since  there  are 
three  branches  and  four  possible  combinations  of  levels  for  each  branch. 
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Also  found  In  Part  Two  is  the  calculation  for  the  line  center  absorption 
coefficient,  KO(I).  The  variable  I  is  incremented  by  1  as  each  line  center 
absorption  coefficient  is  calculated.  (At  this  stage,  KO(I)  does  not  contain 
the  path  length  over  which  the  absorption  takes  place).  In  determining  this 
coefficient  the  three  subroutines  FUPPER,  FLOWER,  and  HONNUM  are  accessed  to 
determine  the  upper  state  energy,  lower  state  energy,  and  the  Honl-London 
factors,  respectively,  for  the  transition.  Also  at  this  time,  the  line  location 
option  is  queried.  If  the  theoretical  line  option  was  selected,  the  line 
center  transition  wavenumber  (WO(I))  will  be  given  by  the  difference  of  the 
upper  and  lower  state  energies  resulting  from  FUPPER  and  FLOWER.  If  the  experi¬ 
mental  line  location  option  was  requested,  then  the  wavenumber  read  with  the 
line  designation  is  used  for  WO(I) . 

D.I.b.3.  Main,  Part  3 


Part  Three  of  the  program  establishes  the  wavenumber  region  in  which  the 
transmission  calculations  are  done.  The  wavenumber  region  will  begin  2.0  Lorentz 
widths  and  2.4  Doppler  widths  (DELTAW)  before  the  line  center  wavenumber  of 
the  first  line  selected  and  will  end  an  equal  amount,  DELTAW,  after  the  line 
center  wavenumber  of  the  last  line  selected.  The  region  defined  by  these 
start  and  stop  wavenumbers  is  divided  into  499  panels.  The  transmission  is 
calculated  at  the  wavenumber  of  each  of  these  panels.  The  transmission  is 

obtained  for  a  panel  from  Beer's  Law, I  =  IQe  '^d.  The  absorption  coefficient 
kv  is  determined  by  summing  the  absorption  coefficients  at  the  panel  from  all 
spectral  lines  which  contribute  significantly  to  the  transmission  at  that 
panel's  wavenumber.  To  determine  if  the  absorption  at  a  panel  from  a  given  line 
is  significant,  its  contributrion  to  kv  is  divided  by  the  total  absorption  coef¬ 
ficient  determined  from  all  lines  previously  examined  at  that  panel.  If  this 
quotient  is  less  than  a  value  given  to  the  variable  ALLOW,  a  counter  (NHOLDER,  NHOLDL) 
is  indexed.  This  counter  is  reset  to  zero  if  the  quotient  becomes  greater 
than  allow.  When  the  counter  value  increases  to  three  (indicating  that  the 
Allow  condition  has  been  met  three  times  consecutively)  the  transmission  cal¬ 
culation  will  advance  to  the  next  panel.  The  consecutive  stipulation  reduces 
the  chances  of  exiting  the  coefficient  loop  when  encountering  an  isolated  weak 
line.  In  addition  to  the  ALLOW  check  for  the  absorption  contribution,  the 
loop  cannot  be  terminated  if  the  absorbing  line's  center  wavenumber  is  within 
DELTAW  from  the  panel  being  examined.  The  absorption  coefficient  must  frequently 
be  evaluated  at  some  distance  from  line  center.  For  the  Voight  profile  used, 
this  value  (kv^)  is  related  to  the  line  center  absorption  coefficient  kv^  by  the 
expression 


kvo  Re 


exp(-(u>i  +  ia ' ) ^)  erfc(Ao^  +  a') 


D1 


D-3 


4  V**/  ' 


where 


Re  =  Real  Part 


2  (vj  -  Vj)  r - — 

^  =  - 1 - li  yin  2 

Avd 

v°  =  line  center  wavenumber 


=  wavenumber 


AvD  =  Doppler  width  of  the  line 


_  AvL 


D2 


a'  =P  /ITT 

Avd 

Avj  =  Lorentz  width  of  the  line 

i  =/r 

erfc  =  complementary  error  function 
The  expression  for  is  evaluated  in  the  subroutine  WFUNC. 


When  the  transmission  calculation  has  been  completed  for  all  499  panels 
the  program  advances  to  Part  Pour. 

D. I. . b. 4^_Main,_Part  4^ 

The  function  of  Part  Four  of  the  main  program  is  to  record  the  inputs  of 
Part  One  and  each  panel  transmission  and  wavenumber  on  a  temporary  file  for 
later  use  by  another  program  (PLOT.).  The  entire  program  MAIN  will  be  repeated 
for  each  NO  concentration.  When  an  end-of-file  is  detected  in  the  input  data 
file  DATA-2,  the  main  program  will  terminate. 

D.I.c.  Input  File  OATA-2 

In  order  to  facilitate  changes  in  the  input  parameters  such  as  path- 
length,  gas  temperature,  and  NO  concentration,  these  values  are  placed  in  a 
data  file  named  DATA-2.  This  file  also  contains  the  spectral  line  information  and 
selections  for  the  various  options.  Since  DATA-2.  is  also  used  as  an  input  to 
a  resonance  absorption  program  (NO-SPECT.),  some  of  the  entries  in  DATA-2.  do 
not  apply  to  NO-ABS.  A  partial  list  of  DATA-2.  with  appropriate  definitions 
follows  in  Table  D-I. 

More  concentrations  can  be  added  to  the  end  of  DATA-2  if  desired.  The 
numbers  appearing  after  the  spectral  line  designation  are  the  center  wavenumber 
of  the  transition  as  reported  by  Deezsi  (1958).  Many  of  these  lines  were  not 
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resolved  by  the  instrument  used;  hence  several  different  lines  may  appear  at  an 
identical  wavenumber.  More  highly  resolved  values  for  the  line  center  wave- 
numbers  are  available  in  Engleman  et  al  (1970) .  Another  input  data  file  with 
Engleman's  experimental  line  locations  exists  in  file  DATAENG.  This  list 
does  not  contain  nearly  as  many  lines  as  DATA- 2 .  and  is  only  a  partial  tabulation 
of  all  the  lines  listed  in  Engleman  et  al.  Choosing  the  theoretical  line 
locations  (i.e.  placing  001  in  line  3  of  DATA-2.),  while  using  DATA-2.  will  very 
nearly  reproduce  all  of  the  experimental  values  available  in  Engleman,  et  al. 

A  partial  list  of  the  difference  between  Engleman's  experimentally  measured 
line  center  values  and  those  derived  theoretically  can  be  found  in  Fig.  Dl. 

The  differences  are  generally  less  than  the  room  temperature  Doppler  width  of 
NO  (.055  A  or  .0005  nm) .  At  high  rotational  values,  the  theoretical  locations 
for  the  ?22  lines  systematically  diverged  from  the  experimentally  measured  values. 
For  highest  accuracy,  lines  should  be  arranged  from  lowest  to  highest  energy 
(increasing  wavenumber) . 


D.I.d.  NO-ABS.  F UPPER 


The  purpose  of  this  subroutine  is  to  calculate  the  upper  state  energy 
of  the  Y(C,0)  transition.  The  equations  and  constants  used  are  the  same 
as  those  described  in  the  text  see  Section  III-B.6.d. 

D.I.e.  NO-ABS.  FLOWER 


This  subroutine  calculates  the  lower  state  energy  for  a  y(0,0)  transition. 
The  equations  and  constants  used  are  the  same  as  those  described  in  the  text 
see  Section  II-B.6.C, 


D.I.f.  NO-ABS .WFUNC 

Subroutine  WFUNC  calculates  the  off-line  center  fractional  absorption 
coefficient  (kv./k9)  for  a  particular  spectral  line.  This  is  accomplished  by 
determining  the  real  part  (Wl)  of  the  function. 


W(z)  =  exp  (-z2)  erf c  (-iz) 


where  i  =  yf- ! 

z  =  wj  +  ia'  (See  NO-ABS  Main,  Part  3). 

The  technique  for  evaluating  this  function  is  taken  from  Abramowitz  et  al.  NBS 
Handbook  of  Mathematical  Functions.  There  are  basically  two  methods  used  to 
determine  Wl.  If  either  ljj  or  a'  are  in  the  intervals 


0  <;  wj  s  3,9 

or 

0  <  a'  <  3.0, 
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THEORETICAL  (EUPPER-EL0WER)  VERSUS  EXPERIMENTAL  LINE  LOCATIONS 


(so-n 


the  value  W1  is  obtained  by  a  two-dimensional  linear  interpolation  from  the 
nearest  wj  and  a'  values  found  in  the  block  data  table  BLOCK-DATA  (see  Davis, 
et  al  (1976)  for  listing).  This  table  consists  of  a  31  by  40  by  2  array  whose 
elements  are  labeled  by  z(Al,  Bl,  Cl)  where  Al  corresponds  to  a’,  B1  corres- 
sponds  to  a) j ,  and  Cl  corresponds  to  the  real  or  imaginery  part  of  W(z)  i.e., 

W1  or  W2. 

The  element  z(ll,  3,  1)  is  the  real  part  of  the  function  u>(z)  evaluated 
for  z  «  0.2  +  1.0  i,  while  the  element  z(ll,  3,  2)  is  the  imaginery  part  of 
the  same  function  and  argument.  Conversion  of  the  arguments  of  z  to  a  table 
element  is  accomplished  by  taking  the  integer  portion  of  ten  times  the  argu¬ 
ment  value  plus  one.  This  procedure  yields  31  real  values  for  a',  and  40 
real  values  for  u>j  since  each  is  evaluated  in  increments  of  0.1  over  the 
regions  previously  defined.  For  example,  to  locate  the  real  part  of  w(z) 
for  z  =  0.2  +  1.0  i,  the  element  Al  is  found  by  multiplying  1.0  by  ten, 
adding  1  and  taking  the  integer  portion  of  this  sum.  Thus  Al  for  a'  equal  to 
1.0  is  11,  and  similarly  Bl  for  ok  equal  to  0.2  is  3.  Therefore  the  real 
part  of  W(z)  will  be  found  in  the  table  element  z(ll,  3,  1).  When  the  argu¬ 
ments  of  z  =  Wj  +  a'i  are  nonzero  beyond  the  first  decimal  place,  this  tech¬ 
nique  automatically  results  in  the  nearest  table  element  less  than  the  desired 
value.  The  nearest  table  element  greater  than  the  desired  value  is  obtained 
by  adding  one  to  the  values  Al  and  Bl  just  found.  From  the  table  elements 
found  at  these  element  locations,  the  two-dimension  linear  interpolation 
approximates  the  desired  value  of  W(z) . 

For  values  of  uij  and  a'greater  than  3.9  and  3.0  respectively,  polynomial 
approximations  for  W(z)  are  used.  For  coj  and  a’  in  the  intervals 

3.9  <  wj  <  6.0 
or 

3.0  <  a'  <  6.0 


W(z)  is  given  by 


W(z)  =  iz  (Tx  +T2  +T3)  +  e(z) 
where  Tj  «=  . 4613135/ (z2  -  .1901635) 

T2  -  .09999216/(z2  -  1.7844927) 
T3  *  ,002883894/(z2  -  5.5253437) 

and 

I e(z) |  <  2  x  10-6. 
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When  either  Uj  or  a'  is  greater  than  6.0,  W (?)  is  given  by 

W(z)  =  iz  (T1  +  T2)  +  n(z) 

where  T1  =  .5124242/(z2  -  .2752551) 

T2  =  ,05176526/(z2  -  2.724745) 

and 

| n(z) |  <  10 

The  last  part  of  this  subroutine  is  devoted  to  evaluating  W(z)  when 
either  or  both  arguments  of  z  are  negative.  In  practice,  the  evaluation  for 
negative  argument  is  not  needed.  Likewise,  for  this  application  all  the  ele¬ 
ments  of  the  block  data  table  labeled  z(Al,  Bl,  2)  are  never  used  since  these 
correspond  to  the  imaginary  part  of  W(z). 

D.I.g.  NO-ABS.  HONNUM 

Subroutine  HONNUM  calculates  the  normalized  Honl-London  factor  for  deter¬ 
mining  the  intensity  of  a  given  transition.  The  equations  are  normalized  such 
that  the  summation  of  the  Honl-London  factors  over  all  the  upper  state  J  values 
(J')  for  a  given  lower  state  J  value  (J")  equals  4(2J"+1).  In  general,  inten¬ 
sities  of  absorption  transitions  originating  in  the  J  =  1/2  rotational  level 
must  be  considered  separately  (Earls  (1935)).  These  intensities  are 

Q21  =  Qll  = 
r21  =  R11  =  2/3 

All  other  branches  are  identically  zero.  In  nitric  oxide,  the  spin-orbit 
coupling  constant  is  so  large  that  the  intensities  of  lines  connected  to 
the  J  =  1/2  level  are  correctly  predicted  to  at  least  eight  significant  digits 
by  the  general  expressions  given. 


D.II.  Program  Description  of  Resonant  Absorption  of  NO  Radiation 

The  program  NO-SPECT. predicts  the  fractional  transmission  of  radiation 
through  a  gas  mixture  containing  a  known  amount  of  nitric  oxide.  The  inci¬ 
dent  radiation  is  obtained  from  excited  nitric  oxide  molecules  at  low  pres¬ 
sure.  The  program  construction  and  description  of  NO-SPECT. is  very  similar 
to  the  program  NO-ABS. previously  discussed.  The  differences  are  contained 


mainly  in  defining  the  spectral  distribution  of  the  incident  radiation.  For 
NO-SPECT.  the  incident  intensity  is  obtained  from  an  experimental  measurement 
of  individual  spectral  lines  emitted  by  the  source  and  from  an  assumed  Doppler 
broadened  line  shape.  The  transmission  is  calculated  on  a  line  by  line  basis. 
Each  source  spectral  line  in  the  selected  wavelength  region  is  divided  into 
100  panels  starting  2.3  Doppler  widths  (DELTAW)  before  line  center  and  ending 
an  equal  amount  beyond  line  center.  The  incident  intensity  at  a  given  panel 
is  found  by  assuming  a  Gaussian  lineshape  of  known  full  width  at  half  maxi¬ 
mum  (FWHM)  centered  on  the  center  wavenumber  of  the  transition.  The  relative 
maximum  intensity  of  the  source  line  is  obtained  by  linear  interpolation  from 
a  table  of  experimentally  derived  intensities  that  have  been  normalized  by 
the  respective  Honl-London  factor  of  the  measured  line.  This  table  is  con¬ 
tained  in  the  data  file  DATA-2  in  lines  4,  5,  and  6.  See  discussion  of  FUNC 
for  more  information  on  source  line  strength.  The  absorbing  line's  strength  at 
each  panel  is  determined  in  a  manner  identical  to  that  described  in  the  discus¬ 
sion  of  NO-ABS . 

At  the  completion  of  NO-SPECT.,  the  integrated  intensity  of  each  lime 
leaving  the  absorbing  gas  is  recorded  on  a  temporary  data  file.  The  line 
intensities  are  recorded  with  each  line's  center  wavenumber  in  groups  for  each 
concentration  of  NO  listed  at  the  end  of  the  data  input  file  DATA-2.  This 
information  is  read  into  another  program,  PLOT,  that  convolves  these  lines 
with  a  spectrometer  slit  function  and  displays  the  results  graphically. 

An  additional  section  exists  at  the  end  of  NO-SPECT .MAIN  Part  Three  that 
outputs  the  integrated  fractional  transmission  of  very  closely  spaced  line 
groups.  This  option  cannot  be  accessed  by  the  input  data  file  DATA-2.  To 
obtain  output  from  this  option  the  mandatory  GO  TO  statement  (GO  TO  560)  must 
be  te.vioved  from  the  program.  These  cluster  transmissions  are  required  for 
determining  source  characteristics. 


All  subroutines  used  for  NO-SPECT.  are  identical  with  those  described  in 
NO-ABS.  An  additional  subroutine  (NO-SPECT.  FUNC)  is  required  for  source  line 
intensities  (see  D.II.a  below). 

Additional  variable  definitions  for  variables  found  only  in  NO-SPECT.  and 
a  figure  (Fig.  D-2)  further  describing  the  model's  operation  are  found  in  the 
variable  definition  section. 
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D.I.f.  NO-ABS  FUNC 


To  predict  the  amount  of  transmission  of  resonant  radiation  through  a  gas, 
the  relative  source  intensity  at  each  line  must  be  known.  Many  times  there 
are  overlapped  lines  or  lines  which  cannot  be  resolved  by  the  spectrometer  and 
whose  intensities  are  not  directly  measurable.  This  subroutine  linearly  inter¬ 
polates  all  line  intensities  from  a  least  squares  fit  of  normalized  resolved 
lines  whose  intensities  can  be  directly  measured.  For  electrical  discharge 
lamps,  a  log-linear  plot  of  the  measured  spectral  line  intensities  divided  by 
their  respective  Honl-London  factors  versus  the  line’s  upper  state  energy  will 
typically  result  in  a  curve  shown  in  Fig.  20.  For  our  purposes  this  curve 
was  fit  with  two  straight  lines.  The  intensity  of  any  spectral  line  can  now 
be  approximated  from  these  two  straight  lines  once  the  upper  state  energy  and 
the  Honl-London  factor  for  the  trasition  are  known. 


D.III.  List  of  Symbols 


MAIN  (NO-ABS) 

G(n)  -  Vibrational  energy  of  nth  vibrational  level 

QVIB  -  Vibrational  partition  function 

ALLOW  -  Ratio  of  absorption  coefficient  for  a  given  line  to  total 

absorption  coefficient  from  all  previous  lines  at  a  parti¬ 
cular  wavenumber 

LINST  -  Line  number  of  first  spectral  line  defining  start  of  interval 
to  be  examined 

LINEND  -  Line  number  of  last  spectral  line  to  be  examined 

NWOTHE  -  Conditional  branch  variable  to  select  theoretical lv-computed 
line  locations  (=1),  or  use  experimentally-determined  line 
locations  (=0) 

XDATA(I)  -  Wavenumber  array  of  resolved  upper  level  (“I)  source  spectral 
lines 

YDATA(I)  -  Normalized  logarithmic  intensity  array  of  resolved  source 
spectral  lines 

NP  -  Number  of  source  lines  used  to  characterize  source  output 

HEADS  -  Version  label 

AP  -  Broadening  parameter  (a') 

TA  -  Absorber  temperature 

TE  -  Emitter  temperatura  (applies  only  to  NO-SPECT.) 

EL  -  Path  length  of  absorption 

YHGT  -  Records  option  taken  for  variable  NWOTHE 

BVEQ0  -  Spectroscopic  constant  for  the  zero  vibrational  state 

AFL  -  Spin  orbit  coupling  constant  for  lower  level 

2 

C1S12  -  f  number  for  transitions  connected  to  the  n,  /„  level 


C1S32 


-  f  number  for  transitions  connected  to  the  tt 


level 


3/2 

C2  -  hc/kTA 

C0NST1  - 2/Jii(2) 

C0NST2  -  Doppler  width  of  source  divided  by  line  wavenumber 
(2/2(in2>  kTE/msC2) 

C0NST3  -  Doppler  width  of  absorber  divided  by  line  wavenumber 
(/2Un2)  kTA/msC2) 

C0NST4  -  2.3  •  (C0NST2) 

J  -  Spectral  line  position  in  array  of  spectral  lines  input,  (not 

the  rotational  level) 

RJPP  -  Rotational  level  (truncated  integer) 

JPP  -  Rotational  level  (truncated  integer) 

BRANCH  -  Rotational  branch  (P,Q,R) 

ICODE  -  Rotational  branch  (=1,  P,  branch,  =  2  Q  branch,  '  3  R  branch) 
NUP  -  Spin  split  state  in  upper  level  (2E ) 

NLO  -  Spin  state  in  lower  level  (  it) 

WO(J)  -  Line  center  wavenumber  of  line  number  J 

ISPIN(J)  -  Array  of  lower  spin  states  for  all  spectral  lines  in  input  list 

S(J)  -  Honl-London  factor  of  the  Jth  line 

NLINES  -  Total  number  of  spectral  lines  read  from  input  list 

vfjo  -  Nitric  oxide  concentration  (molecules/cmJ) 

KO(I)  -  Line  center  absorption  coefficient  of  Ith  lines 

NINT  -  Wavenumber  range  of  interest  is  divided  into  this  many  panels 

(max.  499) 

DWJ  -  Doppler  width  of  source 


DWL 


-  Doppler  width  of  absorber 

DELTAW  -  2  Lorentz  widths  +  2. A  Doppler  widths  absorber  (cm-^) 

WST  -  Wavenumber  at  beginning  interval  of  interest 

WND  -  Wavenumber  at  ending  interval  of  interest 

AZJ  -  Lower  wavenumber  limit  of  mandatory  computation  for 

given  spectral  line 

BZJ  -  Upper  wavenumber  limit  of  mandatory  computation  for  a  given 

spectral  line 

DEL  -  Incremental  wavenumber  of  panels  in  region  of  interest 

W(I)  -  Center  wavenumber  of  Ith  panel 

E(l)  -  Source  intensity  of  Ith  panel  (constant) 

JLEFT  -  $  when  examining  lines  of  larger  wavenumber  than  present 

panel-  =  1  when  examining  smaller 

SUML  -  Total  absorption  coefficient  from  all  lines  of  wavenumber 

greater  than,  or  equal  to,  nearest  spectral  line 

SUMR  -  Total  absorption  coefficient  from  all  lines  of  wavenumber 

less  than,  or  equal  to,  nearest  spectral  line 

NHOLDR  -  Number  of  consecutive  lines  whose  contribution  to  the  total 
absorption  coefficient  is  less  than  minimum  required 

NHOLDL  -  Number  of  consecutive  lines  whose  contribution  to  the  total 
absorption  coefficient  is  less  than  minimum  required 

NLESR  -  Number  of  consecutive  lines  whose  contribution  to  the  total 
absorption  coefficient  is  less  than  minimum  required 

NLESL  -  Number  of  consecutive  lines  whose  contribution  to  the  total 
absorption  coefficient  is  less  tb_n  minimum  required 

BETA  -  2/in2  (vj-Vj^/boppler  width  of  absorber) 

WFUNC  -  Complex  complementary  error  function  to  determine  absorption 

coefficient  at  arbitrary  distance  from  given  spectral  line 
(fraction  of  line  center  absorption  coefficient) 
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TERM 

TNEW 

II 

CAY (I) 

DAT ( I ) 

TAUTOP 

TAUBOT 

TAU 

ALPHA 

W02(I) 

TJ2(I) 

FLIPPER 

IT 

J 

M 

T 

G 

WE 

WEXE 

BE 

ALPHAE 

DV 


Absorption  coefficient  at  arbitrary  distance  from  given  spectral 
line  (fraction  of  line  center  absorption  coefficient) 

Absorption  coefficient  at  arbitrary  distance  from  given 
spectral  line 

Rotational  level 

Total  absorption  coefficient  for  the  Ith  panel  from  all 
spectral  lines  considered 

Intensity  of  Ith  panel  after  passing  thru  absorbing  gas 
Total  transmitted  intensity 
Total  incident  intensity 

Ratio  of  transmitted  intensity  incident  intensity 

Ratio  of  absorbed  intensity  to  incident  intensity  ( 

Reordered  panel  wavenumber,  largest  to  smallest 
Reordered  panel  transmissions  (see  W02(I)) 


Vibrational  level 
Rotational  level 

2 

Spin  split  level  in  upper  state  (I  level) 
Term  value  of  electronic  level 
Term  value  of  vibrational  level 


w  x 
e  e 


BV  -  Bv 

GAMA  -  y 

FN  -  Term  value  of  rotational  level 

F  -  Term  value  of  upper  state  (^E  level) 

FLOWER 

V  -  Vibrational  level 

J  -  Rotational  level 

M  -  Lower  spin  state  (^rr  level) 

USS  -  Spin  split  level  of  upper  state  (“I  level) 

T  -  Term  value  of  electronic  level 

G  -  Term  value  of  vibrational  level 

DV  -  Dv 

A  -  Lover  state  spin-splitting  constant  level) 

BE  -  Be 

BV  -  Bv 

WE  -  W 

e 

WEXE  -  <*>exe 

IB  -  Rotational  branch 

NSIGN  -  Direction  of  lambda  doubling  shift 

LMDADB  -  Magnitude  of  lambda  doubling  term  in  lower  state  ( 

F  -  Term  value  of  lower  level  (n) 

WFUNC 

XF  -  Same  as  BETA  (MAIN . NO-ABS) 


n  level) 
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YF  -  a'  (Broadening  parameter) 

W1  -  Interpolated  real  part  of  the  complex  function  W(z)  where  z  =  XF  +  YF  i 

WZ  -  Interpolated  imaginary  part  of  the  complex  function  W(z) 

Z(DUM1,DUM2,  1)  -  Real  part  (W(DUM2  +  (DUMl)i)) 

Z (DUM1 , DUM2 ,  2)  -  Imaginary  part  (W(DUM2  +  (DUMl)i)) 


X(DUM) 

-  Real  Fart  of  W(z) 

(Table  Value) 

Y (DUM) 

-  Imaginary  part  of 

W(z)  (Table  Value) 

I  FLAG 

-  Conditional  branch 
<f>  ->.when  ■+  XF  >  YF 

variable  whose  value  is 
>  4> 

1  when  XF  <  <(.,  YF  >  4> 

4  when  XF  <  (fi ,  YF  <  <(> 

II  -  Converts  YF  into  an  integer  -  1  (Lower  bound  in  table  look-up 

for  YF) 

JJ  -  Converts  XF  into  an  integer  -  1  (Lower  bound  in  table  look-up 

for  XF) 

III  -  II+l  upper  bound  in  table  look-up  for  YF  interpolation 

•JJ1  -  JJ+1  upper  bound  in  table  look-up  for  YF  interpolation 

Z1  -  Linearly  interpolated  value  for  XF  with  Y=CONST=II 

Z2  -  Linearly  interpolated  value  for  XF  with  Y=C0NST=II1 

W(l)  -  Interpolated  real  value  of  W(XF+YF  i) 

W(2)  -  Interpolated  imaginary  value  of  W(XF+YF  i) 

ZC  -  z  =  XF+  YFi  (complex)  =  argument  of  W  function 

HONNUM 

NUP  -  Upper  vibrational  state  (v') 

NUPP  -  Lower  vibrational  state  (v") 

MEGAN  -  Upper  spin  state 
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MEGAM 


-  Lower  spin  state 


IB 

J 

A 

BE 

ALPHAE 

BV 

RNUM 

DENOM 

HONL 

MAIN 

XDATA(I) 

YDATA(I) 

TE 

CONST2 

CONST4 
EO(  J) 
NINT 

W(I) 

e(i) 

nisi 


-  Rotational  branch 

-  Lower  rotational  level 

-  Lower  state  spin-splitting  constant 


-  Numerator  of  Honl-London  factor 

-  Denominator  of  Honl-London  factor 

-  Honl-London  factor  for  transition 
(NO-SPECT) 

-  Lower  state  wavenumber  of  source  spectral  line 

-  Natural  logarithm  of  resolved  source  spectral  line 
intensity  normalized  by  the  Honl-London  factor  of 
the  transition 

-  Source  Temperature 


-  Doppler  width  of  source  divided  by  line  wavenumber 
(2v/2(£n2)k(TE)"/mc2) 

-2.3  times  C0NST2 

-  Line  intensity  of  Jth  source  spectral  line 

-  Number  of  calculations  (panels)  performed  over  a  single 
source  spectral  line 

-  Wavenumber  of  Ith  pmel  of  a  given  source  spectral  line 

-  Intensity  of  Ith  panel  of  n  given  source  spectral  line 
(Gaussian  shape  assumed) 

-  ‘<.6  Source  Doppler  widths  FWHM 
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TAUTOP 


-  Integrated  intensity  of  lines  passing  through  the 
absorbing  gas 

TAUBOT  -  Integrated  intensity  of  lines  incident  on  absorbing  gas 

FUNC 

X  -  Wavenumber  of  upper  level  transition  of  line  of  interest 

XDATA  -  Wavenumber  array  of  resolved  upper  level  transitions  for  source 

YDATA  -  Natural  logarithm  array  of  normalized  resolved  source 
line  intensities 

NP  -  Number  of  points  used  to  characterize  source  line  intensities 

NP1  -  NP  minus  1 

FUNC  -  Linearly  interpolated  normalized  natural  logarithm  of  the 

intensity  of  unresolved  source  line 
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D.IV.  Program  Description  PL0T.PL0T5 


D.IV.a. Introduction 


The  program  PL0T.PL0T5  convolves  the  spectral  transmission  calculated  in 
programs  NO-SPECT.  and  NO-ABS.  with  a  spectrometer  slit  function  to  aid  in 
reducing  experimental  data.  Results  of  the  convolution  are  presented  graphi¬ 
cally. 

The  program  is  divided  into  four  parts.  A  discussion  of  the  function 
of  each  part  is  presented  in  the  following  section. 

D. IV. b.  Discussion 


D.IV.b.l.  Part  One 


Part  One  of  the  program  reads  the  transmitted  intensity  and  wavenumber 
at  each  wavelength  present  In  the  radiation  source.  For  resonant  radiation 
(NO-SPECT.),  these  wavenumbers  are  the  center  wavenumbers  of  each  line. 

The  transmissions  at  these  wavenumbers  are  the  integrated  intensities  of 
individual  lines  (i.e.,  the  transmission  verses  wavenumber  is  nonzero  only 
at  the  center  wavenumber  of  the  transition;  no  direct  information  on  line 
shape  is  contained  in  these  transmissions).  For  absorption  from  a  continuum 
source  (NO-ABS.),  the  wavenumbers  read  are  uniformly  spaced  throughout  the 
region  of  interest.  In  this  case,  the  original  lineshape  of  the  absorbing 
line  is  retained.  Original  input  parameters  for  NO-SPECT.  or  NO-ABS.  are 
read  for  documentation  purposes  on  the  plots.  Options  are  presented  for 
plotting  the  data  as  read  (zero  slit  plots)  in  addition  to  the  convolution 
plots  generated  by  PL0T.PL0T5.  At  this  time,  the  user  is  also  queried  for 
the  slit  width  of  the  convolution.  The  slit  width  is  defined  as  the  apparent 
full  width  at  half  maximum  transmission  (FWHM)  of  a  resolved  line  in  the 
experimental  spectra.  The  slit  width  must  have  units  of  angstroms  (X).  Part 
One  continues  with  a  conversion  of  the  original  wavenumbers  to  wavelengths. 
Since  line  order  (increasing  wavelength)  is  assumed  in  the  construction  of 
later  parts  of  the  program,  Part  One  contains  a  sorting  loop  (ISORT)  which 
looks  twenty  lines  ahead  of  the  present  line  location  to  resort  lines 
and  intensities  for  continually  increasing  wavelength. 

D-IVj.b^2i.parI  Two 

Part  Two  produces  plots  of  the  data  as  received  after  the  line  sort 
operation  described  in  Part  One.  A  vertical  line  is  drawn  at  the  location 
of  each  wavelength  initially  read.  The  relative  heights  of  these  vertical 
lines  are  proportional  to  the  spectral  intensity  read  at  each  wavelength. 
These  plots  are  referred  to  as  zero  slit  plots,  since  no  instrument  broaden¬ 
ing  has  been  included.  It  should  be  reemphasized  that  plots  for  resonant 
absorption  (NO-SPECT.)  do  not  retain  line  shape  information.  Zero  slit 


plots  for  resonant  absorption,  and  absorption  of  continuum  radiation  are 
shown  in  Figs.  D3  and  DA,  respectively.  Part  Two  also  sets  the  convolution 
interval  (VIDE)  for  Part  Three. 

D. IV. b. 3. Part  Three 


Part  Three  computes  the  resultant  intensity  verses  wavelength  when  the 
individual  intensities  are  examined  by  a  measuring  intrument  (spectrometer) 
of  finite  resolution.  The  transfer  function  of  the  spectrometer  can  be 
determined  by  scanning  and  examining  the  resultant  profile  of  a  spectral  line 
whose  bandwidth  is  very  narrow  relative  to  the  spectrometer's  bandpass.  For 
the  spectrometer  used  in  our  measurements,  the  transfer  function  can  be 
approximated  as  a  triangular  function  for  bandpassesQgreater  than  0.5  A  and 
by  a  Gaussian  function  for  bandpasses  less  than  0.5  A  (Figs.  D5  and  T6) .  To 
produce  the  convolved  spectral  plots,  a  wavelength  region  starting  twice  the 
convolution  interval  (2'WIDE)  before  the  wavelength  of  the  first  spectral 
line  and  ending  an  equal  amount  after  the  last  spectral  line  is  divided  into 
2499  panels.  At  each  panel,  a  convolution  interval  is  defined  (WIDE)  whose 
width  is  twice  the  spectrometer  slit  function  for  the  triangular  case  and 
six  times  the  slit  function  for  the  Gaussian  case.  This  convolution  interval 
is  centered  on  the  panel  at  which  the  calculation  is  being  performed.  Beyond 
this  interval,  the  contribution  to  the  convolution  is  considered  to  be  zero. 

A  loop  is  set  up  to  determine  the  first  spectral  line  number  (ISAVE)  and 
total  number  of  spectral  lines  (NOPTS)  within  the  convolution  interval. 
Starting  at  spectral  line  number  ISAVE,  the  intensity  of  each  of  the  next 
NOPTS  lines  is  weighted  with  a  Gaussian  or  triangular  function  according  to 
the  relation. 


a 

TRIADD  =  HT(JBI.)  exp  (-  ( (XLAM- ANG(.JBL) ) /DLAM  •  1.66512]“) 


or 


TRIADD  =  HT(JBL)  [DLAM-  [  xLAM-ANG (JBL) j  ]  /DLAM 

where  TRIADD  is  the  weighted  intensity,  HT(JBL)  is  the  original  line  inten¬ 
sity,  DLAM  is  the  FWHM  of  the  spectrometer  in  angstroms,  XLAM  is  the  panel 
wavelength  at  which  the  convolution  is  taking  place  and  ANG(JBL)  is  the 
wavelength  of  the  spectral  line  being  weighted.  Figure  D2  presents  symbol 
definitions  graphically.  Weighted  intensities  at  the  edge  of  the  convolu¬ 
tion  interval  for  the  Gaussian  case  are  down  by  more  than  10“  ^  relative 
to  the  same  intensity  at  the  center  of  the  interval.  The  weighted  intensities 
over  the  entire  convolution  interval  are  summed  under  the  variable  name 
SUMMER.  When  the  sum  is  completed,  the  value  of  SUMMER  is  placed  in  the 
array  PLOTY  (KKOUNT)  and  the  center  wavelength  of  the  convolution  interval 
is  placed  in  the  array  1’LOTX  (KKOUNT)  where  KKOUNT  is  the  panel  number  at 


ZERO  SLIT  WIDTH  PLOT  FOR  RESONANT  LINE  SOURCE 


FIG  D4 
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which  the  calculation  was  performed.  The  weighting  and  summing  process  is 
repeated  for  all  2499  panels.  The  maximum  value  in  PLOTY  (KKOUNT)  is 
assigned  to  the  variable  SPMAX  and  the  first  value  of  SPMAX  (no  nitric 
oxide)  is  assigned  to  the  variable  DlV.  The  variable  DIV  is  divided  into  the 
intensities  of  each  panel  for  each  concentration  to  normalize  the  maximum 
peak  height  to  unity  for  plotting  purposes.  The  remaining  sections  of  Part 
Three  write  the  values  of  PLOTX  and  PLOTY  into  a  temporary  storage  file  and 
reread  this  file  for  display  purposes  when  all  spectra  for  each  nitric  oxide 
concentration  have  been  calculated.  This  file  is  also  used  for  wavelength 
scale  expansions;  no  new  calculations  are  done  when  a  scale  expansion  is 
requested.  Sample  convolved  outputs  for  a  continuum  and  a  resonant  absorption 
case  are  shown  in  Figs.  D7  and  D8.  Numerical  values  for  features  in  the  con¬ 
volved  spectra  can  be  obtained  utilizing  Part  Four  of  PL0T.PL0T5. 

D.  IVib_14._Part  Four 

Part  Four  produces  the  numerical  value  of  extrema  on  each  concentration 
curve  for  a  given  wavelength  range.  This  option  is  accessed  by  a  carriage 
return  after  the  completion  of  the  convolution  plots.  The  wavelength  of 
the  start  of  the  extrema  interval  is  assigned  to  the  variable  XSTR  and  the 
end  to  XSTP.  The  minimum  and  maximum  values  and  their  wavelengths  within 
this  interval  are  determined.  The  extrema  of  the  first  case  (no  nitric 
oxide)  are  used  to  normalize  the  extrema  of  following  cases.  This  results  in 
the  fractional  value  of  the  extrema  for  each  curve  in  the  specified  interval. 
The  extrema,  fractional  extrema,  extrema  wavelengths,  and  concentrations  are 
presented  followed  by  a  request  to  determine  if  another  region  is  to  be 
examined.  A  response  of  '  NO'  (include  blank  before  word  NO)  will  terminate 
the  program.  A  sample  of  this  option  is  included  in  Fig.  D9. 
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PLOT. PLOT  5 


D.IV.d. 

NCASE 
I START 

DAT  El 

TIME1 

STMAX 

XMIN 

XMAX 

DLAM2 

IPLT1I 

PLTMIN 

PLTMAX 

YHGT 

NLINES 

ENO 

TE 

TA 

EL 

APRIME 

TAU 

HEADS 


List  of  Symbols 

counter  for  concentration  case 

flag  for  no  NO  concentration  (equals  1  for  NO  =  0.^ 

0  for  NO  #  0. 

date  of  run 
time  of  run 

maximum  line  intensity 

minimum  wavelength  to  be  plotted 

maximum  wavelength  to  be  plotted 

slit  function  of  spectrometer  in  angstroms  (FWHM) 

option  for  plots,  1  for  zero  slit  plots 

0  to  delete  zero  slit  plots 

same  as  L1NST 
same  as  LINEND 

indicates  method  of  obtaining  line  locations  for  data 
(YHGT  equals  1  for  theoretically  computed  line  locations 
and  0  for  experimentally  measured  line  locations) 

total  number  of  spectral  locations  read 

NO  concentration  (molecules/cm  ) 

Source  temperature  (Kelvin) 

absorber  temperature  (Kelvin) 

path  length  of  absorption  (cm) 

gas  broadening  constant 

integrated  fractional  transmission 

program  version  code 
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I PLOT 1 


same  as  IPLT1I 


DLAM 

LINST 

LINEND 

LINT 

ANG(I) 

HT(I) 

151 

152 
WIDE 

II 

XPLT(l) 
XPLT(2) 
YPLT(l) 
YPLT(2) 
I  SAVE 

KKOUNT 

DELLAM 

PP1 


same  as  DLAM2 

starting  spectral  line  number 

final  spectral  line  number 

total  number  of  spectral  lines  read 

wavenumber  of  Ith  spectral  line  (later  converted  to 
wave 1 eng th) 

intensity  of  Ith  spectral  line 

line  number  of  start  of  wavelength  sort 

line  number  of  end  of  wavelength  sort  (IS1  +  20) 

width  in  angstroms  of  spectral  region  to  be  convolved 
at  a  given  wavelength 

same  as  LINT 

line  center  wavelength 

same  as  XPLT(l) 

0 

integrated  intensity  of  line  at  XPLT(2) 

line  number  of  first  spectral  line  in  convolution 
interval 

panel  counter  (1-2499) 

total  wavelength  interval  divided  by  2499 

wavelength  of  first  panel 


l£ 


*-32 


SPMAX 


maximum  convolved  intensity  for  a  given  concentration 


PI 

P2 

XLAM 

I 

NOPTS 

IPN1 

JBL 

SUMMER 

TRIADD 

PLOTX(I) 

PLOTY(I) 

DIV 

TRNCAS  (NCASE.I) 

DECIS 

XSTR 

XSTP 

J1 


wavelength  of  start  of  individual  convolution  interval 

wavelength  of  end  of  individual  convolution  interval 

center  wavelength  of  convolution  interval 

counter  to  determine  which  spectral  lines  lie  within  the 
convolution  interval 

same  as  KOUNT 

line  number  of  last  spectral  line  in  convolution  interval 

counter  for  lines  within  convolution  interval 

integrated  intensity  of  all  lines  within  the  convolution 
interval 

line  intensity  weighted  by  slit  function 
wavelength  of  convolved  line  transmissions  (Ith  panel) 
intensity  of  convolved  line  transmissions  (Ith  panel) 
maximum  convolved  intensity  for  no  NO  present 

array  of  intensities  for  each  concentration  (NCASE)  at  Ith  panel 

option  for  extrema  in  wavelength  region  specified 

starting  wavelength  for  extrema  option 

ending  wavelength  for  extrema  option 

counter  for  concentration  case 


ST LOOK 


XSTR  minus  one  percent  of  extrema  interval 


SMALL 

BIG 

DIFl 

DIF2 

LAMIN 

LAJIAX 

MINSAV 

MAXSAV 

MINRAT 


minimum  transmission  over  extrema  interval 
maximum  transmission  over  extrema  interval 
minima  check 
maxima  check 

wavelength  of  minimum  transmission 
wavelength  of  maximum  transmission 
minimum  transmission  with  no  NO  present 
maximum  transmission  with  no  NO  present 
minimum  fractional  transmission 


MAXRAT 


maximum  fractional  transmission 
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